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THE DISTRIBUTION OF Fet*+ AND Mg** IN COEXISTING 
OLIVINES AND PYROXENES' 


HANS RAMBERG AND GEORGE DEVORE 
University of Chicago 


ABSTRACT 


THE JOURNAL OF GEOLOGY 


May 1951 


The Mg/(Mg + Fe) ratios of coexisting olivines and orthopyroxenes from different types of rocks have 
been determined by means of indices-of-refraction measurements. Some data from literature have also been 


added. 


The distribution function of Mg and Fe in olivine-pyroxene pairs is discussed. It is shown that this func- 
tion may be used as a geologic thermometer when some additional thermochemical data are obtained. The 
distribution function also indicates that orthopyroxene is a nonideal mixed crystal, the intermediate mem- 
bers of the series being more stable than one would expect. This nonideal character of the orthopyroxene 
series is indicates by deviation from linearity of the indices of refraction and the unit-cell dimensions when 
these quantities are plotted against the mol fraction of enstatite in the pyroxenes. The @ and } parameters 
of the unit cell of different orthopyroxenes are measured. The aluminum content in orthopyroxenes is dis- 


cussed briefly. 


INTRODUCTION 


The stability and equilibrium condi- 
tions of the minerals are most directly 
revealed by syntheses under controlled 
physicochemical conditions of the min- 
erals in question. For metamorphic min- 
erals this method meets severe practical 
difficulties because of the sluggishness of 
silicate reactions below the solidus tem- 
perature, and it is likely that we shall 
have to wait a long period of years before 
the complex equilibrium conditions of 
the rock-making minerals are explored 
thoroughly. In the meantime, field ob- 
servations combined with theoretical 
thermodynamics may provide informa- 
tion about the equilibrium relations of 
metamorphic minerals. The present 
study of the distribution of Fe+*+ and 
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Mg** in coexisting olivines and orthopy- 
roxenes represents this latter type of 
approach. 


COMPOSITION OF COEXISTING PYROXENE 
AND OLIVINE 


Table 1 and figure 2 give the composi- 
tion of coexisting olivines and orthopy- 
roxenes found in natural rocks, which, 
with two exceptions, are quartz-free. In 
figure 1 are compiled coexisting olivines 
and calcium-free clinopyroxenes from ex- 
perimental work of Bowen and Schairer 
(1935, figs. 8, 9, 10, 11, 12, 13, 14). Be- 
cause we are interested only in the dis- 
tribution of Mg and Fe between coexist- 
ing rhombic pyroxene and olivine, inde- 
pendent of other phases which may be 
present, the complete Bowen and 
Schairer diagram is not reproduced. 

There are some striking differences 


TABLE 1* 
COEXISTING OLIVINES AND ORTHOPYROXENES IN ROCKS 


Rererence or Locatity 


and G. W. DeVore from the Laramie Ra‘ 


once. Material from the literature is included only if it was definite that the two minerals were in the same 


I 
(2V= 


° 
52°54") 
. 730 
.726 
73° 
.726 
-732 
-735 
-737 
.728 
-735 
730 
-744 
.728 
-735 
.726 
-735 
-745 


Sapphirine rock 
Pyroxenite 
Amphibolite sill 
Amphibolite sill 
Olivine diopside 
pyrope rock 
Dunite 
Diabase dike 
Gabbro picrite 


Picrite 


Anorthosite 
Gabbro 


Anorthosite 
Lamprophyre 
Anorthosite 
Anorthosite 
Norite 

Basic pegmatite 
Anorthosite 
Ti-Fe ore 
Lamprophyre 
Norite 
Lamprophyre 
Anorthosite 
Lamprophyre 
Lamprophyre 
Syenite 
Fe-Ti ore 
Dark lava 
Dark lava 
Syenite 
Norite 

Fe-Ti ore 
Lamprophyre 
Ti-Fe ore 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 
Syenite 


Ramberg (1948, p. 17) 
Greenland 

Greenland 

Greenland 

Eskola (1921, p. 22) 


Greenland 

Wyoming 

Wager and Deer (1930, 
159) 

Wager and Deer (19309, 
p. 68) 


Wyoming 
Wager and Deer (1939, 
p. 68) 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Barth (1930, p. 10) 
Barth (1939, p. 10) 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming 
Wyoming . 
Wyoming 


* Greenland specimens were collected by ¥. my ; Wyoming specimens were collected by W. H. Newhouse, A. F Hagner 
. From t 


1 mie Range complex, many duplications of the above listed determina: 
tions were observed. In the mineral pair Xo = 0.65-0.63 with Xp» = 0.65-0.67 there are twelve additional examples. In the pair 
Xo = 0.61-0.62 with Xp = 0.66-0.67 there are seven additional examples. Most of the other determinations are duplicated at least 


nd specimen and if the 


analysis or index of refraction was listed. All composition values were determined from the index of refraction curves in Kennedy 


(1047, pp. $64 and 567). 
t +0.004. 
talpha + 0.004. 


Mou 
FRACTION OLIVINE PyYROXENE 
+ 0.002 +0.002 
None.........--]| 0.85 1.684 1.682 
-94 .78 1.668 1.693 
.78 1.668 1.693 
| 
.go -93 1.672 1.671 
.88 -79 1.676 1.690 
i 66 -72 1.725 1.703 
62 63 .736 1.715 
q 61 .67 - 737 1.707 
.60 64 - 738 1.713 
6-80-Ha........ .56 63 -748 1.715 
IMIs5-148...... .55 .67 1.707 
Tr11D-7B...... 54 .66 -756 1.709 
51 59 . 760 1.719 
12-37-Kb...... 42 .59 -774 1.719 
Tro-23A...... 63 .778 1.715 
41 -64 .778 1.713 
5-24-F... 41 .778 
q 12-37-KD ... -37 54 . 786 
j 6-25-378....... .36 48 
} .35 51 .792t 
6-25-30B .590 .797t 
12-37-Ke..... 24 
5-24-Fa......... .19 49 .826t 
12-37-Kd....... 19 43 
5-24-Fd...... 16 -47 .795¢ 
5-25-A.... 14 -798f 
32-37-G...... 13 .40 .800f 
4-30-D.........| .10 | .47 
6-44-Cl....... .806t 
.07 .812t 
| 0.07 0.39 
104 
\ 


' 
p™ 


4 


Mol % 40 
Fic. 1.—Compositions of artificial olivines coexisting with Ca-free clinopyroxenes (from the work of 


Bowen and Schairer, 1935). The break between the curves is caused by the existence of a melt, with the 
minerals on the left half of the diagram and with only solid phases on the right half. The dashed lines, 


p-p’, p’’-p’”", o-0’, and o’’-0""’, represent compositional changes with temperature if free silica is absent. 


i 
80 6 Mol % 40 
Fic. 2.—-Compositions of natural olivines coexisting with orthopyroxenes of table 1, plotted so that the 


Mg Siicates 
mean of each coexisting pair falls on the line AB. Line AB forms an angle of 45° with the base of the diagram. 
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between the Mg-Fe distribution in the 
natural olivine-orthopyroxene pairs and 
the artificial pairs in the phase diagram 
of Bowen and Schairer. The artificial oli- 
vine is always somewhat poorer in Mg 
than the coexisting pyroxene, with a dif- 
ference between the mol fractions of Mg 
silicate (denoted as X, for olivine and X, 
for orthopyroxene) of not more than 1o- 
15 per cent. In the natural assemblages 
the olivines also contain less Mg-silicate 
than does the coexisting orthopyroxene 
as long as the average mol fraction, 
(X, + X,)/2 is less than about 0.65. The 
difference between X, and X, of the 
natural coexisting minerals is, on the 
average, considerably greater than in the 
artificial system within the more Fe-rich 
portion of the diagram. Within the Mg- 
rich part of the diagram the situation is 
reversed, and in some cases the natural 
olivine is even richer in Mg than the co- 
existing natural pyroxene (fig. 2 and 
table 1). 


THERMODYNAMICS OF THE PYROXENE- 
OLIVINE ASSEMBLAGES 


There are several possible reasons for 
the difference between the Mg-Fe dis- 
tributions in natural and in artificial oli- 
vine-pyroxene pairs. The artificial min- 
eral pairs are formed at higher tempera- 
ture and lower pressure than most of the 
natural assemblages. In the artificial sys- 
tem it is a question of coexistence of oli- 
vine and Ca-free clinopyroxene, whereas 
the natural pyroxenes here considered 
are all orthorhombic. Bowen and 
Schairer (1935) have concluded that the 
difference in Mg/(Fe + Mg) ratio be- 
tween the rhombic and monoclinic Ca- 
free pyroxene series is negligible. At con- 
stant pressure (P?) and temperature (7) 
one should, therefore, expect to find the 
same Mg-Fe distribution between olivine 
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and clinopyroxene as between olivine 
and orthopyroxene. Natural olivine-py- 
roxene assemblages occur together with 
several minerals not present in the arti- 
ficial system. This fact, however, will not 
disturb the Mg-Fe distribution between 
olivine and pyroxene, inasmuch as 
neither of the minerals tends to dissolve 
appreciable amounts of “‘foreign’’ ions; 
and therefore the presence of other min- 
erals does not affect appreciably the 
thermodynamic properties (free energy) 
of either olivine or pyroxene. Thus the 
difference of the Mg-Fe distribution in 
the artificial and in the natural olivine- 
pyroxene pairs should be mainly due to 
temperature and/or pressure differences, 
although some of the natural coexisting 
olivine-orthopyroxene pairs may cer- 
tainly not represent equilibrium. 

To explain the interaction between oli- 
vine and pyroxene in quartz-free 'rocks, 
consider olivine 0’ and pyroxene p” in 
figure 1 and assume that every trace of 
free silica or tridymite is removed and 
that the temperature is lowered. The 
composition of olivine and pyroxene 
cannot follow the olivine and pyroxene 
curves, respectively, because to do so 
would require a reaction between olivine 
and silica to form pyroxene. The only 
possible reaction, if any, will be that a 
certain amount of Mg*+ goes from the 
one mineral to the other at the same time 
that an equal amount of Fe** goes in the 
opposite direction. The dashed curves 0” 
and p’’—p’” may, for example, repre- 
sent the changing composition of olivine 
and coexisting pyroxene when free silica 
is absent. : 

The curves and p’’-p’” hold 
only for a certain ratio of amount of 
pyroxene and olivine present together. 
In case the olivine and pyroxene are pres- 
ent in any other relative amounts, the 
change of composition with temperature 


| 
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Fe*+ AND Mg** IN COEXISTING OLIVINES AND PYROXENES 


will follow curves somewhat different 
from and 

In a quartz-free assemblage beneath 
the olivine (or pyroxene) curve in figure 
1, the composition of olivine (or pyrox- 
ene) can change at a given temperature 
from pure forsterite (or enstatite) to the 
composition determined by the olivine 
(or pyroxene) curve at this given tem- 
perature. If both mineral series are ideal 
mixtures, there is a simple relation be- 
tween the compositions of pyroxenes and 
olivines which are in equilibrium at a 
given temperature and pressure. In our 
following discussion we shall consider 
the olivine and pyroxene series to be ideal 
mixed crystals, as Sahama and Torgeson 
(1949) concluded, and under this pre- 
liminary assumption determine the char- 
acter of the distribution constant and its 
relation to temperature and pressure. 

Consider first the reaction: 


Mg,SiO, + 2 FeSiO,; 
forsterite ferrosilite 


Fe,SiO, + 2Mgsio, . 
fayalite ) 


enstatite 


The standard free energy of this reaction 
at a given temperature, 7” K. is: 


AF’ = Figy + 2Fen — Fior — 2Fiere » (2) 


where Fj,,, etc., are the molal free 
energies of fayalite, enstatite, etc., in 
their pure state at the temperature and 
pressure in question. The olivine, 
(Mg,Fe),.SiO,, and the pyroxene, 
(Mg,Fe)SiO,, are mixed crystals; there- 
fore, the actual reaction cannot be so 
simple as that expressed in equation (1). 
The fayalite molecules that form will 
dissolve in the forsterite present, and 
the enstatite molecules formed will dis- 
solve in the original ferrosilite. By this 
homogeneous mixing the partial molal 
free energy of forsterite in the olivine will 
decrease, whereas the partial molal free 
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energy of fayalite in olivine will increase. 
Similarly, the partial molal free energies 
of enstatite and ferrosilite in the pyrox- 
ene will increase and decrease, respec- 
tively, when reaction (1) goes to the 
right. 

Conditions for chemical equilibrium 
between olivine and pyroxene at a con- 
stant temperature 7”° K. are that the 
sum of the partial molal free energies of 
forsterite in olivine and of ferrosilite in 
pyroxene is equal to the sum of the par- 


‘ tial molal free energies of fayalite in oli- 


vine and of enstatite in pyroxene (eq. [3]). 
The two minerals are known to be ionic 
crystals in which the exchangeable con- 
stituents Fe++ and Mg** occur in the 
same state in both minerals, namely, as 
individual ions. Under these conditions 
one may consider the reactions connected 
with one mol of the exchangeable ions: 


+ FeSiO;— FeSi,/sO2 + MgSiO; . 


Equilibrium between olivine and pyrox- 
ene at 7” K. exists when the following 
relation is satisfied : 


where F/,,, etc., are now the partial 
molal free energies of fayalite in olivine, 
enstatite in pyroxene, etc., at 7’ K. The 
general relationship between the partial 
molal free energy, F, and the composition 
of an ideal mixed crystal is 


F,=F,+RTWX;, 


where F ; is the free energy of « .mponent 
i in its pure state. Hence equation (3) can 
be expressed as follows: 


+RT’ In x; —RT’ \n X’ (4) 
— 


torr In (1 X}) =O, | 


198 


where R is the gas constant, X, is the mol 
fraction of MgSi,,,0, in olivine, X; is the 
mol fraction of MgSiO, in pyroxene, 
1—X°% is the mol fraction of FeSi,,,0, in 
olivine, and 1 —X} is the mol fraction of 
FeSiO, in pyroxene.’ 

Rearrangement of equation (4) gives: 


+ Freer 3Fiay — Fen 
RT’ 


=e (6) 


where AF’ is the standard free energy of 
reaction (1). 

Because AF’ is constant at a constant 
temperature (and pressure), it follows 
that the term 


is constant at constant P and 7, provided 
that the olivine and the pyroxene are 
both ideal mixed crystals. Hence, by 
knowing the X, and X, for one pair of 
stably coexisting olivine-pyroxene assem- 
blages at a given temperature, the dis- 
tribution constant, 


X, 1-X,_ 


K, 


is known and will give corresponding 
compositions of any stable pair of co- 
existing olivine-pyroxene assemblages at 
that temperature, provided, of course, 
that the two mineral series are ideal mix- 
tures. 

Applying the considerations above at 


? It is here assumed for simplicity that the par- 
tial molal free energy of SiO, in the lattice is con- 
stant at varying Mg/(Mg,Fe) ratio in the olivine 
and pyroxene. 
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any temperature, 7, different from 7”, we 
find that the following relation is valid: 


where 
xX, 1-X, 
Zs 
is the distribution constant at any tem- 
perature, 7; 


=K 


i-xX 
P 


is the distribution constant at a reference 
temperature 7”; and AF’ and AF are the 
standard free energies at 7” and 7, re- 
spectively. Now, assuming that the en- 
tropies of enstatite, forsterite, ferrosilite, 
and fayalite are constant within the tem- 
perature range 7’ —T, the standard free 
energy difference can be expressed as 
follows: 


AF’ — AF = AS (T’ —T), (8) 


where the entropy difference is 
AS = Stor + 2 Storr (Stay + 2S.n) 
This can be introduced into equation (7): 


i-X, X. (9) 


K 


If the pyroxenes and olivines are ideal 
mixed crystals and AS were known, these 
minerals would provide a usable geologic 
thermometer. By using X and X, from 
the Bowen and Schairer experiments, we 
can get experimental values of the con- 
stant 


x, 


P 


il or 
, , x 1 
F _X, 1-X, 
1— Xx?" 
| 
\ 


at known temperatures. Any stable natu- 
ral assemblage of olivine and pyroxene 
would give a new distribution constant, 


K 


which will enable us to measure the cor- 
responding temperature at which the 
natural assemblage was ‘‘quenched.”’ 

However, from the data gathered in 
table 1 and plotted in figure 2, it can be 
concluded that the rhombic pyroxenes 
and/or the olivines are not ideal mixed 
crystals because K is not constant at con- 
stant T. The rocks compiled in table 1 
and in figure 2, in all likelihood, were 
formed at different temperatures because 
they are partly lavas, partly plutonic ig- 
neous rocks, and partly metasomatic and 
metamorphic rocks. There is little reason 
to believe, though, that there is a pro- 
nounced relation between the tempera- 
ture which the mineral represents and 
the average 


X 


in the olivine-pyroxene pairs in such a 
way, for example, that a high X,, repre- 
sents high temperature or vice versa. So 
far, we have reason to believe, because of 
similar facies (granulite) existing across 
the diagram (fig. 2), that the Mg-rich 
minerals from Greenland are formed at a 
temperature similar to the temperature 
represented by the iron-rich minerals 
from the Laramie Range, Wyoming (ta- 
ble 1). We therefore think it is true that 
the temperature interval inside which 
the rocks (in table 1) have formed (or 
recrystallized) is independent of the aver- 
age X,,, of the coexisting olivine-pyroxene 
pairs. On the other hand, we see that, sta- 
tistically, the X,/(1—X,) (1—X.)/X. 
values of our rocks change regularly with 
the average X,,, of the assemblages, the 


Fe** AND Mg** IN COEXISTING OLIVINES AND PYROXENES 199 


iron-rich assemblages having much 
greater K-values than the Mg-rich as- 
semblages (table 2). The reason must be 
that the minerals concerned are not ideal 
mixed crystals, for which equations (7) 
and (9) should be valid. For this reason 
equations (7) and (9g) have to be re- 
stricted to a certain narrow (X, + X,)/2 
range when used as a temperature indi- 
cator. That is, only olivine-pyroxene 
pairs with a rather similar mean X should 
be used when the A’s are compared in an 
attempt to determine the temperature of 
“quenching.” 


TABLE 2 


DISTRIBUTION CONSTANTS FROM 
NATURAL ROCKS 


Specimen No. Xe. Xp | Xm Kx 

3835 0.906 | 0 8&3 0°. 89 ©. 205 
12-56-C . 88 -79 83 0.402 
4087 .68 .7° .69 1.099 
5-37-D .60 63 61 1.135 
12-37-Kb. .. .42 59 1.987 
122 39 1.626 
6.25~36b. 32 45 3.059 
12-37-Ke. .. 24 54 30 3.716 
12-37-G 13 .40 26 4.400 
6-44-Cl ©.09 0.54 0.31 11.870 


Before discussing the matter further, 
we shall consider the effect of pressure 
upon the distribution constant, K. At a 
constant temperature, 7’, and a given 
pressure, P’, we have (eq. [5]) 
AF’ 
DRT’ 
Keeping 7’ constant but changing pres- 
sure to P will change AF’ to AF and 
hence also K’ to K,: 


In K’ = — (10) 


AF 


(11) 


InK,= — 


Equation (11) minus equation (10) gives 


Kk’ 
p=K e( 
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where AF’ — AF is the standard free en- 
ergy difference caused by the pressure 
difference P’—P. Therefore, 


AF’ —AF=AV(P’-P), 
‘where AV is the mol volume difference: 
Viayt 2 ( Vitor + 2 Vterr) . 


The effect of pressure upon the distribu- 
tion constant, therefore, is 


K = K’ 


We know approximate values of AV from 
the following values of densities, d, and 
mol volumes V: 
Forsterite (Poldervaart, 1947): 
d = 3.19 g/cc; V = 44.10 ce. 
Fayalite (Poldervaart, 1947): 
= 4.14 g/cc; V = 49.23 cc. 
Enstatite (Poldervaart, 1947): 
d = 3.18 g/cc; V = 31.57 cc. 
Ferrosilite (Birch, Schairer, and Spicer, 1942): 
d = 3.97 g/cc; V = 33.23 ce. 
AV = 1.81 cc. 

Although the density data are not en- 
tirely reliable and we do not know the 
compressibilities of the minerals, we shall 
see whether high pressure can have an 
important effect on the A-values. As- 
sume that P = 10,000 atm. (more than 
30-km. depth) and P’ =o atm. The 
value of 7” is, for example, 800° K. 
Therefore, 


K= K’e 10,000 /2°82. 1° 800, 


where R = 82.1 atm. cc/degree. There- 
fore, 
K =0.87K’, 

which means that the pressure sensitivity 
of the equilibrium constant is rather neg- 
ligible, assuming that nothing serious 
happens to the mol volume under the in- 
creasing pressure. It is therefore likely 
that the main reason for the discrepancy 
between the A-values in the natural oli- 
vine-pyroxene assemblages and _ those 
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from the Bowen and Schairer experi- 
ments is due primarily to temperature 
differences. 

Although it appears that the complete 
olivine or pyroxene mixed crystal series 
cannot be ideal mixtures, one may rea- 
sonably assume that the pyroxenes be- 
have like an ideal mixture within the 
rather narrow range X, = 0.2-0.5, and 
olivine within X, = 0.05-0.30. Within 
these compositions, equation (g) should 
apply fairly well and is thus a useful 
geologic thermometer. However, the im- 
portant quantity, AS, is not known with 
the accuracy necessary to calculate the 
temperature of the natural assemblages. 
For the sake of order of magnitude, an 
attempt to estimate AS is made from the 
entropy data compiled in table 3. The 
resulting AS is 1.0 cal/mol/degree at 
1,200° K. and A is 0.5 cal/ mol/ degree at 
298° K. These values are certainly sub- 
ject to rather great probable errors and 
do not indicate much more than that AS 
is close to zero or perhaps has a small 
positive value of the order of magnitude 
1-3 cal/ mol, degree. 

There exists another possibility of cal- 
culating AS. From figure 1 the corre- 
sponding X, and X, are experimentally 
known over a temperature range. Pro- 
vided that the data compiled from 
Bowen and Schairer in figure 1 are suffi- 
ciently accurate, the different K’s which 
can be computed from corresponding X » 
and X,,at different temperatures will give 
AS within this range of temperature. 
Table 4 gives corresponding X,’s, X.’s, 
K’s, and the computed AS at the differ- 
ent temperatures. It is interesting to note 
that the AS thus calculated has a positive 
value of a magnitude rather comparable 
with that of the AS calculated from the 
experimental data in table 3. 

In figure 3, K and In K are plotted 
versus temperature. One could now plot 


yy 

q 

3 

; 

q 
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TABLE 3 
DATA AVAILABL 


Enrropies tin Cat/ Deorer/ Mou 


At 208°1 K. | Between 298°1 K. and 1,200° K. Total at 1,200° K. 


.50+0. 36 
.g8 +0. 30 
.95 +0. 10 
.72+0.41 
.27+0.35 
.25 
59 est|| 
73.61 est(q); 73.28 est(t) 
72.95 +(q); 72.95 +(t) 
93.1t 
60.5 +est(q); 60.5 est(t) 
| 62.03 test(q); 62.09 est(t) 
63.6+est(q); 63.68 est(t) 
FeSiO, (av. est)... ‘ | 62.04+(q); 62.09 +(t) 


* From Kelley (1041, 1943)- 

t Storsgg 1,200 = + Senstage 1,200° 
t*‘est(q)’’ means that the entropy of quartz is used. 
§‘‘est(t)’’ means that the entropy of tridymite is used. 


Stet z9:1-200 = Senstoyg 1.200 ~ 


= S802 996-1, 200° 
** Stere = Stay —S FeO. 
tt Sterre = 4( Stay + Ssio,) 


Sterr = Ssio, + Sreo 


tributi t K 


4 16 18 2.0 22 24 26 
Distribution Constant K 


Fic. 3.—Diagram of temperature versus the distribution constants A and the natural logarithms of K 
from the composition interval of X_ = 0.05-0.30 and Xp = 0.20-0.50 of the artificial system of Bowen and 
Schairer (1935). Data are listed in table 4. 


| 
| 
2 
1300 
|200 
1000 ° ° 
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the K’s from the natural olivine-pyrox- 
ene pairs (table 2), with (XY, + X.)/2 
values comparable with those of the syn- 
thetic minerals in table 4, on the exten- 
sion of the K versus temperature curve in 
figure 3. Nos. 12-37-Ke, 12-37-G, 5-24- 
Fd,and 6-44-Cl in table 2 have the proper 
(X, + X.)/2 values. It may be seen that 
their A-values—3.716, 4.46, 4.655, and 
11.87, respectively—will correspond to 
very low temperatures indeed, if the K or 
In A versus temperature functions are 
linear. 
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X, (and X,) in the natural rocks. In par- 
ticular, when X, increases, K decreases 
strikingly, as is shown in table 2 

We have seen that there is every 
reason to believe that the tempera- 
ture (and pressure) represented by the 
large number of rocks here compiled does 
not change in any regular manner with 
the mean X. Therefore, the regular 
change of K with mean X must be due 
to deviation from ideality of one or both 
of the two mixed crystals. 

Sahama and Torgeson (1949) have 


TABLE 4 
DISTRIBUTION CONSTANTS AND AS FROM THE BOWEN AND SCHAIRER DIAGRAM 


AS/2 
(Cal/ Degree/ Mol) 


Ink 


©. 205 


©. 3075 + 
2.43 


+2.22 
+2.30 
+1.25 
+0.66 


+1.20 


-4955 
-4947 
- 
.6678 


.6678 


0.7419 


There is little reason to believe, how- 
ever, that either K or its logarithm is a 
linear function of te.nperature. Further- 
more, the coexisting Y, and X, values 
found in the diagrams of Bowen and 
Schairer may not be accurate enough to 
allow a correct calculation of K as a tem- 
perature function. Therefore, the presen- 
tations here should be regarded as only 
qualitative, inasmuch as too many data 
are lacking at present to apply the oli- 
vine-pyroxene assemblages as an abso- 
lute geothermometer. 


DEVIATION FROM IDEALITY OF THE 
ENSTATITE-FERROSILITE 
MIXED CRYSTAL 


Let us now return to the interesting 
feature that K is changing regularly with 


demonstrated that the olivine series is 
nearly ideal because it shows no heat of 
mixing. These authors also investigated 
the pyroxene series but, unfortunately, 
measured only compositions close to the 
end members. Although they concluded 
(1949) that the orthopyroxenes appear 
to represent an ideal mixed-crystal series, 
their results are not conclusive. The pres- 
ent workers did not have the opportu- 
nity to undertake thermochemical meas- 
urements of the intermediate pyroxenes, 
but some additional facts indicating that 
the orthopyroxenes are not ideal mixed 
crystals will be emphasized. 

The Fe-Mg distributions indicate that, 
as is usually suspected and was also 
demonstrated by Sahama and Torgeson, 
the iron-rich members of the pyroxenes 


| 1,577 0.464 0.381 1.36 
| 1,323 .426 1.50 
£,900.......] 3,473 395 .285 | 1.63 
1,100... | 1,373 .320 -195 1.95 
1,050... 1,323 . 260 -153 1.95 
q 
| 1,273 0.110 2.10 
. 


3 


Refractive index ny 


s 
21: L69 
£ 
a 


i i i A. iL 


MgSi0z 80 60 40 
Mol % 


Fic. 4.—Indices of refraction plotted against composition of orthopyroxenes, showing deviation from 
linearity. The straight line is drawn for comparison. 
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are considerably less stable as related to 
iron-rich olivines than are the magne- 
sium-rich pyroxenes as compared with 
magnesium olivines. Therefore, when the 
average Mg/(Mg + Fe) ratio is less than 
a certain value (some 0.6-0.5), olivine is 
able to absorb much more iron than py- 
roxene; at larger average Mg (Mg + Fe) 
values, pyroxene will take as much iron 
as olivine or more. These relations indi- 
cate that pyroxenes of intermediate com- 
position should have a relatively greater 
lattice energy or less molal free energy 
than that which corresponds to the same 
mol fractions in ideal mixtures. 

One should expect that a relatively 
great stability in certain portions of a 
mixed-crystal series would reflect itself 
qualitively on the mol volume, so that 
there would be a tendency toward con- 
traction and hence a negative deviation 
from Vegard’s law, which states that the 
mol volume (or unit cell) is a linear func- 


tion of the compositions. Again, the re- 
fractive indices may indicate such a con- 
traction so that they should show a posi- 
tive deviation from linearity. The indices 
of refractions of the clinoenstatite-clino- 
ferrosilite series, measured by Bowen 


and Schairer (1935, p. 199), plotted 
against weight percentage are straight 
lines. This means that plots of n, and n, 
versus mol fraction will be curves with 
the convex side up. The same character- 
istics are shown by the n, and n, indices 
of the rhombic pyroxenes (fig. 4). It is 
well known, on the other hand, that the 
ideal olivine mixed crystals have indices 
of refraction which vary linearly with 
mol fraction. 

To check the indication of negative 
deviation of the orthopyroxenes from 
Vegard’s law as suggested by the indices 
of refraction, X-ray measurements of a 
series of orthopyroxenes were made. A 
5.7-cm.-radius Phillips powder camera 
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and FeA, radiation were used. The a and 
b parameters were measured by using the 
good reflections (060), (650), (12, 0, 0), 
and (14, 5,0). The (14, 5, 0) reflection is 
strong, with an angle @ varying from 
65°55 to 67°40. The parameters a and } 
were measured by (060), (650), and 
(12, 0, o), Corrections for film shrinkage 
and rod thickness were made.’ The Bragg 
angles of the (060) and (12, 0, 0), as well 
as the (650) reflections, are approxi- 
mately the same; thus these reflections 
should be subjected to the same syste- 
matic errors due to the rod thickness and 
uneven film shrinkage. It is therefore 
assumed that the ratio between a and 6 
measured from (060), (12, 0, 0), and 
(650) is nearly correct, even if the abso- 
lute values of a and b may be less accu- 
rate. By means of this ratio, |a (12, 0, 0)]}, 
|b (060)], the values of a and b were re- 
calculated so as to fit the good back- 
reflection (14, 5,0). The value of Fe,, is 
taken as 1.93507 A. The maximum error 
of the a and 6 recalculated on the basis 
of (14, 5, ©) is assumed to be 0.002 A. 
Table 5 and figure 5 show the results. 
The chemical compositions were first 
determined by the index of refraction, 
but the large scattering of the X-ray 
data soon indicated that index of refrac- 
tion was insufficient to show the chemical 
composition of the pyroxenes investi- 
gated by X-ray. We, therefore, subse- 
quently analyzed those of our samples 
which could be prepared for reliable 


5 The half-circumference, C/2, of the shrunken 
film was measured as the distance between the 
middle point of several corresponding front reflec- 
tions and the middle point of several corresponding 
back reflections. The film was supposed to have 
shrunk evenly. The distances Y on the film corre- 
sponding to the double deviation angles, 46, for 
the several! reflections were measured, and a cor- 
rection, AL = r(1 + cos 20), was subtracted. The 
symbol ¢ is the radius of the rod. The corrected devi- 
ation angle was then taken as 

26= degrees . 


| 
q 
| 


TABLE 5 


@ from 6 from 
(060) | (14, 5, 0)* | (a4, 5, 0)* 
(A) | (A) | 


.225 
. 256 
279 
.279 
. 292 
.317 
18.408 


Mol % 
Fic. 5.—The a and b parameters of the unit cell of the orthopyroxenes. Data are listed in table 5. Speci- 
men numbers are given at the bottom of the figure. 


UNIT CELL DIMENSIONS OF RHOMBIC PYROXENES 
| @ from 4 
No. xX | (12, ©, 0) 
(Ay | 
37303-........| | 28.032 | 8.832 | 8.828 
Pa eK .720 18.278 | 8.850 | | 8.852 
37006.........| | 28.247 | 8.816 | | 8.821 
| «18.278 | 8.875 | | 8.876 
6-82-B. ...| «656 | 18.286 | 8.857 | 8.847 
37850 .§90 | 18.283 | 8.898 | 8.900 
-467 | 18.313 8.922 8.924 
| 18.305 | 9.0a7t | 9.033 
t Mean 6 from (060) and (650). 
18.42 
° 
304 / 18.38 
° 
900 B34 
od ° 
od 
£ 
2 & = 
$896 BO 
w 
2 oP 
20 
B26 
? 
888 
° 
8.84 
BS 
80 50 40 is 
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TABLE 6 


ALUMINUM CONTENT OF CHEMICALLY ANA- 
LYZED ORTHOPYROXENES CALCULATED AS 
ATOMIC PER CENT Al OF ALL METALS. AND 
THE DISTRIBUTION OF ALUMINUM BETWEEN 
THE SILICON AND MG, FE POSITIONS OF THE 


Per Cent Al 


in Si 
Position 


Atomic 


I. 
4- 
I. 
4. 
4. 
2. 
6. 
2. 
6. 
5. 
4. 
6. 
°. 
I. 
3. 
6. 
8. 
° 
6. 
4 
3. 
2. 
2. 
1. 
1. 
I. 


* Nos. 1-21 are analyses 1—21 in Poldervaart (1947, p. 168); 

A, no. 37393, this paper; B, Bushveld (Hess and Phillips, 1938, 
p. 453, anal. no. 3); C, Labrador (Washington and Merwin, 1923, 
p. 64, anal. no. 2); b, India (Srirama Rao, 1947, p. 138); £ 
no. 37096, Greenland, this paper; F, Heinola, Finland (Lokka, 
1943, p. 29); G, no. 344, Greenland, this paper; H, Heinola, Fin- 
| kka, 1043, 29); 1, Adirondacks (Buddington, 1939, 
. 32); J, no. 6-82-B, Laramie Range, Wyoming, this paper; 
K, ganda (Groves, 1935, p. 152); L, Minnesota (Bayley, es, 
. 5); M, Turku, Finland (Hietanen, 1047, P. 1037, NO. 1059); 
& India (Rajagopalan, 1946, p. 317); O, nia’ Rajagopal . 
1946, p. 315); P, Terra del Fuego (Lokka, 1943, p. 32); R, no. 
37218, Greenland, this paper; and S, no. x Z, Greenland, this 


paper. 


chemical analyses. Unfortunately, some 
of the already X-rayed samples could not 
be analyzed because of some weathering. 
The noticeable decrease in the size of the 
b edge of the unit cell given by points 
37096 and 6-82-B is accompanied by a 
high aluminum content in these pyrox- 
enes. Aluminum makes up 6.44 per cent 
of the total metals in 37096 and 6.07 per 
cent in 6-82-B (see fig. 5 and tables 6 
and 7). Even with the high content of 
aluminum in the mineral, the index of 
refraction falls on the curve established 
by aluminum-poor or -free orthopy- 
roxenes. 

The cell contraction along the b edge 
is probably due to the small Al ion sub- 
stituting for the larger Mg**+ and Fet* 
ions. The simultaneous exchange of Si by 
Al apparently does not produce a large 
enough expansion to compensate for the 
contraction. The variation of the a edge 
of the cell is more puzzling. The Al-rich 
pyroxene no. 37096 shows a definite con- 
traction of the edge a, whereas no. 344, 
which also is rather rich in Al (tables 6 
and 7), has a relatively large a. These ir- 
regularities may be explained by the data 
in table 6, which show that aluminum 
does not distribute itself evenly between 
the Si and the (Mg,Fe) positions in the 
pyroxenes. Therefore, pyroxenes with 
the same percentage of aluminum may 
have very different fractions of the alu- 
minum content in the Si and the (Mg, 
Fe) positions. Because aluminum in the 
(Mg,Fe) positions causes a contraction, 
and aluminum in the Si positions most 
likely causes an expansion, of the unit 
cell, it is natural that the unit-cell di- 
mensions should show some irregularities 
as long as the actual positions of the alu- 
minum content of the pyroxenes are not 
determined. 

It is significant that, even when only 
the aluminum-poor orthopyroxenes are 
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MINERAL 
No.* Xp Per Cent 
| | | 
gas | . 896 75 100.0 
ces 888 75 62.5 
.817 68 53.8 
814 57 94.2 
q 809 31 53-9 
781 74 69.8 
-758 39 100.0 
-756 03 56.9 
D.. .| 66 13.6 
.740 36 50.0 
44 46.2 
_ | .722 61 80.9 
-720 13 44.1 
H.. pe 71.2 
‘ 87 85.0 
O81 52 56.6 
.664 71 77-7 
614 23 0.0 
. 580 95-3 
12 42.8 
. 500 19 40.5 
J -453 19 9.1 
.238 86 73-9 
434 82 30.4 
“ . 166 30 100.0 
a Exe) 0.0 
o.118 | 12 57-4 
| 
\ 
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considered, the dimensions of the unit 
cell (or at least the parameters 6 and a) 
versus composition are still curved lines 
which one should expect if the orthopy- 
roxenes are not ideal mixtures between 
Mg and Fe silicates. 

It is also interesting to note that the 
aluminum content in orthopyroxenes 
seems to depend upon the Mg/ Fe ratio, 
as is seen from table 6 and figure 6. A re- 
gion of composition (X = 0.8-0.4) in- 


IN COEXISTING OLIVINES AND PYROXENES 


TABLE 7 
ANALYSES OF RHOMBIC 
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field below the curve 75 per cent of the 
time, assuming normal distribution of 
data within the class intervals. The prob- 
ability of additional data following the 
same distribution is probably higher, 
since there should be some definite lat- 
tice limitations of the orthopyroxenes to 
the acceptance of aluminum. It is no- 
ticed that nine analyses are plotted in 
figure 6 that are not included in the sta- 
tistical treatment. These analyses were 


Specimen No. 


6-82-B* 


* Analysis by H. k Wiik. 
No. 37303: 


Yo. 37096: (Cas oFe ~~ le. 
No. 344: oo Mng Ors. 


t wee by B. Bruun. Na,O, K.O, and P.O, n.d. Minerals calculated to the standard formula 


No. XYZ: (Cas. 


cludes the most aluminum-rich ortho- 
pyroxenes. In order to get a statistical 
expression for the validity of the sample 
and the possibility that new data would 
fall outside present data, the curve on 
figure 6 representing the mean plus the 
tolerance factor (Eisenhart ef al., 1947, 
p. 102) times the standard deviation of 
each class interval is drawn. The toler- 
ance factor is such that at least 95 per 
cent of any additional data will fall in the 


uM 


199) (Sip. 


rejected because they do not recalculate 
to the standard formula. 

We thus think that our data with ref- 
erence to the distribution of aluminum 
in the orthopyroxene series are represent- 
ative. Aluminum is mainly found in py- 
roxenes with X ranging from about 0.8 
to 0.4. However, this does not necessarily 
mean that the pyroxene lattice within 
this range of X is particularly susceptible 
to accepting aluminum. It may also 


| 
37393* | 37096* | | 37218t | XYZt 
SiO,. .. 56.56 | 51.45 | 52.64 50.06 | 49.63 | 47.51 
0.03 0.16 69 | 0.28 
ALO,.. ..| ©@.43 5.98 3.82 5.57 2.16 1.20 
Fe(,...... ©. 32 1.52 | 0.16 0.88 | o.58 1.07 j } 
5.33 15.12 | 16.85 | 19.30 | 30.69 | 44.03 
MnO... 0.08 | °. 26 0.22 0.32 | 0.30 | 0.05 
MgO 36.80 24.56 24.54 21.42 | 15.22 3.990 
0.16 0.00 0.98 1.68 | 1.16 | 0.92 4 
H,0+ 0.31 ©. 34 ©.29 0.45 | nd. n.d. 
HO-. 0.04 0.12 0.05 | 0.05 nd. | n.d. 
Total... .| 100.06 | 99.51 | 99.71 | 100.42 | 100.11 | 99.95 
©.9109 | 0.725 | 0.720 0.656 | 0.466 | 0.135 
m,t0.001...| 1.668 1.702 | 1.701 | 1.711 1.734 | 1.770 ' 
i 
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mean that the rocks in which Mg-rich 
or Fe-rich pyroxenes develop have little 
aluminum and thus the chance for find- 
ing aluminum in such pyroxenes is small. 
On the other hand, if pyroxenes in the 
interval X = 0.8-o0.4 are especially sus- 
ceptible to accepting aluminum in the 
lattice, then this is only another indica- 
tion that the orthorhombic pyroxenes do 
not constitute an ideal mixed crystal 
series. 
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SUMMARY 

The distribution of Fe and Mg in 
natural pyroxene-olivine pairs is investi- 
gated and compared with the distribu- 
tion of the same elements in artificial 
olivine-pyroxene pairs. It is found that 
the difference in the distribution con- 
stant is mainly due to different tempera- 
tures of formation of the rocks, these be- 
ing at lower temperatures than the syn- 
thetic samples. It is shown that accu- 
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© Specimen Value 


e imen not used in the 


Qu 
MgSi03 80 FeSi0z 


Mol % 
Fic. 6.—Atomic per cent of aluminum in orthopyroxenes plotted against the Mg /(Mg + Fe) ratio. Statis- 


tical test for the validity of the distribution is represented by the curve. Aluminum data are listed in table 6, 
statistical data are listed in table 8. 
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rately measured entropy data of the sys- 
tem would enable us to calculate tem- 
peratures of formation from the distribu- 
tion of Mg and Fe in olivine-pyroxene 
pairs in rocks. 

The orthorhombic pyroxenes do not 
represent an ideal mixed-crystal series. 
This is shown by the Fe-Mg distribution 
as well as by the index of refraction and 
the dimension of the unit cell across the 
series. 
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siderably. The effect of Al upon the a 
parameter is not clear because some Al- 
bearing pyroxenes do have normal or 
even large a, while others show contrac- 
tion. This may be due to the fact that Al 
can distribute itself in any proportion be- 
tween the Si and the Mg,Fe positions in 
the lattice. 
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TABLE 8 


Class | No. of 
Interval* | Samples | 
(C.1.) 


Standard 
Deviation 


N FIGURE 6 


Tolerance | Sample 
Factor 


1.0 -0.81 
©.90-0.71 
0. 80-0. 61 
©.70-0. 51 
©.60-0.41 
©.50-0. 31 
©. 40-0. 21 
©.30-0.11 
©.20-0.01 


* Class interval of 0.19 includes all samples that have Mg/(Mg + Fe) ratios within this com- 


position range. 


t The tolerance factor is such that, with a risk of 0.75, at least 95 per cent of additional data will 
be included within the limit curve. Values of the tolerance factors were taken from Eisenhart, 


Hastay, and Wallis (1947, p. 102). 


The aluminum content in pyroxenes is 
discussed briefly. It is found that the 
chance of finding aluminum in rhombic 
pyroxenes is greater for pyroxenes of in- 
termediate to Mg-rich composition (X = 
0.8-0.4). Al does not noticeably affect 
the index of refraction but makes the } 
parameter in the unit cell contract con- 
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SIZE VERSUS SORTING IN SOME CARIBBEAN SEDIMENTS! 


JOHN C. GRIFFITHS 
Pennsylvania State College 


ABSTRACT 


The relationship between average size (phi median) and sorting (phi percentile deviation) in sediments is 
explored empirically by means of statistical analysis of data from over twelve hundred individual samples 
taken from Miocene and Oligocene horizons in southern Trinidad, the Eocene Scotland series from Barbados, 
and the Pliocene Superficial Gravels of British and Dutch Guiatia. The association between size and sorting 
is estimated by means of the correlation coefficient, and the relationship is reduced to linear equations of the 
form PDg = bMdg + a. Deviations from the trends are measured in terms of the standard error of estimate 
(Sy). On the basis of the statistics it is shown that, in water-deposited sediments, size varies with sorting 
over the range o¢ (1.00 mm.) to 6¢ (0.0156 mm.) in median diameter and o-4.5¢@ in phi percentile deviation. 
Best-sorted sediments (PDg < 1@) possess median diameters around 2.5¢ (0.177 mm.), and sorting becomes 
poorer with both increasing and decreasing size from the fine-sand rarge (2-3). Departures from the 
trend of over 1 Sy unit indicate that in finer-grained sediments (clays and fine silts) flocculation had taken 
place during grain-size analysis; whereas in coarser-grained sediments (very fine sands and coarse silts) the 


departures imply events of geological significance which interfered with norma! sedimentation. 


INTRODUCTION 


The treatment of measured data re- 
sulting from investigations into the 
petrography of sediments has been sys- 
tematized by Krumbein (19344, 6, 1936a, 
b, 1938, 1939, 1940; Krumbein and 
Aberdeen, 1937; Krumbein and Petti- 
john, 1938); the principal emphasis was 
perforce devoted in the early stages to 
the reduction of grain-size distributions 
to a few representative statistics. Subse- 
quently, Krumbein illustrated how such 
statistics could be further used in inter- 
relating size with other measured char- 
acteristics of grains (Krumbein and 
Pettijohn, 1938) and in following changes 
in size within a depositional unit (Krum- 
bein and Aberdeen, 1937). One of the 
aims of Krumbein’s work was to study 
the interrelationships of the statistics 
among themselves (1940) and to eluci- 
date this interrelationship in terms of 
sedimentation. It is important to realize 
that if, as the earlier pioneers of sedi- 
mentary petrography believed, the varia- 
tion of the several parameters—size, 
shape, etc.—in sediments is infinite, 


* Manuscript received August 7, 1950. 


then the problem of studying measured 
characteristics of sedimentary rocks be- 
comes insuperable. Krynine, however, 
approached the problems of sedimentary 
petrography on the implied, if not stated, 
basis that, although it is theoretically 
possible for this infinite variation to oc- 
cur in, for example, textural characters 
and mineral composition, in practice it 
appears that there are very definite limi- 
tations and that within single rock 
groups there is a clear tendency for a pre- 
ferred texture to be associated with a 
preferred composition (1943). If Kry- 
nine’s approach is correct, then the solu- 
tion to most problems of sedimentary 
petrography not only becomes possible of 
achievement but also can be accom- 
plished with a relatively few representa- 
tive samples when once the interrela- 
tionships between the significant vari- 
ables of sedimentary rocks have been 
established. The present investigation is 
aimed at the elucidation of the interrela- 
tionship between size and sorting in 
sedimentary rocks. 

A relationship between size and sort- 
ing of a sediment may be deduced on 
a priori grounds from the behavior of a 
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transporting medium. Indeed, as current 
velocity ultimately determines both ca- 
pacity and competence, it will also exer- 
cise a definite size selection on the ma- 
terial available for transportation. The 
most obvious relationship is that be- 
tween average current velocity and the 
most abundant particle diameter of a 
sediment—the mode. In addition to the 
mode, however, if the exact nature of 
the transport is restricted to either sus- 
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in Trinidad, British West Indies, and 
contiguous areas and is a preliminary 
description of the method of approach 
and results obtained by empirical testing 
of the relationship between size and sort- 
ing in water-borne sediments. 


ANALYTICAL PROCEDURE 


Two statistics form the basic data 
upon which the subsequent computa- 
tions are performed; they comprise the 


UNIQUE FREQUENCY CURVE, 
COMPUTED (2). 


Mo@ = 3.'5, graphical from U.F. curve. 
= 3.25, graphical from cumulative} 


PERCENT BY WEIGHT 


frequency curve. 


M® =3.8/6, by method 
of moments (17). 


MO 


2 4 


5 6 8 


® Scole 


Fic. 1.—Cumulative and unique frequency curve of typical Miocene sand, southern Trinidad 


pension, saltation, or traction, then the 
range of sizes is also closely determined 
by the current velocity. It appears from 
this that size and sorting are controlled 
by similar factors and hence that there 
should be a traceable, if not direct, rela- 
tionship between them. 

On this premise the author has for 
some years been compiling grain-size 
analyses both from published work and 
from his own investigations, and it is 
clear that the relationship can be traced 
fairly simply for water- and air-borne 
sediments; whereas, as would be ex- 
pected, gravitational and glacial deposits 
do not follow the same pattern. This dis- 
cussion is limited to the author’s results 


median diameter and percentile devia- 
tion, both measured in terms of the phi 
scale (Krumbein, 1936). The phi median 
(Md,) is obtained by reading the phi 
diameter at the 50 per cent point in the 
size distribution. It is obtained by direct 
reading from the cumulative curve and is 
a commonly used “‘average’’ diameter 
based upon geometrical position in the 
center of the distribution. It is not af- 
fected by the extremes of the range and 
is therefore particularly suited to “‘open- 
ended” distributions. The phi medians of 
the sediments used for this investigation 
approach the modal values more closely 
than does either the arithmetic or the 
geometric mean (fig. 1) and hence are a 
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closer approximation to this genetically 
important diameter than is any other 
average. It is, of course, much simpler 
to derive than the mode which can be 
satisfactorily obtained only from the 
unique frequency curve (Krumbein, 
1938; Brotherhood and Griffiths, 1947). 
The phi percentile deviation (PD,) 
is analogous to the more orthodox phi 
quartile deviation (Krumbein, 19360), 
but the phi quartile diameters are re- 
placed by the phi diameters at the 10 
and go percentiles. The formula for der- 

ivation becomes 
— Pwe 


PDs= 


(1) 
This statistic is used as a measure of 
sorting and represents half the spread be- 
tween the ninetieth and the tenth percen- 
tile diameter values. It has several ad- 
vantages over the phi quartile deviation 
while still retaining its ease of derivation. 
The percentile deviation encompasses 
some 8o per cent of the distribution, as 
compared with 50 per cent for the quar- 
tile deviation. It is equally applicable 
to open-ended curves and, as would be 
expected, is more sensitive to fluctua- 
tions in the fine (P,.,) and coarse (P,.,) 
ends of the distribution. It has, of course, 
a greater range in absolute value than 
the quartile deviation has, and this is also 
some advantage. It is obtained by read- 
ing the phi diameters at the tenth and 
ninetieth percentiles, respectively. 

It is important to emphasize here that 
the three diameter values—median, 10, 
and go percentiles—are mutually inde- 
pendent. Variation in one does not 
necessitate or imply variation in either of 
the others. Thus at the outset the two 
statistics—-median (average size) and 
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percentile deviation (sorting)—are both 
mathematically and petrographically in- 
dependent. Any relationship established 
between these two measures is therefore 
not implicit in the manipulations. The 
succeeding mathematical treatment is 
conventional and uses the standard sym- 
bols X and Y for the two variables. It is 
necessary, therefore, to decide which 
statistic shall be represented by the inde- 
pendent variable X and which by the 
dependent variable Y. It seems clear that 
the more accurately determined vari- 
able is Md, based on a single reading in 
the steepest part of most cumulative 
curves, and Md, is, therefore, designated 
X; by the same reasoning PD, is de- 
pendent on two readings in the as- 
ymptotic parts of the curves and hence 
is least accurate and becomes the de- 
pendent variable Y. It is interesting to 
note that this is in agreement with their 
petrographic significance, since the 
modal value, of which the median is a 
close approximation in this case, is much 
more precisely circumscribed than is the 
range by the current velocity. 

The phi medians and phi percentile 
deviations of over one thousand grain- 
size distributions of Tertiary sediments 
in Trinidad and contiguous areas were 
compiled, and these form the basic data 
upon which the subsequent manipula- 
tions were performed. In order to 
evaluate empirically the relationship be- 
tween size and sorting, scatter diagrams 
were constructed of average size (Md,) 
versus sorting (PD,), and a discernible 
trend was displayed in most cases (figs. 
2-12). The measure of correlation be- 
tween size and sorting can be usefully 
expressed by the correlation coefficient. 
The raw-data formula was used (Rider, 


1947, p. 48): 


(2) 
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where r is the correlation coefficient and 
> is the sign for summation, X and Y 
are Md, and PD,, respectively, and N is 
the number of samples. 

In this computation random selections 
of 15-100 samples were made from vari- 
ous sets of data, and, using r and the 
number of samples, the statistical sig- 
nificance of the correlation coefficient 
may be directly tested from Fisher's 
table (1948, p. 209). Most of the correla- 
tion coefficients computed during the 
present work proved to be statistically 
significant (P much less than 0.01 and in 
many cases less than o.0o1). Literally, 
this test is carried out by assuming that 
the “parent-population” of Md, and 
PD, values is uncorrelated and then 
computing how likely an r of the value 
obtained from the sample of Md, and 
PD, represented by the data is to arise 
from this uncorrelated population. In the 
present work this likelihood (P?) was 
found to be less than i in 100 (P < 0.01) 
or less than 1 in 1,000 (P < 0.001), and 
hence it is assumed that a correlation 
exists in the population from which the 
sample tested was drawn. Those samples 
which failed this test are discussed where 
relevant in the text. 

The presence of association once estab- 
lished, the next step is to compute the 
relationship between the variables. This 
was accomplished by the “method of 
least squares”’ (Mills, 1938; Fisher, 1948) 
and the line of “best fit” was computed 
by means of the normal equations 

=Na+bzX, 


where the notation is similar to that in 
equation (2) and a and 6 are the un- 
known coefficients to be determined. 
The spread of the values around this 
line of best fit was measured by the 
standard error of estimate (Sy), and 
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this was computed from the following 
formula (Mills, 1938): 


_ 3(¥?) (¥)—b2(X VY) 
(Sy)? = 


(4) 


The notation is the same as for equa- 
tions (2) and (3). 

The line of best fit and the estimated 
spread (+Sy) are plotted on the scatter 
diagrams and used as a standard of 
reference for investigating unusual de- 
partures from the linear relationship. On 
the basis of the implied assumptions that 
the relationship between size and sorting 
is linear, the trends may all be expressed 
in the following form: 


PDs= bMdgt+a. (5) 


DESCRIPTION OF THE RESULTS 


In order to shorten the description, it 
will be sufficient to outline one example 
in detail and trace the steps followed and 
the assumptions involved. Where subse- 
quent treatments are identical, only the 
computed trends need be dealt with dur- 
ing the description of the results. In some 
cases, however, deviations from the 
trends become important, and the distri- 
bution of the observations over the 
scatter diagram is then given in full, but 
the mathematical steps involved are only 
briefly stated. 


APPLICATION OF THE PROCEDURE 


Well A in Forest Reserve Field, south- 
western Trinidad, was cored through 
part of the Miocene sequence of sedi- 
ments. Sixty-five samples from these 
cores were subjected to grain-size anal- 
ysis, and the phi median and 10 and 
go percentile diameter values were read 
from the cumulative curves. The phi 
percentile deviations were obtained by 
substitution in formula (1). 

The values of Md, and PD, were then 
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plotted in a scatter diagram with Md, as 
abscissa (X) and PD, as ordinate (Y). 
The resulting trend is shown in figure 2. 
A sample of 25 pairs of values (Md, and 
PD,) was randomly selected from the 65 
for purposes of computation. Random- 
ness was insured by using standard 
tables of random numbers (Snedecor, 
1946; Fisher and Yates, 1948). The 25 
pairs of values were then tabulated, and 
the necessary sums (ZX, TY), sums of 
squares ({X?, Y¥*), and sum of products 
(ZXY) of the values were obtained on 
an automatic calculator. these 
values the correlation coefficient 7, the 
least-squares trend, and the standard 
error of estimate were computed by using 
equations (2), (3), and (4), respectively. 
In this case the sample of 25 paired 
values gave a correlation coefficient, r= 
0.9408. Assuming that the set of 25 
values is actually taken from an infinite 
number of such pairs of values (popula- 
tion) in which there is no correlation be- 
tween Md, (X) and PD, (Y), the chance 
of obtaining an 7 of this magnitude is 
much less than 1 in 1,000, i.e., P < 0.001 
(Fisher, 1948). In this sense, therefore, it 
appears that the measure of association 
between size (Md,) and sorting (PD,) is 
statistically significant. 

Computing the line of best fit by the 
method of least squares (eq. [3]), the 
following relationship results: PD,= 
0.6190Md, — 0.3510; and this is plotted 
in the scatter diagram (fig. 2). As a 
measure of the scatter of the data around 
this line, the standard error of estimate 
computed from equation (4) gives Sy = 
0.27; statistically, this implies that if the 
values belong to normal distributions, 
some 68 per cent of all the values should 
lie around the line of best fit and within 
1 Sy unit on either side of it. Parallel 
lines are inserted at 1 Sy unit on either 
side of the line of best fit (see fig. 2), and 
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it can be seen that, of these 25 pairs of 
values, more than two-thirds and nearer 
three-quarters fall within the specified 
limits. From the entire 65 values, how- 
ever, far more than one-third fall outside 
these Sy limits. 

A second set of 20 pairs was selected 
from the original data after the previous 
25 had been eliminated from considera- 
tion. The correlation coefficient derived 
from this second set of values gives r = 
0.9085, and the statistical significance is 
equally high (P < 0.001). The least- 
squares trend is PD, = 0.806Md, — 0.954, 
and the standard error of estimate 
is Sy = 0.49. It appears from this that 
the Sy value calculated from the first set 
of pairs is fortuitously low. 

From this one example it can be seen 
that the sorting varies with size, and 
from the computed trends it is also clear 
that, as Md, increases, PD, also in- 
creases. In millimeter units this signifies 
that sorting becomes poorer as the sedi- 
ments become finer over the size and 
sorting limits represented by these 
values, i.e., between the Wentworth 
grades of fine sand and coarse silt (i.e., 
Md, = 2.c-6.0). It appears convenient 
at this point to emphasize some of the 
basic assumptions involved in the an- 
alytical procedure which leads to this 
conclusion. It is not possible to establish 
these assumptions as valid at this stage, 
but at least a brief discussion will reveal 
the steps in the analysis necessary for 
their ultimate establishment. Three prin- 
cipal questions arise: Is it correct to as- 
sume that the trend is linear? What is 
the nature of the sampling error in- 
volved? What is the significance of the 
standard error of estimate as a measure 
of the spread of the values around the 
trend line? 

Examination of the data in the scatter 
diagram clearly indicates that, if only 
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those points with medians less than 4.5¢ 
and percentile deviations less than 3.0¢ 
are considered, the trend, if linear, would 
be much steeper than either of those 
computed. Such a trend is confirmed by 
a single value from the 65 beyond the 
limits of the scatter diagram (Md, = 
5.150; PD, = 4.60). The decrease in 
the slope of the line is, therefore, primari- 


straight lines. Without a much larger 
number of points in the region Md, = 
4-8¢@, it is not possible to decide between 
these alternatives, and, hence, for the 
present it is convenient to assume a 
linear trend as a first approach. 

The second question concerning the 
variation in trend line due to sampling 
error can be tested empirically and cer- 


PD 


*=First set of 25 samples 
PD = 0.6190 Md@ - 0.3510 


e=Second set of 20 samples 
= 0.803 Md - 0.940 


©=Samples not included in sets 


05 3 


5 6 7 


Md > 


Fic. 2.—Size versus sorting relationship, well A, Miocene, southwestern Trinidad 


ly due to those points lying between 
Md, = 5-8@ and PD, = 3-5¢. Several 
alternative explanations exist for such 
occurrences. The points beyond Md, = 
5¢ may be in error;? the trend is actually 
curvilinear; or the trend in the two 
ranges may be better expressed by two 


? A simple assumption would be to consider these 
points as representative of clayey sediments which 
had suffered flocculation. This possibility is dis- 
cussed more fully on p. 222. 


tain practical limits on the number of 
samples used for computation thereby 
established. Sixty-nine subsurface sam- 
ples from well B, Forest Reserve Field, 
southwestern Trinidad, were analyzed, 
and the resulting Md, and PD, values 
were plotted in a scatter diagram (fig. 3). 
Successive sampling in sets of 15, 20, 25, 
30, and 69 pairs of values was performed 
and the statistics computed in turn. The 
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results are summarized in table 1. It may 
be seen that the correlation coefficient r 
varies between 0.4054 and 0.5077, a 
small range compared with the change in 
the number of samples. In the case of the 
set of 15, however, a value of r (0.4524) 
would arise because of errors of random 
sampling from an uncorrelated popula- 
tion about 1 in 10. This can hardly be 
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considered significant in a_ statistical 
sense and would certainly not be evi- 
dence of a valid association in a geo- 
logical sense. With a set of 20 pairs of 
values, the position improves; but such a 
correlation as that listed in the table 
could arise about once in 20 times from a 
set of 20 pairs drawn from an uncorre- 
lated population. There is little real 
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trend 
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Fic. 3.—Size versus sorting relationship, well B, Miocene, southwestern Trinidad 


TABLE 1 
SIZE VERSUS SORTING RELATIONSHIP IN 69 SUBSURFACE CORES FROM 


WELL B, MIOCENE, TRINIDAD 


No. of 
Samples 


ivast-Squares Trend 


P De =0.6806Mdg—o. 2263 
PDo=2.9214Mds—8. 2680 
PDe=0.5859Mdg—0.0177 
PDoe=0.6584Mdg—o. 2745 


Standard 
Error of 
Estimate 


Correlation 
Coefficient 


©.4524 
©.5°77 
©. 4320 
0.4054 
°. 401060 


|| 
2 | 
a Line 
569 
25. °.90 <0.05 
: 0.40 | <0.05 
69 59 <0.001 
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change up to 30 samples; but, with 69 
pairs and a value of r = 0.4616, such a 
value would arise much less than 1 in 
1,000 from an uncorrelated population, 
and hence the association is statistically 
valid. 

Examination of the trend lines in 
figure 3 indicates the amount of variation 
to be expected solely from random sets of 
samples from the same set of data; lines 
1, 4, and 5 are in reasonably close agree- 
ment; line 2 is somewhat aberrant, but 
line 3 is so markedly different that it is 
clear that the sampling is “biased.” 
Analysis of the trend reveals that the set 
of samples (see fig. 3, samples ®) con- 
tained an unusually large number con- 
fined to Md, of 3-4¢ but with a very 
large variation in PD, (1.125-4.03¢ 
units). As will be shown later (p. 223), 
such a combination is a reflection of two 
events of geological importance combin- 
ing to obscure the average trend of 
change in sorting with grain size. Sta- 
tistically it may be considered that this 
peculiar set of values is from a different 
parent-population. 

The third question concerns the mean- 
ing of the standard error of estimate, Sy. 
Statistically, Sy is a measure of the de- 
parture of the values from the trend, or, 
more specifically, it measures the aver- 
age discrepancy between values of PD, 
computed from the least-squares trend 
equation and those actually observed. 
One of the more obvious reasons for 
such variation is due to “experimental er- 
ror,’’ wherein this error includes sampling 
error between individual samples and 
arises from choice of sedimentation units 
(Apfel, 1938; Otto, 1938). This has not 
been separately evaluated in the present 
case but would principally affect lami- 
nated sediments, in which, owing to the 
nature of the samples, many sedimenta- 
tion units would be included in one 


JOHN C. GRIFFITHS 


sample. This would imply that “silty”’ 
sediments with Md, between 4 and 6¢ 
would be the most likely types to show 
large variation. As such an expectation is 
not borne out in any of the sets of 
samples tested, it is apparent that the in- 
dividual sampling error is not very im- 
portant. It will, however, mean that from 
this cause, i.e., mixed sedimentation 
units, the value of Sy can never be re- 
duced below a certain unknown mini- 
mum. A second feature which adds to 
the comprehensive “experimental error”’ 
is variation due to analytical technique. 
This has been evaluated in the case of 
two different sand types in the Miocene 
sediments of southern Trinidad and gave 
as maximum experimental error 0.07¢ in 
Md, due to this cause, whereas in PD, 
the maximum so far evaluated is 0.24¢. 
Experimental error is here given in terms 
of the standard deviation determined by 
replicate analysis (Krumbein, 1934; Rit- 
tenhouse and Connaughton, 1944) and 
comprises errors due to analytical tech- 
nique, smoothing during construction of 
the cumulative curves, and reading 
errors. This can be broken down a little 
further by noting that the error in the 
case of the median is not nearly so large 
as that in the case of the percentile devia- 
tion. Here again the errors can be split in 
the latter case into those associated with 
Prog and those associated with Pog. 
From experience with the Trinidad Mio- 
cene sediments it can be stated that the 
errors associated with P,,, are not very 
different from those associated with the 
median, whereas P,.4, as can be seen 
from figure 1, is very greatly affected by 
very small differences in smoothing and 
reading in the flat part of the curve. This 
alone probably accounts for the greater 
part of the experimental error in PD, 
and, in turn, therefore, with the variation 
in Sy. It seems clear from this discussion 
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that Sy will always be large in sediments 
with open ends and particularly where 
the sorting over the whole distribution is 
poor when compared with the central 
sorting. 

Statistically, Sy can be precisely cir- 
cumscribed only in the case of ‘“‘normal”’ 
distributions, and in such cases it can be 
shown that 68 per cent of the samples 
should fall within +1 Sy unit from the 


THE MIOCENE SEDIMENTS OF 
SOUTHERN TRINIDAD 


WELL SAMPLES FROM SOUTHWESTERN 
TRINIDAD 


Further sets of surface and subsurface 
samples from the same and different 
areas of Miocene sediments in southern 
Trinidad maintain a similar relationship 
between average size (Md,) and sorting 
(PD,). The generalized relationship al- 
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C-G ondL. Subsurface samples, Forest Reserve 
H. Subsurface samples, Central Trinidad 


K. Surface samples, S.E. Trinidad 


H 


Field, S.W. Trinidad 
S.E. Trinidad 
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Fic. 4.—Size versus sorting relationship, Miocene, southern Trinidad 


trend lines, 95.4 per cent of the samples 
within +2 Sy units from the trend, and 
only 3 in 1,000 outside the limits of + 3 Sy 
units from the trend. In the case of well L 
(fig. 5) these prerequisites are not fulfilled, 
principally because the distributions of 
both Md, and PD,, but chiefly the 
former, are not “normal.”’ This is one of 
the reasons why the strict mathematical 
treatment, using, say, 2 Sy units, cannot 
be adopted in the present analysis. 


ways indicates that sorting become 
poorer with decrease in grain size. Some 
results based on several sets of samples 
are illustrated in figure 4, and the sta- 
tistical data are summarized in table 2. 
The trends, however, indicate that the 
rate of change varies widely. ‘Two effects 
can be partially isolated; first, an over- 
all change in average size of a set of 
samples shifts the trend line across the 
graph but maintains approximately the 
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same slope. This would account, for ex- 
ample, for the changes in trends C, E, F, 
and K. Second, a change in slope is 
generally due to “bias”’ in the sampling, 
a larger number of sands giving a steeper 
trend than a set of samples with many 
clays. This accounts for the divergence of 
trends G and J. Trends D and H are, as 
can be seen from the table, not statisti- 
cally significant and could represent sets 
of samples from uncorrelated popula- 
tions; it appears most likely, however, 
that these aberrant trends result from 
experimental error and, in the case of D, 


that, on the basis of 531 samples, such a 
value of r would arise much less than 1 in 
1,000 by chance from an uncorrelated 
population. It therefore appears valid to 
deduce that in the Miocene sediments of 
southern Trinidad the percentile sorting 
becomes poorer with decrease in average 
grain size. The relationship computed as 
a trend line gives the following equation: 
PD, = 0.5574Md, + 0.1754; Sy = 0.77. 
This trend is inserted in figure 4, trend 
line L, and the general parallelism with 
trends C, E, F, and K is clearly estab- 
lished. 


TABLE 2 


SUMMARY OF TRENDS IN SIZE VERSUS SORTING IN MIOCENE 
SEDIMENTS OF SOUTHERN TRINIDAD 


Least-Squares Trend 


PDe=0. 5567Mde—0. 7799 
PDe=0.1034Mde+o. 9905 
PDe=0. 5019 Mdg —0. 6296 
PDo=0.4283Mdg—0.6776 
PDo= 1.0931 Mdg— 2.0829 
PDe=0. 2812Mdg+o. 9688 
PDe=1.1440Mdg — 2.2752 
PDe=0.5797Mdg—o. 5311 

* Not statistically significant. 

t Surface samples. 


from sampling bias (see pp. 222 and 
227). 


WELL L FROM SOUTHWESTERN TRINIDAD 


The complexity of the variation can be 
illustrated by a study of a set of 531 
samples from well L, Forest Reserve 
Field in southwestern Trinidad. This 
well was continuously cored from the sur- 
face (approximately 500 feet) to 5,500 
feet, a stratigraphic range of some 5,000 
feet. The scatter diagram of size (Md,) 
versus percentile sorting (PD,) is shown 
in figure 5. The correlation coefficient for 
this set gives a value of r = 0.7585; test- 
ing for statistical significance indicates 


In figure 5 the standard error of esti- 
mate (Sy) is inserted as two lines parallel 
to the main trend given above and 1 Sy 
unit away in each case. As a first ap- 
proximation, it is assumed that values 
outside +1 Sy unit are sufficiently dif- 
ferent to demand explanation, whereas 
within 1 Sy unit of the trend the varia- 
tion is due essentially to individual 
sampling and experimental errors (opera- 
tional errors, reading errors, etc.). The 
samples which fall outside the unit Sy 
band fall into four areas, two below the 
band and two above. It will, however, be 
more convenient for our present purpose 
to consider these four groups in order of 


i Well P No. of 
| 
f e y Samples 
Cc <o.o1 | 29 
D °.50 ©.0405 >o.10* | 36 
E 0.70 | 0.7423 <o.oo1 50 
F ©.9207 <0.001 5° 
G 0.47 ©. 6009 <0.001 84 
H 0.48 ©. 3404 >o.05* | 31 
J | 0.48 | 0.6566 <0.001 43 
Kt 0.51 0.6264 <0.001 37 
i) 
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average size, and they, therefore, will be 
grouped as follows: (1) samples below the 
band of ‘‘average trend”’ and with Md, > 
5%; (2) samples above the band of 
“average trend’’ and with Md, > 4¢; (3) 
samples below the band of “average 
trend”’ and with Md, = 3-4¢; and (4) 
samples above the band of “average 
trend” and with Md, = 3.5-4.5¢. 


DEVIATIONS FROM THE TREND IN WELL L 


1. Samples below the band of “average 
trend” and with Md, > 5@.—If the linear 
trend PD, = 0.5574Md, + 0.1754 rep- 
resents the “true’’ relationship between 
size and sorting in the samples from 
well L, then those samples which fall be- 
low the trend are better sorted than 
would be expected from the equation. 
The samples below the trend and with 
Md,> 5@ are therefore fine-grained 
sediments of silt and clay sizes, which are 
better sorted than would be expected. 
Examination of the cumulative curves of 
these samples clearly indicates that this 
better sorting is due to flocculation dur- 
ing analysis, and the better sorting of 
these samples is therefore only apparent. 
It appears most likely that these samples, 
if the size analysis had been more closely 
controlled, would all show poorer sorting 
and perhaps even finer average grain size 
(larger value of Md,) if flocculation had 
not occurred during the analysis. 

2. Samples above the band of “average 
trend” and with Md, > 4¢.—-Again on 
the assumption that the trend represents 
a true relationship, these samples repre- 
sent sediments which contain a large pro- 
portion of silt and clay sizes and with a 
much poorer sorting than would be ex- 
pected. This variation is due to an as- 
sumption made by the author during the 
computation of PD, for these samples. 
In all these samples the cumulative 
curves never intersected the go percen- 


tile, and in the majority they passed the 
10@ ordinate of the grain-size distribu- 
tion curve (approximately 1 yu, the fine- 
grained limit of analysis in the experi- 
ments) running parallel to the abscissa. 
This would imply that they were con- 
tinuing indefinitely below the go percen- 
tile. Pragmatically, the author chose 14¢ 
as the go percentile intersection with the 
cumulative curve. In all these samples, 
therefore, Py, = 14¢. The variation in 
Prog in most of the Miocene samples 
from southern Trinidad is small, and, 
therefore, PD, remains relatively con- 
stant over large changes in Md,. The 
distribution of samples in figure 5 shows 
the marked artificiality of this pro- 
cedure by the linear parallelism between 
PD, = 4.5-5.0¢, whereas Md, varied 
from 4 to 7¢. The peculiarity of this 
distribution is particularly emphasized 
by the rarity of any samples in the range 
PD, = 3-75-4.25¢ and Md, = 5.5-79. 
The following argument can be used to 
develop a correction which appears to be 
petrographically logical and mathemati- 
cally satisfies the trend. It has many 
times been stated that mineral composi- 
tion varies with grain size (see, for ex- 
ample, Pettijohn, 1949), the clay min- 
erals being the dominant phase in the 
finer-grained sizes. In the Miocene sedi- 
ments of southern Trinidad the clay 
minerals comprise illite and kaolinite; 
and hence, with a known mineral com- 
position, the grain size of the finest frac- 
tions can be approximately predicted. 
Marshall (1931) has shown that, over a 
wide range of ionic saturation, clays tend 
to show a median value around 12- 
12.5, and these Miocene sediments are 
unlikely to prove exceptions; hence the 
cumulative curves which run parallel to 
the abscissa through 10¢ probably turn 
toward the 100 per cent ordinate around 
12-12.5¢. If this assertion is correct, then 
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it means that the P,,, for the samples in 
this group are in error some 1.5~2¢ units, 
and hence PD, is in error 0.75-1.00¢ 
units. Applying this as a “blanket’’ cor- 
rection to the samples herein discussed 
immediately brings them within the belt 
of average trend. Hence it appears that 
the majority of these samples are unlike- 
ly to be exceptions to the rule expressed 
in the formula given above, after due 
allowance is made for experimental error 
and the author’s incorrect initial assump- 
tion as to the value of Pog. 

3. Samples below the band of “average 
rend” and with Md, = 3-4¢.—-These 
samples are better sorted than would be 
expected from the formula, but, unlike 
case 1 above, the sediments are very fine 
sands and not clays and silts. Clearly, 
therefore, the same explanation is not 
sufficient to account for their departure 
from the trend. 

Detailed examination of this group of 
sediments shows that, of the total 42 
samples, 20, or some 50 per cent, are 
from a single 250-foot sand body and the 
remainder are stratigraphically widely 
distributed. The 250-foot sand body is 
known to be exceptionally uniform on 
the basis of petrographical character- 
istics and to extend over a considerable 
area in which these characteristics are 
relatively constant. The electrical log 
pattern of this sand is singularly con- 
sistent over the area, as would be ex- 
pected from the petrographic uniform- 
ity, and the sand body is used as a 
marker horizon on this account. It has 
been called a “sheet” sand because it is 
the nearest approach to this type of sand 
body in the Miocene sediments of 
southern Trinidad. It very often occurs 
immediately below an unconformity, 
and, where absent, it is suspected that 
it was cut out during the period of 
erosion represented by the unconform- 
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ity, which, in turn, indicates the biggest 
single break during Miocene sedimenta- 
tion in southern Trinidad. 

On the basis of this evidence, one may 
deduce that this sand represents the top 
sand of a Miocene sedimentary cycle 
and, iurthermore, that the duration of 
the period of sedimentation indicated by 
this exceptionally uniform sediment is a 
period of stillstand, one of the few well- 
marked such periods in the Miocene 
sedimentational history of the area. It 
was, therefore, a period of reworking and 
relative “‘tectonic’’ quiescence which led 
to an extremely well-sorted sediment, as 
indicated by its grain-size distribution. 
The period of stillstand was followed by 
a relatively violent renewal of sedimen- 
tation and its accompanying “‘tecto- 
nism,’’ and hence this is a well-defined 
unit in the sequential history of Miocene 
sedimentation in Trinidad. 

By extending this reconstruction, one 
can suggest that wherever such sands 
are found in the Miocene sediments, a 
period of stillstand and reworking oc- 
curred. By the argument given above, no 
other period is so well marked in the 
historical sequence, and hence the other 
samples of this group are more widely 
distributed stratigraphically (and geo- 
graphically, as can be shown from other 
areas sampled). Hence the divergence of 
this group of samples from the average 
trend is attributed to relatively unusual 
conditions of sedimentation during the 
Miocene sedimentation in southern 
Trinidad. 

4. Samples above the band of “average 
trend” and with Md, = 3.5-4.5¢.—-These 
samples are more poorly sorted than 
would be expected from the average 
trend, and they represent very fine sands 
and coarse silts in size. In general, they 
are also more mixed, containing a rela- 
tively larger amount of coarser and finer 
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sizes than do the average Miocene sedi- 
ments. This, of course, is the factual ex- 
planation of their position on the scatter 
diagram. The geological significance of 
this departure from the trend is also of 
some considerable significance in terms 
of sedimentational history. 

The average Miocene sediment of 
southern Trinidad is a graywacke-type 
sediment (Krynine, 1948) and is relative- 
ly consistent in grain-size distribution. 
The sediments have been characterized 
as deltaic in environmental type (Grif- 
fiths, in press) and possibly represent 
the subsea deposits of the Miocene 
equivalent of the present-day Orinoco. 
Their size and mineral composition are 
typical, and, owing to detailed sampling, 
the petrographic characteristics of these 
sediments have been thoroughly evalu- 
ated. Sediments which belong in the pres- 
ent group contain elements not typical of 
the average deltaic type. They contain 
an admixture of minerals and sizes which 
differ from the true delta type. The ac- 
cessory (heavy) minerals suggest that 
these elements are derived from a local, 
possibly Oligocene and older, source 
within the area of southern Trinidad. By 
a somewhat less direct argument it can 
be shown that the coarser sizes are also 
locally derived. 

The sedimentational characters of 
delta deposits of large rivers are char- 
acteristic of transported materials in 
which the transport has had around 
1,000 miles in which to impress its efiects 
on the materials. The very thoroughly 
investigated Mississippi delta sediments 
are an outstanding example (Krumbein 
and Aberdeen, 1937; Russell and Russell, 
1939). Such sediments are usually repre- 
sentative of the suspension load of a 
large river and, as such, have a closely 
circumscribed upper size limit. Such 
rivers characteristically transport ma- 
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terials which vary from fine sand to clay, 
and the modal size is generally silt 
(Krumbein and Aberdeen, 1937; Russell 
and Russell, 1939; Tercier, 1939). The 
Miocene sediments of southern Trinidad 
fulfil these requirements almost exactly, 
although the amount oi clay is generally 
in excess of that of the typical Mississippi 
deposit. This has been attributed by the 
writer to the fact that the Trinidad de- 
posits represent the subsea facies of the 
Miocene river delta, a suggestion sup- 
ported by the well-marked and oft- 
repeated marine incursions throughout 
the Miocene succession. Hence on the 
basis of this characteristic sedimenta- 
tional history it is very difficult to under- 
stand how coarser materials can be 
added to the deposits by way of the main 
river; the most obvious alternative, con- 
firmed by other features, is the local 
source of such material. The major diffi- 
culty lies in explaining how an area with- 
out marked relief can be the locus of 
deposition and a source area at one and 
the same time. In the case of the Trini- 
dad deposits there is a simple solution 
to this enigma—-the source of the ab- 
normal material lies in the outbursts of 
local mud volcanoes which abound in 
the area of Miocene deposition in south- 
ern Trinidad. Their effects can be 
studied during their present sporadic 
activity, which is both subaerial and 
submarine (Kugler, 1938). The ma- 
terials ejected by these volcanoes com- 
prise a size range from boulders to clay, 
and they generally indicate Miocere, 
Oligocene, Eocene, and Cretaceous sedi- 
ments as provenance. Hence these out- 
bursts supply the correct type of ma- 
terial at the right place and time (con- 
tinuous throughout the Miocene to 
Recent at least). The volume of material 
supplied is not easy to compute, but, in 
the author’s opinion, at least some 10 
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per cent of Miocene sediments could 
have come from such a source. 

The mud volcano ejectmenta supply 
two dominant sizes of material to the 
Miocene sediments: a large amount of 
clay and some much smaller amount of 
sand and gravel; the latter has been 
stressed in this discussion only because it 
is easier to trace petrographically and is 
more obvious to both geologists (pres- 
ence of pebbles and boulders of older 
rocks) and paleontologists (source of the 
so-called “reworked” faunas of Oligo- 
cene, Eocene, and Cretaceous “age” in 


from the scatter diagrams, it appears 
probable that a linear trend would re- 
main over the median size range of 2.5- 
6@ (and perhaps to 8@) in the Miocene 
sediments of southern Trinidad, and for 
this reason such an approach is con- 
sidered justifiable by the writer. Varia- 
tions from this statistically and petro- 
graphically valid trend are of some im- 
portance, and, when experimental varia- 
tion is removed, any deviation that re- 
mains is capable of explanation on a 
geological basis similar to that outlined 
in the above cases of groups 3 and 4. It is 


TABLE 3 
STRATIGRAPHIC SAMPLING PATTERN IN MIOCENE SEDIMENTS OF WELL L 


No. of 


Member Samples 


Upper Morne L’Enfer........ 16 
Lower Morne L’Enfer. . . 172 
Forest sands........ 153 
{ Upper... 58 
Cruse; Middle... .. 75 
\Lower..... $7 


Total Miocene............] 531 


clays at the base of Miocene sedimentary 
cycles). It is probable that modern tech- 
niques of analysis for clay-size minerals 
could confirm the presence of Oligocene 
clay detritus in reasonably large amounts 
at certain horizons in the Trinidad Mio- 
cene sediments. 

This deviation from the trend is also, 
therefore, of very great geological im- 
portance. These sediments are less well 
sorted because of volcanic incidents 
which occurred during their sedimenta- 
tional interval, and they therefore indi- 
cate interference with the normal pro- 
cedure of deposition after fluvial trans- 
port of some 1,000 miles. 

When all those deviations are removed 


Remarks 
Incompletely represented 
Well represented 
Well represented 
50 per cent represented 
Well represented 
Well represented 
Not represented 


Approximately two-thirds 
represented 


now possible to examine the variation in 
trends in space and time over the area of 
southern Trinidad throughout the Mio- 
cene. 


STRATIGRAPHICAL VARIATION IN TREND 


The scatter diagram illustrated in 
figure 5 includes samples from five dis- 
tinct members of the Miocene succession, 
and the exact stratigraphic range of 
sampling is summarized in table 3. It 
can be seen that, in terms of stratigraphic 
members, five are adequately repre- 
sented, and the sampling interval runs 
approximately one sample to ro feet, ir- 
respective of member, i.e., the sampling 
is statistically random and has no strati- 


ia ve 
| 
Feet Feet per 
88 5-5 
1,628 9.5 
1,357 8.9 
435 7-5 
650 8.6 
815 14.2 2 
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graphic bias. The linear trends for the 
five members adequately sampled from 
well L are displayed in figure 6, and the 
statistics are tabulated in table 4. From 
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ably consistent and that the variations 
are sufficiently small to be neglected. The 
Upper Cruse samples, however, show a 
trend (no. 3 in fig. 6) which, although 


parallel to the other trends, is moved 
toward the finer sizes. No sand samples 


the figure and the table it can be seen 
that trends 7A, 1B, 2, and 4 are reason- 


Stratigraphical variation, Well L 


Line 
1A. Morne L Enfer, lower: Ist set 50 
iB. " i2nd" 
2. Forest Series 
3. Cruse, upper 
3A. 49): 
4. middie 
5. Cruse, lower 
6. All members of Well L: 


| 


4 5 7 
Md 


Fic. 6.—Size versus sorting relationship, Miocene, southwestern Trinidad. Stratigraphical variation 
in well L. 


TABLE 4 
STRATIGRAPHIC VARIATION IN WELL L 


‘ | No. of 


Member Least-Squares Trend 
Sa 


Morne L’Enfer Lower. . 
Forest series 


PDe=o. 7832Mdy— °. ». 

| PDe=0.7790Mdg—0.8360 
PDo=0.8617Mdg—0.9823 
PDo=0.7681Mdg—1.6918 
P Do = 127.06Mdg — 479.7864 
PDg=0.6651Mdg—0.0751 
PDo=0. 209Mde+t. 


1754 0.77 


| Lower.... 


All horizons... . 


| 
| 
Cruse) widdle 


*NC, not computed. 
t Twenty-three selected samples of Upper Cruse sands with Mdy <4@ 


= 
6 
3A 
Ze 
2 46 
| 
| ‘ Dies 
q 
3 3A 
i 
P 
0.75 | 0.8722 | 50 | 
| 0.80 0.7831 | 50 | <0.001 
0.73 | 0.7176 | 50 | <o.0o1 
| | | <o.01 
| NC | NC 23 
| 0.84 7668 50 <0.001 
0.98 ©.3112 | 40 <0.05 
| 0.7585 | 531 | <0.001 
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such as those present in the remainder of 
the succession occur at this horizon, and 
so the size effect is marked in this trend. 
The only sands which are represented in 
these samples are from the peculiarly 
uniform body of sand already discussed. 
Further to illustrate the effect of this 
sand on the trend, all samples from this 
interval with Md, greater than 4¢ were 
eliminated, and the remainder (23 
samples) were used to calculate a trend 
line (no. 3A, fig. 6). It is clear that 
mathematically there is no association 
between size (Md,) and sorting (PD,) in 
this sand. Petrographically the sand 
shows a tremendous range in PD, for a 
very small change in Md,, and this is, of 
course, emphasized in the peculiar trend 
line. The only component varying in this 
case is P,.,, and small changes in the 
amount of clay-size material would, of 
course, have a very great effect on the 
Pyeg, Which could, in part, be the ex- 
planation of the peculiar variation. It 
was also found, however, that in this 
sand the value of experimental error far 
exceeded that of variation from sample 
to sample in terms of Md,, and hence, 
by the method of size analysis used, no 
variation in Md,, oF Pog could be 
established as a valid change in size 
from sample to sample. This, of course, 
serves to emphasize the uniformity in 
this unusual sand body. 

Another exceptional trend line is 
shown by samples from the Lower Cruse 
member, and here the slope of the trend 
line (no. 5, fig. 6) has changed consider- 
ably. By analogy and from our earlier 
discussion, this implies the presence of 
clay because the slope is less steep than 
the average trend (no. 6, fig. 6). It is no 
surprise, therefore, that in well L most 
of the Lower Cruse samples comprise 
silt and clay. The larger part of the re- 
duction in slope in this case is also il- 
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lusory because the clays had in many 
cases flocculated, falling well below the 
sorting measure predicted for their size. 
Hence this trend line is almost entirely 
due to experimental error in terms of 
lack of control in the experimental 
technique. 

The trend line for the complete se- 
quence is given in line 6, and the com- 
bined effects of trend lines 3 and 5 have 
resulted in only a slight change in slope 
in trend 6 as compared with lines 7A, 
1B, 2, and 4. As would be expected, there 
is very little displacement due to size. 
On this basis it is probable that trend 
lines 1A and 7B are nearer the “true” 
relationship between size and sorting 
than are any other trends computed. 
This is borne out by the fact that less 
flocculation occurred in the clayey sam- 
ples from this horizon.’ 

It is possible to carry this analysis of 
stratigraphic variation in trend a step 
further, as a number of other sets of 
samples have been taken from similar 
horizons in other wells from the same 
field. In this way it is also possible to 
illustrate the variation within members 
over a very restricted geographical area. 

The trend lines of size versus sorting 
in the Forest series from two wells and 
the Middle Cruse from three wells are 
illustrated in figure 7, and the relevant 
statistics are summarized in table 5. Two 
other trends are inserted for comparison, 
the first representing the average trend 
for all horizons in well L and the second 
trend based on the average of two sets of 
samples from the Lower Morne L’Enfer 
horizon from well L, the “true” trend 
line. It can be seen from the scatter dia- 
gram that, in terms of size (i.e. position 


3 Owing, no doubt, to the base saturation of the 
clays and their associated fluids, as these samples 
for the larger part lie in the ground-water horizon 
rather than in the highly saline oil-bearing horizons. 
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on the graph), there is no significant less steep. It would, of course, be tempt- 
variation between trends, but they can _ ing to infer that the trends for the Forest 
be grouped into two sets on the basis of series are different from those in the 
slope. The steeper set—lines 1,2,and7—- Middle Cruse of these wells; but, by in- 
are sandier sediments and approximate _ serting lines 6 and 7, it can be seen that 
the “‘true”’ trend or relationship, whereas _ the difference is probably much less sig- 
the others—lines 3, 4, 5, and 6—are all nificant, being probably attributable to 


Line 
Forest series 


3.Middle 


6. All horizons::- ‘L 
*7 Morne LEnfer, 


* Average of 2 sets of 50 samples. 


6 
Md 


Fic. 7.--Size versus sorting relationship, Miocene, southwestern Trinidad, variation within horizons 


TABLE 5 
STRATIGRAPHIC AND LIMITED GEOGRAPHIC VARIATION WITHIN MEMBERS 


| No. of | 


Horizon Least-Squares Trend | Sa maples, 


Forest series PDe=o. 806 Mde— 0.954 

Forest series | 
Middle Cruse PDg=0.5019Mdg—0.6296 
Middle Cruse |§ PDg=0.4283Mdg—0.6776 
Middle Cruse | PDg=0.6651Mdg—0.0751 
All horizons | PDg=0.5574Mde+o.1754 
Morne L’Enfer 8820 

Lower* 


* se two sets of 50 


i 
6 5 
3 
q 
a E 
+ 
| 
j 
2 3 4 
= — = = = = — — — — — — 
Well | P 
A.. | 0.49 | 0.9085 | 20 <0.001 
L. | 0.73 0.7176 | 50 <0.001 
| 0.70 | 0.7423 | 50 <0.001 
| ©.51 | 0.9207 | <0.001 
L. | 0.84 07668 § 50 | <o.001 
| 0.77 | 0.7585 531 | <0.001 
| 0.78 | 0.8276 | 100 | <o.oor1 
| 
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more clays and more flocculation of the 
clays from the Middle Cruse than those 
from the Forest series. 

In summary, then, any variation 
from horizon to horizon stratigraphical- 
ly or within horizons and over geo- 
graphic limits of 1-2 miles is not brought 
out in these examples, and, as far as the 
present analysis is concerned, all strati- 
graphic variations are to be attributed to 
other than purely geological factors. 
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two are based on surface samples and 
the third is from a well. The geographic 
range from east to west is about 50 miles. 

The trend lines computed from these 
sets of data are exhibited in figure 8, and 
the statistical parameters are summa- 
rized in table 6. Trends 2, 3, and 5 in 
figure 8 are parallel, but, from the posi- 
tions of the trends, 3 and 5 are from 
samples which are finer grained than 2. 
In terms of geographic position the sur- 


|. Lower Morne L'Enfer Well L 
2.All members Welll 
3.Kalmaposs Beach Cliff 


4. Well J, S.E, Trinidad 
5.Galeota Beach Cliff 
6. Well H, S. Central Trinidad 


Dien 


5 6 7 


Fic. 8.—Size versus sorting relationship, Miocene, southern Trinidad, geographical variation 


GEOGRAPHIC VARIATION IN TREND 


The Miocene sediments of southern 
Trinidad extend from the extreme west 
to the eastern land limits of the island. 
The examples so far discussed are nearly 
all from a single field (Forest Reserve, in 
southwestern Trinidad). One example 
has also been computed on the basis of 
subsurface samples from south-central 
Trinidad and three from the extreme 
southeastern part of the island. Of these, 


face samples from southeastern Trini- 
dad are finer grained than those from 
well L in’ southwestern’ Trinidad, 
but the relationship between size and 
sorting is otherwise almost identical. 
Trend 1, which represents the selected 
“true” trend from studies of horizons in 
well L, is steeper than all three of these 
trends, and, as it nearly coincides in 
position with trend 2, it also contains 
coarser-grained sediments than 3 or 5. 
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This change in grain size from east to 
west is confirmed by comparing sets of 
samples from the two areas in terms of 
their cumulative curves. 

Trend 4 in figure 8 represents the rela- 
tionship between size and sorting (PD,) 
based on 43 subsurface samples from 
well J in southeastern Trinidad. It is 
somewhat steeper than the other trends 
in the figure, but its general position co- 
incides with trends 7 and 2. On the basis 
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cases the correlation coefficient is low 
and not statistically significant. 

On the basis of these examples it can 
be seen that, while the majority of 
samples from the southeast of Trinidad 
are finer grained than those in the south- 
west of Trinidad, the relationship be- 
tween median size and percentile sorting 
remains consistent. Any variations from 
the general trends are due to either 
sampling or experimental error or both. 


TABLE 6 


GEOGRAPHIC VARIATION IN TREND IN MIOCENE 
SEDIMENTS, SOUTHERN TRINIDAD 


Least-Squares Trend Sy 


PDo=0.5574Mdo+0.1754 | 9.77 | 0. 7585 
©. 8276 


PDo=0.7811Mdg—o. 8820 


PDe=0.5797Mde—o. 5311 
PDe=1.144 Mdg—2.2752 | 


| 
PDg=0. 2812Mdg+o. 9688 


°.78 


0.48 


0.40 | 
0.48 | 
| 


* Average of two sets of 50 samples. 
t Surface samples from southeastern Trinidad. 


0.6264 | 


| 
| No. of | 


Remarks 
Samples 


r 


All horizons well L 

Lower Morne L’Enfer 
well L* 

Kalmapass Beach Cliff 

Well J, southeastern 
Trinidad 

Galeota Beach Cliff 

Well H, south-central 
Trinidad 


} Correlation coefficient not statistically significant; P nearly equal to, but slightly greater than, 0.05. 


of earlier comments, such steep trends 
are usually indicative of sampling bias; 
in general, it can be stated that, if few 
clays occur in the set, the trend will be 
steep, and, if many clays occur, the 
trend will be much less steeply inclined. 
On the basis of individual grain-size dis- 
tribution, it can be stated that the 
samples from well J are all very fine sands 
or coarse silts. 

Trend 6 is much less steeply inclined 
and indicates quite clearly the presence 
of a large number of clay samples in the 
set, and, by analogy with the discussion 
on well L, these clays have flocculated 
during analysis. Examination of the indi- 
vidual cumulative curves confirms this 
conclusion. It is noteworthy that in such 


In summary, then, the average trend ex- 
pressing the relationship between size 
and sorting in the Miocene sediments of 
southern Trinidad remains consistent 
stratigraphically over a range of some 
7,000 feet and geographically over a dis- 
tance of some 50 miles. 


THE OLIGOCENE DEPOSITS OF 
SOUTHERN TRINIDAD 


The Oligocene deposits of southern 
Trinidad have been subdivided into two 
facies, the Nariva “formation,” a near- 
shore sediment, and the Cipero ‘‘forma- 
tion,’ representing a shallow-water 
marine environment (Stainforth, 1948). 
The subsequent analysis refers solely to 
the Herrera sands and associated clays 


|| 
| | | 
Well | 
100 | <0.001 | 
| 

Kt <0.001 | 
H.. 0.3404 >o.ost 
| | 

| 
\ 
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from the Cipero facies of the Oligocene. 
The Herrera sands are quartzose gray- 
wackes with much less matrix than the 
Miocene sands. In addition, they are 
coarser grained and contain a relatively 
large proportion of gravel (conglomer- 
ates). The most striking contrast be- 
tween Miocene and Oligocene deposits is, 


Miocene series. The least-squares trends 
of a number of Miocene and Oligocene 
samples are illustrated in figure 9, and 
the relevant statistics are summarized 
in table 7. Superficial examination of the 
Oligocene trends shows their similarity 
to the average trend from well L in the 
Miocene and their even closer identity 


5 


MIOCENE 
| Lower Morne L'Enfer Well L 


2 All members Well L 
3. Galeota Coast section 


OLIGOCENE 


4. Subsurface samples 


5. Subsurface samples 


Md 


Fic. 9.—-Size versus sorting relationship, Miocene and Oligocene, southern Trinidad 


however, the dominance of silt in the 
Miocene and the paucity of silt in the 
Oligocene. These two groups of sedi- 
ments are clearly representative of two 
very different environments of deposi- 
tion. 

For the purpose of comparison, the 
first set of Oligocene samples has been 
chosen from sequences which do not con- 
tain gravel grain sizes in any quantity. 
Despite the selection, the Oligocene sedi- 
ments are still, on the average, coarser 
than their lithologic equivalents in the 


with the sandier sediments of the ‘‘true’”’ 
trend. The two Oligocene trend lines 
(lines 4 and 5, fig. 8) are displaced toward 
coarser average sizes, and their steepness 
reflects the preponderance of sands in 
the sets of samples chosen. Here again, 
then, over the range of size from ap- 
proximately 2.5@ (0.177 mm.) to 6@ 
(0.0156 mm.) median diameters the per- 
centile sorting becomes poorer with de- 
creasing grain size. 

Taking into account the sets of 
samples with gravel fractions shows an 


| 
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entirely new type of distribution (fig. 
10). In this scatter diagram it is clear 
that no single linear trend can adequate- 
ly represent the relationship. Continu- 
ing the assumption of linearity as a first 
approximation, two trends are neces- 
sary, and they have been computed and 
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inserted in figure 10; the relevant statis- 
tics for the four sets of Oligocene data 
are given in table 8. 

From these data it can be seen that 
over a range in median diameter of 2.5— 
6¢@ the percentile sorting becomes poorer 
(PD, increases in absolute value) with 


TABLE 7 


COMPARISON OF MIOCENE AND OLIGOCENE RELATIONSHIPS OVER 
THE SAME SIZE RANGE 


Trend* Least-Squares Trend Sy 


0.77 
°.78 
©.49 
49 


PDe=0.5574Mdo+o.1754 
PDe=0.7811Mde—o. 8820 

PDe=0.5194Mdg—0.6056 | 
PDg=0.7083Mdo—c. 3420 | 
9063 | 


No. of 


‘Samples Horizon 


Miocene SW. Trinidad 
Miocene SW. Trinidad 
Miocene SE. Trinidad 
Oliogcene 
Oligocene 


* See fig. o. 
+t Surface samples, Galeota coast section = M of table 6. 


Line 


Oligocene “Sands” 
7. PD= 2.4791 -0.5371MdQ@;Olig. "Graveis" 


Md 


Fic. 10.—Size versus sorting relationship, Oligocene, southern Trinidad 
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decrease in grain size (increase in value 
of phi median); from median diameter 
2.5-1.0, however, this relationship is re- 
versed, and it can be seen that, with in- 
creasing median diameter (decreasing 
phi value), the percentile sorting be- 
comes poorer. It seems apparent, 
therefore, that sediments of median di- 
ameter around 2.5@ (0.177 mm.) are 
the best sorted, whereas the sorting be- 
comes poorer with both increasing and 
decreasing size from this diameter. This 
conclusion refers solely to the Oligocene 
sediments analyzed and is confined to a 
grain-size range of 1-6@ (0.500-0.0156 
mm.). 
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ment which is coarse grained and con- 
tains gravel grain sizes in considerable 
proportion. Superficial examination of a 
few thin sections of this material sug- 
gests that it is also a graywacke. This is 
another marine sediment, but one which 
suggests a very different environment 
from either Miocene or Oligocene of 
Trinidad. Through the courtesy of the 
late Dr. A. Senn, the author was enabled 
to collect a representative series of sam- 
ples from the Lowes Ridge type section 
(Senn, 1940, p. 1557) of the Scotland 
series in Barbados during a vacation 
visit in 1942. A plot of these samples in a 
scatter diagram is illustrated in figure 11, 


TABLE 8 


Least-Squares Trend 


PDe=0.7083Mdg—o. 3420 
PDe=0.8500Mdg—0. 9063 
PDe=1.2107Mdg— 2.3221 
PDo= 2.4791 —0.5372Mdg 


No. of 


Samples | P 


Remarks 

<0.001 
<0.001 
<o.10t 


Oligocene sands 
Oligocene sands 
Oligocene sands 
Oligocene gravels 


* See fig. 10. 
t Not statistically significant. 


The reversal in association between 
phi median diameter and phi percentile 
diameter is shown by the change in sign 
of the correlation coefficient. It must be 
emphasized, however, that the level of 
statistical significance (P < 0.10) at- 
tached to the negative correlation co- 
efficient (— 0.3788, trend line 7) indicates 
that such a relationship would arise by 
chance from random sampling of an 
unrelated population about 1o times in 
100. Statistically, therefore, this set of 
samples is insufficient to establish the re- 
versal of trend. 


THE EOCENE SEDIMENTS OF BARBADOS 


The Scotland series of Eocene age in 
Barbados (Senn, 1940) is another sedi- 


and here again the distribution cannot 
logically be fitted to a single linear trend. 
The two trends, together with the two 
Oligocene trends, are inserted in the 
figure. Here again the sediments around 
2.5-2.75@ (0.177-0.149 mm.) appear to 
be the best sorted, and sorting becomes 
poorer when it proceeds toward coarser 
and finer sizes. The relevant statistics 
are summarized in table 9, and it can be 
seen that the trend is reversed and that, 
while still not statistically significant at 
the level chosen for this report (P < 
0.02), such a value of r would arise by 
chance from an uncorrelated population 
less than twice in 100 examples. The 
similarity of the Trinidad Oligocene 
trend to that of the Barbados Eocene 


TRENDS IN THE OLIGOCENE OF SOUTHERN TRINIDAD 

| 
0.55 | 0.5913 | 30 
°.49 0.7092 30 (| 
0.44 0.8115 28 | 
| 0.31 —0.3788 | 20 | ; 

| 
j 
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serves further to confirm the reality of 
the reversal of trend. Increase in the 
number of coarser-grained sediments 
would almost certainly bring this level 
of significance into line with the others 
tabulated in this report.‘ 


THE SUPERFICIAL GRAVELS OF 
THE GUIANAS 


The Superficial Gravels of both Dutch 
and British Guiana have been penetrated 
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depositional environment most nearly 
resembles that of the Barbados Eocene 
sediments. The Guiana gravels are con- 
sidered to have been “dumped” after a 
short period of transport; the grain size 
varies from gravels to clays, and the 
size-sorting relationships on the basis of 
two sets of samples, one from each well, 
are displayed in figure 12. The relevant 
statistics are given in table 10, to- 
gether with representative data from the 


3q _ 40, 


3a. Gravels 


4 Barbados Eocene “Sands* 
40. 


e Borbados Samples 


© Intersections * best 
sorted sediments 


0 


3 5 


Fic. 11.—Size versus sorting relationship, Oligocene, Trinidad, and Eocene, Barbados 


in two deep wells, and cores from both 
wells have been analyzed for grain size. 
The sediments from the British Guiana 
well have been described in detail 
(Griffiths, 1942); on the basis of mineral 
composition they are arkoses, but their 


4This problem is being actively pursued at 
present, and initial results offer confirmation of these 
conclusions. 


Oligocene and Miocene of Trinidad and 
the Eocene of Barbados. The size versus 
sorting relationship in the British Guiana 
gravels is somewhat similar to that in the 
Miocene in slope but is displaced mark- 
edly toward coarser average grain sizes, 
as would be expected. The general rela- 
tionship suggests that percentile sorting 
becomes poorer with decrease in grain 
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size over the whole range tested (o-5¢, 
1.0-0.0312 mm.). The average grain size 
of the Guiana gravels is very much 
coarser than that of the Miocene and Oli- 
gocene sediments of Trinidad but is 
reasonably similar to the average size of 
the Barbados samples. No reversal of 
trend around a median diameter of 
2.5@ (0.177 mm.) is noted in the Guiana 
gravels. This suggests that the deposi- 
tion of the Guiana gravels was different 
from that of the Barbados Eocene de- 
posits, despite their similarity in average 
grain size. 
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The most outstanding difference be- 
tween the Guiana samples and those of 
Trinidad and Barbados is in the value of 
the standard error of estimate (Sy). The 
Guiana gravels have much larger Sy 
values than any of the other trends, and 
this appears to be true of both Dutch 
and British Guiana samples. In a geo- 
logical sense it appears reasonable to 
infer that these departures are also sig- 
nificant. Considering only the trends 
over the sand fractions from the Miocene 
of southern Trinidad reduces the value of 
Sy (Sy generally less than 0.50), and 


TABLE 9 


COMPARISON OF BARBADOS EOCENE AND TRINIDAD OLIGOCENE RELATIONSHIPS 


Least-Squares Trend 


..| PDg=0.7083Mdg—0. 3420 
.... PDe=o.8500Mdg—o. 9063 
3..... PDg=1.2107Mdg—2. 3221 
3A 
4..... PDg=0.6193Mdg—0. 5596 
4A. PDo=2.6474—0. 5608Mdg 
* See fig. 11. 
t NC, not computed. 


The association of average size (Md,) 
and percentile sorting (PD,) in terms of 
the correlation coefficient is of the same 
order in the British Guiana samples as 
that in the Miocene and Oligocene of 
Trinidad and the Eocene of Barbados 
over the size range of 2.5—5¢@. It is sta- 
tistically significant at the 1 per cent 
level (P < 0.01). In the Dutch Guiana 
samples the level of significance is lower 
(P < 0.05), but this is probably due to 
the fact that fewer samples are available 
(16 as compared with 28, respectively). 
The association is therefore statistically 
significant and, in addition, by analogy 
with the Trinidad and Barbados samples, 
geologically significant. 


Ne. of 
Samples 


Remarks 

30 Oligocene “‘sands” 
30 | Oligocene “sands” 
28 Oligocene “sands” 
20 | Oligocene “gravels” 
34 «| Barbados “‘sands” 
10 Barbados “gravels” 


the paucity of clays in the sets of Oligo- 
cene samples also exemplifies this fea- 
ture (Sy = +0.5). The same is true of 
the Barbados gravelly deposits (Sy = 
0.50). In other words, sandy deposits 
show a size-sorting relationship, with de- 
partures from the trend measured by Sy 
of less than 0.5. The median diameters of 
the Guiana gravels vary between o and 
4¢, and yet even in these sandy sedi- 
ments the Sy is approximately twice as 
large (Sy =+ 1.0@), indicating that the 
departures from the trend are twice as 
large. The significance of this contrast 
cannot be categorically stated, and its 
interpretation will therefore be treated in 
the discussion of the results (p. 240). 


| 

Trend 

©.5013 
©.49 ©. 7092 
©.44 | 
0.31 —0.3788 | { 
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DISCUSSION OF THE RESULTS 


INTERPRETATION OF THE TREND 


The relationships between size and 
sorting in the sediments herein described 
lead to the following general conclusions: 

1. Over the range in average size 
(Jid,) of o-6¢ (1.0-0.0156 mm.) the per- 
centile sorting (PD,) varies with size. 

2. Between size limits of o and 2.5¢ 
(1.0-0.177 mm.) sorting becomes poorer 
(PD, increases in value) as the average 
size becomes larger (Md, decreases in 
value). 
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cern sediments which are water-borne 
and water-deposited. From preliminary 
testing, similar relationships exist in 
wind-borne and wind-deposited sedi- 
ments, but glacial and gravitational de- 
posits show very different relationships 
(usually vague to none) between size and 
sorting. 

It would appear, therefore, that a 
generalized relationship between size and 
sorting exists in water-borne and water- 
deposited material, and even on the 
basis of this preliminary testing the con- 
clusions appear to be amply justified. It 


TABLE 10 
COMPARISON OF THE GUIANA GRAVELS WITH OTHER CARIBBEAN SEDIMENTS 


Least-Squares Trend 


7523Mdgt+o. 46081 
5300Mdg+ 2.0788 
5574.Mdg+o.1754 
7811 Mdg—o. 8820 
7083Mdg—o. 3420 
2107 Mdg — 2. 3221 
4791 —0. 5372Mde 
6474—0. 5008 Mdg 
6193Mdg—o. 5596 


PDe=0. 
PDg=0 
PDe=0 
PDe =o 
PDe=0 
PDe=1 
PDg=2 
PDe=2 


I. 
2 
4 
6 
8.. 


t NC, not computed. 


3. Between 2.5 and 6¢ (0.177-0.0156 
mm.) sorting becomes poorer (PD, in- 
creases in value) as the average size be- 
comes smaller (Md, increases in value). 

4. best-sorted sediments lie 
around 2.5@ (0.177 mm.) in average 
diameter and possess percentile sorting 
values generally less than 1¢@ (i.e, 
PD, < 

Such conclusions are confirmed by 
Trask’s observations on size and sorting 
based on empirical studies of sediments 
(1932) and on Inman’s analytical ap- 
proach based on the hydrodynamics of 
water transport (1949). It must be re- 
membered that all these remarks con- 


—0.7020 


° 


| | 
| No. of | 
| Samples 


Remarks 


British Guiana 
Dutch Guiana 
Miocene, Trinidad 

Miocene, Trinidad 
Oligocene, Trinidad 
| Oligocene, Trinidad 


| 


Oligocene, Trinidad 
Eocene, Barbados 


5465 Eocene, Barbados 


also seems reasonable to assume a linear 
trend over restricted size intervals, but 
only so long as the limitations are fully 
recognized and taken into account in any 
conclusions which may be drawn from 
the data so tested. By using millimeter 
and semilogarithmic plots, Inman (1949, 
fig.1) obtained a curvilinear trend for 
various kinds of deposits; and, if the 
trends displayed in the present work are 
replotted on a semilogarithmic scale in 
terms of millimeter units, similar curvi- 
linear relationships are, of course, ob- 
tained (fig. 13). The use of the phi scale, 
therefore, tends to simplify the analysis 
and does not distort the conclusions. 


|| 
Line* | | 
= - - — 
| 1.14 | | 
| 0.05 | 0. 5238 
| 0.77 | | 5 
| | 78 0.8276 I 
.| | 0.55 | 9.5013 
| 0.44 o.8115 
NCt | 10 <0.02 
* See fig. 12. 
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TRANSFORMATION OF THE LINEAR TREND 


FROM PHI TO ITS EQUIVALENT IN 
MILLIMETERS*® 
Krumbein derived the phi scale from 
the relationship £ = 2-*, where & = 
diameter in millimeter units (1938). 
Therefore, 


this may be similarly transformed to 
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it therefore follows that 


a 22 


2°M 


Therefore, in millimeter units, PDg = 
1/2°Md?. Commencing with PD, = 
bMd,+a results in the millimeter 
equivalent, PDg = 2°/ Md}. 

The generalized equations in terms of 
millimeter units are exhibited as actual 
trends in figure 13, using the empirically 


_ (—logio Po) / (logw 2) + (ogi Pio) / (logy 2) 
PDs= 


which may be reduced to 


— Pw) log PDg 
log io 2 


2 log 1 2 


where PDg is the geometric percentile 
deviation in millimeter units and PDg = 
V Pyo/ Poo (see analogous transformation 
of QD, to So = QDg; Krumbein and 
Pettijohn, 1938, p. 234). Therefore, 


log 1 


But the geometric percentile deviation 
PDg = V Pyo/ Poo (millimeter units), and 
it follows that log V Pyo/ Po = 0.301 
PD, or Pye = 10°? % = 
Hence PDg = 2°"%; substituting this 
relationship in the linear trend PD, = 
bMd, — a, we get 


Px z 


— loge 10 -log Md;; 


$I am indebted to Professor E. Johnson for this 
mathematical operation. 


derived values for two sets of samples 
from the Miocene sediments of southern 
Trinidad, two from the Oligocene, and 
two from the Eocene of Barbados. It 
will be readily appreciated that the use of 
phi units and the assumption of linear 
trends over restricted size ranges serve to 
simplify materially the resulting an- 
alytical treatment and its subsequent in- 
terpretation. 


SOME IMPLICATIONS OF THE TREND IN 
SIZE VERSUS SORTING 


Commencing, now, with the linear 
trend PD, = bMd, — a, we can trans- 
form this trend to PD, = bMd, by 
changing the origin to some value of PD, 
equal to —a. This equation can then be 
further reduced to 


= ot Prop—Pwe= 20M dg; 


but, as Md, is a close approximation to 
the modal diameter in the sediments 
which gave rise to this relationship (e.g., 
the Miocene sediments of southern Trini- 
dad), then Pyog — Prog = kMo,, where k 


‘ 
‘ 
and, as 
— 
but 
Md, — 
logy 2 
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is some constant approximately equal to 
2b. In simple terms, Poy — Prog « 
f(Mo,), or the spread between the phi 
percentiles varies as a function of the 
phi modal diameter. For purposes of 
exploration let us assume that the cur- 
rent velocity remains constant and that, 
therefore, Mo, is reasonably constant. 
Then any variation in P,., must be ac- 
companied by a similar variation in 
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velocity and bearing a sand-size material 
as its preferred load (Mo, around very 
fine sand size in, say, the Mississippi 
River) there is a sudden influx of clay- 
sized material into the stream, then 
some of the coarser sand sizes must be 
deposited. This is, of course, a much 
simplified and idealized case but offers, 
in part at least, a possible explanation 
for local deposition of fine and medium 


|. 531 Somples, Well L 
2.100 Samples, Morne UEnfer Well 


3. Oligocene, 


Fine Sands and SGilts 
Coarse Sands 
Fine Sands and Silts 
Coarse Sands 


5. Borbados, 
6. 


sediments) 
3 and 4+ Oligocene, Trinidad 
5 and 6= Eocene, Barbados 


e@ intersection of trends (best sorted 
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Fic. 13.—Size versus sorting relationship, Md, versus PD, in millimeter units 


For example; suppose is in- 
creased from 2 to 4¢; then P,., must also 
be increased 2¢ units to maintain the 
relationship. In terms of sedimentation, 
increasing P,., means decreasing the 
diameter of P,. in millimeters, i.e., 
adding clay size to the material. Simi- 
larly, increasing P,., means decreasing 
the diameter value of P,, in millimeters. 
Such an implication suggests that if in a 
stream traveling with uniform average 


sand in streams by collapse of clay-sized 
stream banks. 

The converse of this operation is not 
so easily visualized, i.e., suppose simi- 
lar initial conditions and an influx of 
coarse sand. Then by similar reasoning 
the P,, in millimeters increases and the 
Prog decreases in value. For uniform 
conditions to hold, Py. must also de- 
crease in value, and some clay therefore 
must drop out. It seems more reason- 
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able to suppose that the stream would 
fail to pick up the coarser sand and, 
therefore, that such additions would have 
little ultimate effect on the size distribu- 
tion of the transported material. In ac- 
tual practice it seems even more prob- 
able that sudden local additions of ma- 
terial would affect the average velocity 
and, therefore, the value of Mo,; but, as 
the size of the materials ultimately de- 
posited by such large streams as the 
Mississippi is relatively constant (and 
such restricted distributions can be 
demonstrated for the Mississippi and the 
Miocene sediments of southern Trinidad, 
for example), then local conditions play 
a very small part in the ultimate size 
distribution of delta deposits of large 
rivers. 

It is also evident that one of the most 
important loci of deposition is in large 
deltas, and this can be independently 
demonstrated by examination of the 
present-day conditions of sedimentation 
and those of the geologic column. Hence 
this depositional type is of very consider- 
able significance, and, if approximate 
volumetric controls are estimated, the 
rivers must be of the order of 1,000 miles 
long for a delta of some 200 miles across. 
The Amazon, Congo, Ganges, Hwang 
Ho, Mississippi, Nile, Orinoco, and 
Yangtze Kiang are all of this order of 
magnitude, whereas the Indus, Irrawad- 
dy, and Rhine are only a little smaller— 
both in length and size of delta. In all 
these cases it appears that the trans- 
ported material has had a very long time 
interval during which relatively undis- 
turbed conditions have operated on the 
materials to produce sediments of simi- 
lar size distributions. The size-sorting 
relationships in only two have been ex- 
amined, but it is obvious that under such 
physical conditions as those represented 
by these deposits such a relationship is 
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likely to reach its maximum develop- 
ment. Indeed, it appears probable that 
delta deposits represent as near equi- 
librium conditions as we can ever hope to 
expect within the average sedimenta- 
tional environment. The size-sorting rela- 
tionship in these cases is therefore as well 
expressed as we can expect. It may be ob- 
jected that a more efficient agent, acting 
for a shorter time interval, would bring 
about a similar equilibrium, but the 
actual size range would then be very 
different, as can be seen by examining 
beach sands or dune sands (Inman, 1949; 
and author’s observation), and, further- 
more, such deposits are far more con- 
spicuous at present than they are in the 
geologic past as represented in the sedi- 
mentary column. 

On this basis it is considered that the 
size-sorting relationship in any sediment 
represents an estimate of how near to 
equilibrium conditions the sediment has 
approached. Thus water-borne deposits 
will generally be better sorted than 
glacial deposits, owing to different hydro- 
dynamic properties of the transporting 
medium, but within the water-borne de- 
posits themselves, when the size range is 
established, the size-sorting relationship 
reflects the “amount”’ of sorting the de- 
posit has undergone, and in this way it 
may become a very useful tool in predict- 
ing the environment of deposition. This 
aspect of size-sorting relationships as a 
whole is in process of evaluation at 
present. 

In summary, then, it appears likely 
that the existence of size-sorting relation- 
ships in sediments can be used to deter- 
mine the nature of the transporting 
medium, and, once this is decided, the 
perfection of the relationship can be used 
to determine the kind of depositional 
environment to some extent. On this 
basis it is possible to explain the dis- 
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crepancy in size-sorting relationship be- 
tween the Trinidad Miocene sediments 
and the Superficial Gravels of the 
Guianas. Both are deposited in water, 
but the Trinidad deposits indicate a 
much greater perfection in size-sorting 
relationship on the basis of deviations 
from the trend (value of the standard 
error of estimate, Sy) than do the Guiana 
gravels. The gravels have not, therefore, 
undergone transport, either in intensity 
or in time, equivalent to the Miocene 
sediments of Trinidad. The heterogene- 
ity of sizes in the Guiana gravels im- 
plies lack of sorting intensity, and hence 
the transport was short, and reworking 
after deposition was relatively nil. In 
other words, the material was ‘““dumped”’ 
into the basin of deposition and covered 
rapidly by succeeding debris. Such was 
the depositional history deduced from 
the petrography of these deposits (Grif- 
fiths, 1942), and apparently the size- 
sorting relationship amply confirms those 
conclusions. 

Interpretation could be continued 
further, even on the basis of the present 
evidence, but it is prudent to await a 
thorough evaluation of the size-sorting 
relationship in terms of transporting 
fluid and depositional environment by 
empirically analyzing a much larger 
number of cases before extrapolating be- 
yond the limits of the present investiga- 
tion. 


4 


SUMMARY AND CONCLUSIONS 


A study of the average grain size in 
terms of the phi median diameter (Md,) 
and uniformity of size distribution in 
terms of the phi percentile deviation 
(PD,) has shown that there exists a 
definite relationship between size and 
sorting in many different sediments 
which are water-borne and water-de- 
posited. By confining attention to re- 
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stricted grain-size ranges, their relation- 
ship may be expressed in terms of linear 
equations of the general form PD, = 
Md, — a. 

Using this relationship in sets of sedi- 
ments which are known to be water- 
borne, the deviation from the trend be- 
yond that expected due to “‘chance’’ 
causes (+1 unit standard error, or +1 Sy 
unit) can be used to characterize sedi- 
ments which have not followed the aver- 
age relationship. In the sand sizes these 
deviations have been attributed to inci- 
dents in the sedimentational ‘‘tectonics” 
which accompany the deposition of the 
material and which are presumed to have 
caused the marked departures from the 
trends. In the finer sizes such deviations 
are generally a reflection of flocculation 
in the clay fractions. 

The trends are shown to vary with 
both sandiness (number of sand and 
clay samples in the set) and with coarse- 
ness (average size of any set of samples). 
In sets which comprise predominantly 
sand samples, the trends tend to be 
steeper than in those which contain 
many clays, i.e., the slope of the trend 
line is changed. Differences in average 
sizes of entire sets of samples lead to 
translation of the trend without nec- 
essarily affecting its slope. 

Trends in the Miocene sediments of 
southern Trinidad show no great varia- 
tion either stratigraphically over some 
7,000 feet or geographically over some 
50 miles, apart from variations due to 
sampling bias, as explained above. Oligo- 
cene trends are very similar to the Mio- 
cene within the same size ranges, but 
over the entire range of sizes displayed 
in the Trinidad Oligocene the size-sort- 
ing relationship is curvilinear and indi- 
cates that sorting becomes progressively 
poorer with increasing and decreasing 
average size from the best-sorted sedi- 
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ments, which possess diameters around 
2.5 (0.177 mm.). Barbados sediments 
of Eocene age show relationships smiilar 
to those of the Oligocene of Trinidad. 

The Guiana gravels (Pliocene, Grif- 
fiths, 1942) show an unusual trend as 
compared with either the Trinidad or 
the Barbados samples. Here the trend is 
less well defined and is not in accord 
with the size-sorting relationships of 
either Trinidad or Barbados sediments, 
although the average size of the sedi- 
ments is similar to that of the Barbados 
Eocene samples. The size-sorting trend 
line is parallel to those of the Trinidad 
Miocene trends, despite the fact that the 
size range is more in accord with the 
Barbados samples. 

The Sy value in the Guiana samples is 
twice as large as that in the Trinidad 
or Barbados samples. This is considered 
to be due to lack of sorting in the 
Guiana gravels as compared with those 
from Trinidad and Barbados. Such a 
conclusion is in accord with its deduced 
depositional environment based on other 
petrographic evidence. 
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A very brief and hypothetical exam- 
ination of the implications of the size- 
sorting relationship suggests that it may 
prove to be a useful quantitative tool in 
deciding the nature of the transporting 
medium and the depositional environ- 
ment if future investigations now in hand 
serve to confirm the present findings. 
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CAUSES OF THE GLACIAL LAKE STAGES 
IN SAGINAW BASIN, MICHIGAN* 


J HARLEN BRETZ 
University of Chicago 


ABSTRACT 


Because the volume of Glacial Grand River varied by several ten thousands of second-feet between the 
maximum discharge during a retreat of the Saginaw, Huron, and Erie lobes and the minimum discharge dur- 
ing a readvance, the erosional and transportational ability of the river varied. An advance of the ice front, 
accompanied by slackened erosion in the Grand River outlet for ponded water in front of these lobes, was 
marked by essentially static lake levels, which are recorded by shore lines. A retreat was marked by erosion- 
al deepening of the outlet and falling lake levels. Seven successively lower shore lines were made in the Sagi- 
naw Basin while six moraines were being laid down. The seventh moraine, to correspond with the lowest 
shore line, is presumed to lie submerged in Lake Huron. 


INTRODUCTION 


Decades of study of the records of gla- 
_cial lakes in the Huron and Erie basins, 
chiefly by those devoted students, Frank 
Leverett and Frank B. Taylor, have un- 
raveled a complicated succession of 
changing levels and shifting outlets, a 
history discussed briefly in every current 


American text on historical geology. Two 
contrasted causes produced these chang- 
ing lake levels in front of the shrinking 
glacial lobes. One was change in location 
of the ice front: retreats opening new 
lower outlets and thus lowering lake 
levels, or readvances closing outlets 
and producing higher levels determined 
by the next-lowest available discharge 
routes. The other cause, yielding only a 
series of consecutive lowerings, was deep- 
ening of outlets in use. The problem at- 
tacked in this paper is the occurrence of 
pauses in such deepening, the results of 
which were stationary lake levels and 
definite shore lines. These pauses alter- 
nated with times of active outlet erosion 
and sufficiently rapid lowering of lake 
levels to leave no recognizable shore fea- 
tures. We are therefore seeking the rea- 
son for the second cause. 


* Manuscript received July 10, 1950. 


The Saginaw Basin is a branch or limb 
of the Huron Basin, projecting south- 
westward into Michigan from the larger 
and deeper feature. During the recession 
of the Wisconsin ice sheet, this basin de- 
termined a large subordinate lobe, whose — 
oscillatory shrinkage is recorded by 
moraines distributed almost the whole 
distance from the Indiana-Michigan 
state line northeast to Saginaw Bay. 
About 50 miles from the head of the 
bay, retreat began to expose a surface 
which sloped toward the shrinking ice, 
and free drainage of meltwater westward 
was replaced by standing water at the 
ice front. The outlet used by this lake 
water was the Grand River channel, 
which headed on the divide between 
present Lakes Hiron and Michigan and 
led westward for approximately 100 
miles to the contemporary Glacial Lake 
Chicago in the Lake Michigan Basin. 

Seven different shore lines, succes- 
sively lower toward the east, record 
seven pauses in the deepening of this 
channel, while ponded glacial water 
stood in front of the retreating Saginaw 
lobe.2 Why did not the westward-dis- 


2 Only once during this history was there a retreat 
sufficient to open temporarily another lower outlet, 
and none of these seven shore lines belongs to that 
episode. 
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charging water continuously deepen the 
Grand River outlet channel? 

Monograph 53 of the United States 
Geological Survey (1915) contains most 
of Leverett and Taylor’s matured inter- 
pretations of Michigan’s Pleistocene 
geology. Although each later published 
(1939) revisions, corrections, and addi- 
tions to the interpretations in the mono- 
graph, none of these alters the succession 
or the relations set forth in the mono- 
graph for the Saginaw Basin. The only 
explanation for the alternate lowerings 
and still-stands to be found in all their 
work was Taylor’s brief suggestion 
(Leverett and Taylor, 1915, p. 360) that 
a stoping process, twice initiated at the 
western end of Grand River channel by 
two drops in Lake Chicago’s level, pro- 
gressed eastward up the channel, even- 
tually affecting static water levels in 
the Saginaw Basin. 


SALIENT POINTS IN THE LAkE HIS- 
TORIES, FROM LEVERETT AND 
TAYLOR 


Leverett and Taylor grouped the re- 
cessional moraines of the Southern Pen- 
insula of Michigan into four systems, all 
of them represented in the Saginaw lobe 
deposits. The last two, the Lake Border 
and Port Huron systems, the only two 
concerned in this study, are remarkably 
deployed in the Saginaw Basin into sepa- 
rate moraines describing concentric arcs 
rudely parallel with the outline of the 
present bay. Thirteen of these moraines 
are recognized for the Lake Border sys- 
tem and three for the Port Huron. The 
first nine of the earlier system were land- 
laid throughout. The remainder of the 
sixteen, lying east of the divide, were 
built partly or wholly in the ponded 
water whose discharge was westward 
across the first nine to be deposited. 

When the glacial front stood at the 
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Flint (8th) moraine of the Lake Border 
system, Lake Maumee at the western 
end of the Erie Basin had abandoned its 
outlet to the Wabash River at Fort 
Wayne, Indiana, and was using a newly 
uncovered, lower spillway northward 
and westward across the “thumb” of 
Michigan. This, the Imlay River, skirted 
the southern edge of the Saginaw lobe 
and entered the Grand River channel 
near Maple Rapids, just west of the con- 
vex front of the lobe. There was as yet 
no lake in the Saginaw Basin. 

With the retreat from the Flint mo- 
raine came the first exposure of eastward- 
sloping terrain and therefore the first 
possibility of standing water in the basin. 
The lake was described as a small, cres- 
cent-shaped water body, moraine- 
bounded along its convex westward side, 
ice-dammed along its northeast concave 
margin, and discharging at Maple Rapids 
into the Grand River channel. According 
to the interpretation in the monograph, 
Imlay River presumably abandoned that 
lower part of its channel now occupied 
by Hayworth Creek (fig. 1) when ice 
withdrew from the Flint moraine, and 
crossed this moraine 10 or 12 miles east 
of Maple Rapids, entering the lake and 
adding its volume to the outlet river. _ 

When the ice front readvanced to 
build the Owosso (oth) moraine, this 
lake was obliterated, Imlay River was 
turned back into its abandoned lower 
course, and, at the apex of the lobe, the 
new moraine for a few miles of length was 
piled up on the back of the Flint moraine. 

Retreat from the Owosso moraine re- 
established Early Lake Saginaw, and a 
new short-cut for Imlay River was soon 
discovered across both the Flint and 
Owosso moraines northwest of the city 
of Flint. The lake level was held sta- 
tionary long enough for a shore record 
to be made, then fell between 10 and 12 
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feet as Glacial Grand River trenched in- 
to the Owosso moraine at Maple Rapids, 
pausing after that fall in level to consti- 
tute Second Lake Saginaw. 

Four faint water-laid moraines (Hen- 
derson, West Haven, Chesaning, and 
“moraine north of Chesaning’’)? on the 
lake bottom inside the Owosso appear to 
record that many more oscillatory shift- 
ings of the ice front before the Lake Bor- 
der system was complete. The algebraic 
sum of these oscillations was a general 
retreat, and during this time lower land 
became exposed in the “‘thumb,” the re- 
entrant angle between the Saginaw and 
Huron lobes. This finally brought about 
an abrupt drop in Lake Maumee to the 
level of Second Lake Saginaw, the two 
lakes becoming one in the same sense 
that Lakes Michigan and Huron today 
are one. Abandonment of all the Imlay 
River channel marked this fusion, which 
formed Lake Arkona. No change in the 
Saginaw Basin itself occurred immedi- 
ately except a broadening of the lake. 

But Lake Arkona had three different 
levels during its history, the total fall of 
its surface from Second Lake Saginaw’s 
level (with which Arkona’s first level 
coincided) having been about 15 feet. 
Thus there were two pauses in channel 
deepening at Maple Rapids during 
Arkona time, alternating with two times 
of channel erosion. 

The next forward movement of the 
ice front, recorded by the Port Huron 
(14th) moraine, overrode the strait in the 
re-entrant angle and forced the Erie por- 
tion of Lake Arkona to rise 30 feet or 
more to a spillway (Ubly) hard against 


3 Taylor was a bit uncertain in his identification 
of two of these. They have very faint relief on the 
lake plain but are accompanied by boulder belts 
which he interpreted as residual from wave erosion. 
The spacing intervals between the definite moraines 
of the region also seemed to require this interpreta- 
tion of these features. 
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the ice front at the tip of the re-entrant. 
As during Maumee-Imlay times, the dis- 
charge entered the lake in the Saginaw 
Basin and went on down Grand River. 

This dismembered Erie portion of 
former Arkona, called Lake Whittlesey, 
endured long enough for well-developed 
beaches to form on the landward side. 
They are presumably contemporaneous 
with the lowest Arkona beach in the 
Saginaw Basin. But nowhere have these 
beaches been found on the frontal slopes 
of the Port Huron moraine, the marker 
for the lowest Arkona-Whittlesey ice 
dam’s greatest advance. 

The retreat from the Port Huron mo- 
raine was greater than any earlier with- 
drawal of the lobe fronts had been. It 
cleared the region south of Syracuse, 
New York, sufficiently to offer the 
ponded water in Ohio, Michigan, and 
Ontario an escapeway to the Mohawk 
lower than the head of the Grand River 
channel, although that had been trenched 
into perhaps 30 feet since the days of 
Early Lake Saginaw. 

Lake Wayne, this eastward-draining 
successor to Whittlesey and lowest Ar- 
kona, endured only during the maximum 
of the retreatal phase following the Port 
Huron advance. 

The next advance, to the Bay City 
(15th) moraine position, promptly closed 
the eastern dischargeway and re-estab- 
lished westward drainage along Grand 
River channel. The new level in the 
Saginaw Basin, Lake Warren, stood 
about 5 feet below the lowest Arkona 
beach, recording that much channel 
deepening in the Arkona-Warren inter- 
val. Lake Warren left two definite shore 
lines, about 15 feet apart vertically. In 
places a weaker intermediate Warren 
shore has been found. During the life of 
the lake the two remaining members of 
the Port Huron morainic system—the 
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Bay City and the Tawas (16th)—were 
built and record two more oscillations of 
the front of the retreating Saginaw lobe. 

Retreat from the Tawas moraine coin- 
cided with lowering to the last stage of 
Lake Warren. The writer will argue, in 
this paper, that to maintain that level, 
there theoretically must have been an- 
other advance but that the moraine 
marking its culmination was entirely 
water-laid and is still submerged in the 
water of Lake Huron and Georgian Bay. 

Retreat from this theoretical moraine 
reopened the eastern dischargeway, and, 
from that time on, no closure of this 
route to the Mohawk by ice advance ever 
occurred. After Lake Warren's lowest 
stage, the Grand River channel ceased to 
be a glacial lake outlet. 


INADEQUACIES AND INFELICITIES IN 
THE PRECEDING INTERPRETATIONS 


That a stoping process could occur 


along the entire length of the Grand 
River channel, which is wholly in drift, 
seems difficult to accept. That there was 
adequate time for developing “grade”’ 
in the channel floor is a question. The 
most pointed objection, however, is that 
where Lake Chicago experienced two 
lowerings, of about 20 feet each, the 
Saginaw Basin lake series records six 
lowerings, totaling about 50 feet. Tay- 
lor’s suggestion of a causal relation seems 
impossible. 

Although six moraines were laid down 
after the first lake’s record was made in 
the Saginaw Basin and there were seven 
successively lower lake levels with west- 
ward discharge, Leverett and Taylor 
never attempted correlation of either 
stationary or falling levels with retreats 
or with readvances of the ice front. 

The Wayne and Arkona stages were 
pointed out as coinciding with two stages 
of glacial retreat, Whittlesey and War- 
ren, with two readvance climaxes. But 
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Arkona’s history is defined, by implica- 
tion at least, as covering the time of 
three retreats and readvances, and War- 
ren’s as covering two or three. The writer 
holds that the only lake stage properly 
connected with the maximum of a re- 
treat was that of Wayne. Had not an 
eastern outlet been uncovered in this 
retreat, it would logically have produced 
only another, still lower, Arkona stage 
in the following readvance. 

Although the Port Huron moraine 
was recognized as marking the maximum 
ice advance during the life of Whittlesey 
and lowest Arkona lakes, the entire lack 
of their shore lines on that moraine is but 
vaguely explained. ‘‘ Either the ice itself 
formed the shore, or readvance of the ice 
or outwash or the erosion of streams issu- 
ing from the ice obliterated the marks 
that the waves produced on the face of 
the new moraine”’ (Leverett and Taylor, 


1915, Pp. 379)- 
THESIS OF THIS STUDY 


1. Oscillations of the front of the 
Saginaw, Huron, and Erie lobes caused 
the oscillations between times of active 
deepening and of pauses in such deepen- 
ing of the Grand River outlet channel. 

2a. Lake levels remained sufficiently 
unchanged during glacial advances to 
permit the making of shore-line records. 

2b. Beginning at the maximum of an 
advance and the deposition of a moraine 
and continuing during a retreat, the out- 
let river deepened its channel, and lake 
levels steadily fell until retreat ceased 
and another advance began. 

3. The cause of this alternating be- 
havior of the outlet river was alternating 
maximum and minimum volume of dis- 
charge from the melting ice. The river 
volume during a retreat was increased by 
the total additional melting that pro- 
duced the recession. When the ice front 
readvanced into territory previously 
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abandoned, as much water was denied 
the outlet river, per mile of advance, as 
had been supplied in excess during the 
preceding retreat. The assumption of ap- 
proximately equal rates of forward and 
backward shifting of the glacial front is 
justified. Figure 2 is an attempt to por- 
tray graphically the ideas in items 1-3. 

4. During a retreat, the maximum 
volume that caused the channel erosion 
readily carried away the finer constitu- 
ents of the drift being trenched into, 
although a floor of lag cobbles and boul- 
ders gradually accumulated and retarded 
the deepening. When maximum volume 
ceased to arrive, the velocity of the 
shrunken discharge was so decreased 
that the lag material dropped entirely 
out of transit. It thus constituted an 
armor which endured until the climatic 
cycle had come full-turn: until another 
lake stage had been recorded, another 
recessional moraine had been built, an- 
other retreat had begun, excess volume 
was again supplied, the velocity of the 
outlet river had increased, the effective- 
ness of the armor had ceased, channel 
deepening was resumed, and the level of 
ponded water again dropped. 

Because the length of the contributing 
ice front consistently increased during 
the succession of moraines and lakes here 
considered, it is argued that each flood 
stage of Grand River, greater than the 
preceding one, was capable of disposing 
of the boulder armor it found. 

The volume of discharge from such 


melting on the Saginaw, Huron, and: 


Erie lobes as just to balance the rate of 
glacial flow is unknown, the extent and 
rate of retreat unknown, the thickness 
of the marginal belt removed during a 
retreat unknown. The length of the ice 
front whose meltwater escaped down this 
river was from 500 to nearly 600 miles. 
Can reasonable assumptions for the un- 
known factors allow a test of the quanti- 
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tative sufficiency of the proposed mecha- 
nism? 

Taylor was of the opinion that 25 miles 
was a fair estimate for a retreat from 
one of these deployed moraines. His 
figure 6 in the monograph indicates this 
distance between the Port Huron mo- 
raine on the‘ thumb” and the hypotheti- 
cal position of the ice margin out in Lake 
Huron for the following Lake Wayne 
stage, although he also shows a distance 
of 60 miles between these markers along 
the axial line of Saginaw Bay. 

Because the Port Huron retreat 
opened the Syracuse outlet, it was 
greater than retreats from the later Bay 
City or Tawas moraines and may logi- 
cally have been greater than from any 
earlier ones. If 25 miles was its average 
extent, the building of other moraines 
may be assumed to have been followed 
by at least 15 miles of frontal withdrawal. 
Such a retreat would remove all glacial 
ice for from 7,500 to 9,000 Square miles. 
If the thickness of this vanished ice aver- 
aged a fifth of a mile, then between 1,500 
and 1,800 cubic miles of ice became lique- 
fied in each of these frontal back-steps, 
and, because none of the lakes was very 
deep, most of this water disappeared 
down the Grand River outlet channel. 

Another assumption must be made for 
the duration of this excess melting. The 
retreat was across a region of low relief 
and not far from the same latitude as 
north-central Iowa, where Gwynne 
{1942) found faint, close-set, “annual” 
ridgings made during the retreat of the 
Des Moines lobe at the same time that 
some of these Saginaw Basin events were 
occurring. Annual moraines have also 
been found in Quebec (Norman, 1938), 
in Sweden (De Geer, 1940), and in Fin- 
land (Sauramo, 1918). Spacings range 
between 300 and 1,000 feet. Gwynne’s 
average was 15 moraine ridgings per 
mile. 
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Wilson (1943) estimated from incom- 
plete cores of varved sediment at San- 
dusky Bay, south side of Lake Erie, 
that submergence in glacial lake water 
endured here for at least 12,223 years, 
this time beginning with Later Lake 
Maumee and covering but little more 
than the life of all the subsequent glacial 
lake stages south of the Whittlesey out- 
let. During this time the front of the 
Erie lobe retreated approximately 250 
miles, the average rate being 108 feet per 
year or nearly 50 years per mile of re- 
treat. However, if two-thirds of the terri- 
tory abandoned in each individual re- 
treat was retrieved during the next ad- 
vance, then the ice edge itself shifted 
through five times the distance actually 
lost. No matter how many or how few 
the oscillations that occurred during this 
250 miles of total withdrawal of the gla- 
cial front, the actual distance traversed 
in its shiftings was, by our assumption, 


1,250 miles. Hence the actual rate of 
shifting was 1 mile in about ten years, or 
nearly 550 feet per year, a figure which 
agrees well with the rates obtained from 
the annual moraines.‘ 


4 The rate of glacial retreat in New England has 
been computed (Antevs, 1922) from the varve 
record left in numerous small glacial lakes as the ice 
front withdrew through a distance of 185 miles. The 
total time required was 4,100 years. The rates for 
different parts of the distance varied, but the aver- 
age was 24 feet annually. Thus it required nearly 220 
years for 1 mile of average retreat. This computa- 
tion, however, includes time elapsed during what- 
ever readvances occurred but are not in the record. 
If each oscillatory advance repossessed two-thirds 
of the area just abandoned, the edge of the New 
England ice must have melted back and readvanced, 
on the average, at five times the above rate. This 
would be 120 feet annually, or 1 mile in 44 years, a 
considerably slower rate than that in Iowa or the 
Great Lakes. The marked differences in topography 
should be one cause for the difference in rate. 

Sauramo (1923, p. 142), using the varve record 
of southern Finland, obtained retreatal rates of 10.8 
and 13.2 years per mile during his Periods 2 and 3. 
Period 4 saw a greatly accelerated rate because of 
marked decrease in mass of the wasting ice, the mar- 
ginal recession for 75 miles averaging 1 mile in 8 
years. 
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Accepting the interpretation of an- 
nual moraines and using the slower rate, 
a 15-mile retreat would have been con- 
summated in 264 years. This would mean 
that 6} cubic miles of ice, in excess of the 
‘“‘normal” wastage, would be melted an- 
nually during the retreat, which would 
yield more than 29,500 second-feet addi- 
tional discharge to the outlet river. When 
the following advance occurred, the out- 
let river’s ‘“‘normal’’ volume would be 
shrunken by an equivalent amount (al- 
though for a shorter time) because the 
vanished ice must be replaced on a belt 
about two-thirds as wide as that exposed 
previously. Thus, by the figures used 
above, the maximum discharge down 
Grand River exceeded its minimum by 
nearly 60,000 second-feet. This would be 
comparable to adding, during a retreat, 
more than the average flow of the Wa- 
bash River and to subtracting that quan- 
tity during an advance.‘ 

The quantitative sufficiency of the 
mechanism advocated may be challenged 
on another point. Was 264 years long 
enough for the outlet to be deepened 
from 5 to 10 feet? We have no knowledge 
of the volume of meltwater that bal- 
anced glacial flow, hence no estimate for 
the percentage increase or decrease pro- 
duced in Grand River by the glacial oscil- 
lations. But the erosional ability of a 
stream drops in geometrical ratio as 
volume decreases, and a grant of more 
time involves a correspondingly smaller 
volume per year. 

It is conceivable that the recently de- 
posited drift crossed by the outlet river 
farther west still retained much inter- 
stitial water and lacked the degree of 


5A considerably more rapid melting-down of a 
portion of the Scandinavian ice sheet is reported by 
Mannerfelt (1945). If the rate of thinning that he 
computed (3}-4 meters yearly) obtained on the 
Saginaw lobe, the duration of a retreatal stage would 
scarcely exceed 150 years, and the Grand River dis- 


charge would be correspondingly increased. 
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resistance to the river that later compac- 
tion has given it. The erosion that drew 
down the lake levels did occur, and so did 
the pauses in outlet erosion. 


FIELD EVIDENCE SUPPORTING 
THE THESIS 


Monographic treatment may imply, 
but does not guarantee, intensive and 
detailed study. Moreover, the Leverett 
and Taylor investigations were made in 
horse-and-buggy days, without much 
help from road cuts or gravel-pit sections 
and largely without benefit of topo- 
graphic maps. Thus a week’s field work 
in 1950 along the southern side and con- 
vex front of the Saginaw lobe revealed 
several significant items either not men- 
tioned in monograph 53 or at variance 
with descriptions and interpretations 
therein. The four most important items 
are briefly as follows: 

1. A well-developed valley train en- 
ters Imlay River channel from the 
Owosso moraine and is readily traceable 
downstream along Hayworth Creek as 
far as Maple Rapids. It was deposited at 
the Owosso (gth) maximum in an Imlay 
channel of Flint age and became much 
eroded by the last of Imlay discharge 
during Henderson (10th moraine’ time. 

There is no trace of such a deposit in the 
Imlay channel upstream from the city 
of Owosso. 

2. Imlay River never entered Lake 
Saginaw near Duplain, as Taylor thought 
it did (Leverett and Taylor, 1915, pp. 
257-258). Early Lake Saginaw and its 
immediate successor, Earliest Lake Ar- 
kona (= Later Lake Saginaw) discharged 
across the Flint moraine by three con- 
temporaneously functioning spillways: 
Duplain, Ferdun, and Maple Rapids (see 
fig. 1). The Duplain spillway led directly 
into Imlay River, which, throughout 
the entire history shown in fig. 1, dis- 
charged down Hayworth Creek Valley. 
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3. There are five glacial river channels 
unnoted in the monograph but signifi- 
cantly related to those described therein, 
the existence of which modifies and 
makes possible more precise interpreta- 
tions of the region’s glacial history. 

4. An older drift sheet, very bouldery 
in both its till and its washed phases, 
has been reached in the deeper glacial 
and postglacial river channels, and its 
contributions have played a part in the 


later glacial drainage history. 
THE OWOSSO VALLEY TRAIN 

Taylor's idea that the Owosso ice locally 
crowded up on the back slope of the Flint mo- 
raine near Maple Rapids is supported by fea- 
tures of the Flint’s northern slope near Ovid and 
Owosso. For several miles the frontal ridge of 
the Owosso moraine (which ridge Taylor consid- 
ered a part of the Flint) closely conforms to the 
sinuous course of the Flint, and the two are sep- 
arated only by imperfectly developed marginal 
or submarginal channels. This ridge is really 
only a line of partially separated, elongated, 
kamey and eskerine tracts, from whose juxta- 
posed ice a strong flow at Owosso maximum 
crossed the Flint moraine at Ovid, at Lewis 
Drain, and at Owosso. Another kamey, gravel- 
floored discharge from Owosso ice, the Cale- 
donia channel, entered the marginal drainage 
route at Corunna and crossed the Flint by the 
Owosso transection. All told, about 20 miles of 
Owosso ice front appears to have poured its 
meltwater across the Flint at these three places. 
At no other frontal stand along the Saginaw 
lobe’s southern limb in this latitude is there rec- 
ord of any comparable meltwater discharge or 
load of gravel. 

Along the front of the Flint flowed the 
early Imlay River. Its volume and gradient were 
inadequate for carrying away all the load con- 
tributed by the Owosso discharge. The valley 
train that resulted has been traced from Owosso 
to Maple Rapids, 30 miles farther west, its sur- 
face descending about 25 feet in that distance. 
Numerous large gravel pits show a deposit at 
least 30 feet in maximum thickness. The mate- 
rial is well washed and well sorted, has many 
strata of westward-dipping foresets, and is pro- 
gressively finer-textured with increasing dis- 
tance from the moraine. It is to be sharply dis- 
tinguished from contiguous morainal gravel de- 
posits. 
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When ice withdrew from the Owosso mo- 
raine, the valley train was promptly attacked by 
the now underloaded Imlay River, leaving only 
terraces and mid-channel hillocks as the glacial 
riverway was re-excavated. Only the time be- 
tween Owosso and Henderson maxima is allow- 
able for this destruction because Imlay River, 
finding a lower spillway northwestward from 
Flint, Michigan, wnen the Henderson ice front 
retreated, entered Lake Saginaw at Flushing 
and abandoned the entire channel length that 
we are considering in this study. 


IMLAY RIVER’S RELATIONS TO 
LAKE SAGINAW 


Taylor’s idea that Imlay River ever entered 
Lake Saginaw before the above-noted spillway 
was uncovered is in error but, without today’s 
maps and excavations, would be almost impos- 
sible to disprove. Taylor located his earlier de- 
boucnment in the Duplain gap. Here the well- 
marked Imlay channel, now followed by Maple 
River, virtually disappears as it encounters the 
broad embayment of the low, sandy lake plain 
flanked by Lake Saginaw beaches. Here is the 
head of Hayworth channel, correctly interpret- 
ed by Taylor as a part of the early Imlay River 
course, its floor only 5 feet higher than that of 
most of the embayment bottom but 15 feet 
lower than a beach ridge margining the bottom 
land only a mile distant. Here tne Owosso valley 
train continues from the Maple-Imlay channel 
westward across the head of the embayment and 
on down the Hayworth-Imlay channel. Not a 
pebble from the 30 feet of gravel here was ever 
carried into Lake Saginaw, either when the fill 
was made or later when 20 feet or more of dis- 
section by Lake Maumee discharge occurred.‘ 

Owosso ice covered much of the Duplain sag 
when the valley train was deposited, and Imlay 
River had no choice but to follow tne Hayworth 
channel of Flint moraine time. But only ponded 


© The reconnaissance character of Leverett and 
Taylor’s study in this region is shown by the exist- 
ence of an esker associated with the Ferdun channel, 
higher, longer, and more massive than any described 
under “Moraines of the Eastern Limb of the Saginaw 
Lobe” (Leverett and Taylor, 1915, pp. 249-250). It 
is visible for 3 miles in any direction. The failure to 
recognize its existence, or that of the channel, in the 
monograph is a further indication that neither 
author ever traversed the township in which these 
features are subcentrally located. Both esker and 
channel show on the Perrinton topographic map 
(surveyed 1915-1916). Leverett’s 1924 depiction of 
a narrow lake embayment here seems obviously only 
a map interpretation. 
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water stood here during the dissection, and only 
southward discharge from the head of the em- 
bayment can explain the absence of a deltaic 
deposit in it. 

A remnant of the Owosso valley train, cover- 
ing about 2 square miles, lies on the south side 
of Imlay channel at the head of Hayworth and 
the junction of Duplain and Maple-Imlay. 
Taylor was puzzled by its features, describing it 
(Leverett and Taylor, 1915, p. 257) as “more 
like a kame than any otner type—a kame depos- 
ited in a small re-entrant of the [Flint] ice front 
and shaped partly by the ice with which it was 
in contact and perhaps partly by the later [Im- 
lay] river current.” Today the 5-foot contouring 
of the Perrinton and Elsie topographic maps 
(see fig. 3) and the existence of two large gravel 
pits make clear that the lower southern part of 
the triangular area is a shallow channel of the 
aggradation episode. Its two northerly margins 
are walls of deeper channels, one from Maple- 
Imlay and one from Lake Saginaw along the 
Duplain sag. The highest part, in the apex of the 
triangle and along the western side, Taylor de- 
scribed as “‘a prominent gravel ridge—like a 
parapet—{which] closely resembles a beach 
ridge.” It is, in fact, two gravel beaches, the high- 
er at 735+ feet A.T., built of the local gravel 
by Lake Saginaw waves after the episode of val- 
ley-train dissection was completed, i.e., in early 
post-Henderson time. 

Another conspicuous remnant of the Owosso 
valley train is a mid-channel hill about a mile 
long, separating Maple and Little Maple rivers 
where both are using the Imlay channel west of 
Ovid. The summit of the hill is about 25 feet 
above the late Imlay floor on the north side.? 


PREVIOUSLY UNNOTED GLACIAL 
RIVER CHANNELS 


The glacial river channels involved in this 
problem and recognized in monograph 53 are 
the Grand, the Lookingglass-Libhart, the 


7A medial gravel ridge occupies the last mile of 
Maple-Imlay south and southwest of Duplain, with 
the northwest-flowing Maple River on one side and 
a shallower, southeast-draining depression on the 
other. If the ridge is a mid-channel remnant of the 
valley train, the flanking trenches should both open 
out into the Duplain sag. But the shallower one is 
blocked in the village by a transverse sand ridge 
whose crest is above 735 feet A.T. It is the original 
“Duplain beach” of Taylor (1897), made by the 
highest recorded stage of Lake Saginaw. The blocked 
channel seems to be a survival from channel braids 
on the valley-train surface, similar to the one south 
of the “parapet” beach. 


| 


254 


Maple Rapids, the Stony, the Imlay (Maple and 
Hayworth), the Duplain, and the Corunna- 
Owosso. The present interpretation of the Du- 
plain reverses its supposed direction of flow. The 
Corunna-Owosso channel is reinterpreted as a 
discharge route from ice on the Owosso moraine, 
not a detour or accessory channel of Imlay 
River. 

The Swargard channel was added on Lever- 
ett’s 1924 map. The writer here adds the Cook, 
the Ferdun, the Lewis, the Baker, and the 
Caledonia channels, all a part of the shifting 
drainage changes in front of the shrinking Sagi- 
naw lobe. 
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Lyons. The modern Grand now follows that 
course. 

Retreat from tne Lyons (6th) moraine al- 
lowed the Cook channel to form, the swampy 
col in it but little above 710 feet A.T. and the 
actual channel floor probably still lower beneath 
the swamp accumulation. Cook was but a small 
distributary of Stony channel. 

Stony channel carried marginal drainage 
along several miles of the Fowler (7th) moraine 
front. A low place between tne Lyons and Fow- 
ler moraines in the vicinity of Pewamo was 
aggraded to make a broad gravel flat at 720 feet 
A.T. Then, the load apparently decreasing, the 
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Fic. 3.—Vicinity of Duplain. From Elsie and Perrinton topographic maps, U.S. Geological Survey 


During the building of the Portland (5th) 
moraine, marginal drainage from tne southern 
side of the Saginaw lobe used an imperfectly de- 
veloped cnannel now followed by Lookingglass 
River as far as Portland, thence by modern 
Grand River for 3 miles, and thence by Libhart 
Creek to its entrance into the Glacial Grand. 
Tne floor of Libhart channel and terraces at and 
near junction with Glacial Grand (Portland 
topographic map) record a main channel floor 
but little below 740 feet A.T. at that time. 

Retreat from the Portland moraine exposed a 
surface lower than the Libnart floor, and conse- 
quently the Glacial Lookingglass crossed to the 
back of the moraine aLwut 4 miles south of 


Glacial Stony trencned into its own deposits. 
The Cook distributary probably was still func- 
tioning, for its floor is lower tnan the Pewamo 
flats. 

Glacial drainage. along tne full length of 
Stony channel ceased on retreat from the Fow- 
ler moraine. Despite a distributary-like relation 
in ground plan, Stony almost surely never was 
supplied from Imlay River. lhe eastern few 
miles of its length are very weakly developed, 
compared with the rest of its length, and at its 
junction with the Imlay its floor hangs about 20 
feet above that of the main channel. Except for 
later local erosion, there is nardly a notci.ing of 
the top of tne Imlay bluff at the junction. On the 
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monograph map, plate VII, Stony was shown 
as the lower part of a channel 50 miles long, 
called the “Bennington.”’ The 1924 map, issued 
after two of the topographic maps covering this 
Bennington channel had appeared and surveys 
had been initiated on two others, omitted 35 
miles of this supposed channel and depicted 
Stony correctly, as in the accompanying figure 
1. This change appears to be another revision 
from map study. 

Cook and Stony are of interest here because 
they record a deepening in the trunk drainage 
line, the Grand River channel, of at least 20 feet 
from the time of the Portland (5th) maximum to 
that of the Fowler (7th). 

The Swargard channel records a broad, shal- 
low overflow across the Fowler from the earliest 
glacial drainage along the front of the Flint 
(8th) moraine. Swargard may have been an 
Imlay distributary before the trenching which 
followed accession of Lake Maumee’s overflow. 
Although Swargard’s junction with Stony is ac- 
cordant, the trenching in Pewamo flats and 
Cook channel seems best explained as done by 
Swargard discharge. So interpreted, the Grand 
River channel floor had been lowered another 
10 feet by Flint moraine time. 

Imlay’s course in the Saginaw Basin was 
along an intermorainic depression into which 
the river trenched to maximum known depths 
of 30 feet. Following this came the deposition of 
the Owosso valley train, nearly filling this trench 
and spreading a broad sheet of gravel in lower 
Hayworth, where lay a deeper portion of the 
original depression. Flat-topped remnants of 
this valley train, left from the immediately post- 
Owosso dissection, are 720-725 feet A.T. at Ma- 
ple Rapids, and the floor of Glacial Grand at this 
place thus was approximately 715 feet A.T. 
when the Owosso moraine was deposited. 

Imlay River may not have deepened the 
Glacial Grand channel very much before the 
making of the Owosso valley train, for all the 
debris from trenching along the front of the 
Flint had to be carried across the state to the 
Allendale delta in Lake Chicago, west of Grand 
Rapids. 

Taylor depicted a broad gap in the Flint 
moraine at Maple Rapids, into which the Owos- 
so ice pushed and built its own moraine. He 
believed the gap to have been eroded by the 
outlet of a pre-Owosso lake, across the newly 
built Flint moraine. This gap, however, is at 
least four times as wide as that made in the 
Owosso moraine at the same place by all subse- 
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quent outlet erosion. It seems much more 
probable that the pre-Owosso gap was largely 
a sag in the crest of the Flint moraine, a dis- 
continuity similar to the Duplain sag in both 
Flint and Owosso moraines and to others in vir- 
tually every moraine the region possesses.® 

Early Lake Saginaw, the first of the lake 
series in the Saginaw Basin to leave positive 
records, was dammed by the Owosso moraine. 
Because some of its beaches lie on the back slope 
of the Henderson moraine at the same altitude, 
735+ feet A.T.,° as at Duplain, this highest 
recorded lake stage must have endured without 
lowering by outlet erosion during Owosso re- 
treat, Henderson advance, and much, if not all, 
of Henderson retreat. Its waves made the 
“parapet” beach ridge on the triangular rem- 
nant of the Owosso valley train 3 miles west of 
Duplain; hence dissection of the valley train 
was already advanced. 

This duration through two retreats is a seem- 
ing contradiction of the thesis that retreating 
glacial episodes were outlet-deepening, and 
therefore lake-lowering, times. The contradic- 
tion disappears when one remembers that the 
Owosso valley train had nearly filled the early 
Imlay channel and had extended an unknown 
distance down the Glacial Grand. ‘The energies 
of late Imlay River and the contemporaneous 
Glacial Grand were all absorbed in the immedi- 
ate task of restoring the earlier channel depths. 

Ferdun channel’s southern part was origi- 
nally an esker trough across the Flint moraine. 
As a surface route across both Flint and Owosso 
moraines, it was one of three dischargeways for 
Lake Saginaw, Maple Rapids and Duplain be- 
ing contemporaries. Because its floor across the 
Flint is accordant with that of the late Imlay 
trench cut into the valley train in Hayworth, 
Ferdun is known to have functioned at least to 


* One has only to compare the map in the mono- 
graph with the revision issued in 1924 to see what 
difficulties in moraine mapping these discontinuities 
constitute. At the front and along the south limb of 
the Saginaw lobe, Leverett in 1924 abandoned more 
than 100 miles of moraine lengths that were shown 
in 1915, bodily shifted nearly 100 more miles, and 
added go miles of newly depicted linear moraines! 
The moraine least altered was the strong Flint, with 
the Imlay channel hard against the foot of its south- 
ern slope. 

* Altitudes are of beach crests, taken from the 

hic maps. Taylor’s figure was “720 to 725 
feet aneroid,” doubtless measured at the lakeward 
base of a ridge. The maps show that 730 is a better 
figure. 
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the end of that episode of dissection, at which 
time Lake Saginaw had been lowered to 710- 
715 feet A.T. (=Highest Arkona) and Imlay 
had found the Flint-Flushing short cut. 

At this same Later Lake Saginaw time, Du- 
plain discharge was recorded by the 715-foot 
col at the head of Hayworth, essentially the 
channel in the valley train which Imlay aban- 
doned for the route north to Flushing. 

Part of the Lewis channel transects the Flint 
moraine and part appears to be an esker trough 
in the Owosso. Lewis contributed considerable 
gravel to Imlay during the valley-train aggrada- 
tion. Glacial drainage continued to use it dur- 
ing Henderson maximum, when the valley train 
was still being trenched and Lake Saginaw 
still stood at its highest level. 

The Baker channel crosses an Owosso mo- 
raine unit along what originally must have been 
another subglacial drainage route. It appears 
to have been supplied later from Henderson 
meltwater and to have found Lake Saginaw 
still at 725-730 feet A.T. 

From this examination of sequences among 
glacial drainage channels of the region, we con- 
clude that the earliest ponded water back of the 
Owosso moraine had three outlets—Duplain, 
Ferdun, and Maple Rapids—to the same trunk 
stream—Glacial Grand River—and that, up to 
the time of making the Later Lake Saginaw 
beach, none obtained precedence over the 
others. Only after te Owosso valley train had 
been eroded could Grand River, fed by lake 
water alone, begin deepening when a glacial 
retreat was in progress and offer opportunity 
for competition among the three. This was dur- 
ing the West Haven retreat, next after the Hen- 
derson withdrawal. But the three ran neck and 
neck in that first deepening. Then Duplain and 
Ferdun fell out of the contest, and only Maple 
Rapids outlet continued to be deepened when 
the Chesaning retreat occurred and the ponded 
water fell to the intermediate Arkona level. 


THE SUBJACENT BOULDERY DRIFT 


Remnants of a channel-bottom boulder ar- 
mor, suggesting that which figures in the writ- 
er’s proposed mechanism, were reported by Tay- 
lor to lie on “terraces and boulder benches” 
from the general vicinity of Maple Rapids as 
far down the Grand River outlet channel as 
Lyons. Referring to the entire deepening, Tay- 
lor said that “with the floor of the channel heav- 
ily protected by boulders, the lowering was 
slow and required a long time.” But it was also 
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spasmodic, with three times as many episodes 
when the armor was ineffective as Taylor’s ex- 
planation (stoping initiated by Lake Chicago 
lowerings) can account for. 

Because the supply of boulders adequate for 
checking channel erosion during a glacial ad- 
vance is assumed in this study to have been 
swept away by the enlarged and quickened Gla- 
cial Grand River during the next retreatal epi- 
sode, the boulder armor now lying on the late 
Warren channel floor (Leverett and Taylor, 
1915, fig. 2) is not the cumulative result of the 
earlier episodes of deepening. Its impressiveness 
is due largely to contributions from a subja- 
cent bouldery drift that was reached only in this 
latest channel deepening. The channel floors 
along Stony and Hayworth dischargeways, 
although earlier developed, shorter lived, and 
less deeply trenched, show the same extreme 
boulderiness in their downstream portions but, 
traced upstream, possess only boulderless and 
almost cobbleless gravel. 

Taken with the exposures in ravines and pits 
near Maple Rapids of bouldery till and outwash 
beneath the moraine drift and its gravel, this 
seems good evidence that these excessively 
bouldery channel floors are not cumulative ag- 
gregates. The high percentage of glacially 
marked boulders in these exposures also argues 
for only a brief experience at the bottom of the 
outlet river. 


LACK OF SHORE LINES ON THE PORT HURON 
MORAINE FRONT 

The glacial front at which the Port 
Huron moraine was built served also as 
the dam for Whittlesey and Late Arkona 
lakes. But shore lines recording these 
contemporaneous lake levels have not 
been found on this moraine. Taylor 
reported both shores on older moraines 
in several places where they, with the 
Port Huron, outlined narrow bays at 
these lake levels. He reasoned that their 
strong development in such constricted 
situations was impossible and that the 
open expanses of water required for ade- 
quate waves existed during the immedi- 
ately pre-Port Huron retreat and were 
destroyed as the Port Huron advance 
culminated. 
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GLACIAL LAKE STAGES IN SAGINAW BASIN, MICHIGAN 


However, some other explanation is 
needed for the lack of these shores in 
places where the Port Huron moraine 
faced wide stretches of lake water. The 
thesis of this paper supplies a satisfac- 
tory explanation. It requires that a static 
lake level exist only during a glacial ad- 
vance and that this level start dropping 
immediately when a retreat begins. No 
moraine fronts, therefore, should carry 
more than the faintest shore lines of their 
contemporaneous lakes. 


SUMMARY 


Glacial Grand River, flowing west 
across the southern peninsula of Michi- 
gan, carried meltwater directly from the 
Saginaw glacial lobe during the building 
of the Lake Border and Port Huron mo- 
raine systems and from the successive 
glacial lakes impounded in front of this 
lobe. It also received the overflow of 
several successive lakes dammed by the 
Erie and Huron lobes. 

Because the named lakes were largely 
a succession of lowered levels in the same 
basins, discharging through the same 
outlet, some mechanism must have oper- 
ated to produce the sequence of static 
levels alternating with lowering levels. 

The volume of this great glacial river 
varied between maximal discharge, when 
the contributing lobe fronts were retreat- 
ing, and minimal discharge, when their 
fronts readvanced. Because the succes- 
sive lowerings resulted directly from out- 
let deepening, the mechanism proposed 
is (1) that the greater volume and ve- 
locity of the river during a glacial retreat 
produced the episodes of lowering levels 
and (2) that during a readvance the 
weakened river, hampered by a bottom 
armor of lag boulders and unable to 
deepen its channel, allowed the static 
levels time for making shore-line records. 
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The pondings began when the shrink- 
ing Saginaw lobe had retreated east of 
the Michigan-Huron divide. Each mo- 
raine thereafter deposited should theo- 
retically be matched by a lake stage. 
But the earliest of these lake levels, 
Early Lake Saginaw, endured through 
the retreat from the Owosso moraine, 
the advance to the Henderson moraine, 
the retreat from that moraine, and (theo- 
retically) the advance to the West Haven 
moraine. 

This departure from the theoretical 
schedule resulted from a literal clogging 
of the outlet river by a considerable val- 
ley train of Owosso age, brought into 
Glacial Grand River below the lake out- 
lets. Low places in the Owosso moraine 
provided three outlets, all continuing to 
function until the valley train had been 
sufficiently cleared away. After that, 
deepening in one outlet outdistanced the 
other two, and all subsequent lake dis- 
charge from the Saginaw Basin used and 
continued to deepen the Maple Rapids 
outlet channel. 

According to the scheme proposed, 
Early Lake Saginaw endured from the 
retreat from the Owosso moraine to the 
building of the West Haven moraine. 
Early Lake Arkona (=Later Lake 
Saginaw) endured during the advance 
which culminated in the Chesaning mo- 
raine. Intermediate Arkona existed while 
the lobe front was advancing to the posi- 
tion of the moraine north of the Chesa- 
ning. Late Arkona’s existence coincided 
with the advance to and the building of 
the Port Huron moraine. Early Lake 
Warren levels were recorded while the 
ice front was advancing to the Bay City 
moraine’s position. Intermediate Warren 
endured through the episode of advance 
to the Tawas moraine. Lowering to the 
Late Warren stage coincided with retreat 
from the Tawas, the last moraine left on 
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land in extreme eastern Michigan. To 
maintain that lowest Warren level while 
its shore lines were made, theory requires 
another glacial advance. The moraine 
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marking its culmination is assumed to 
have been entirely water-laid and still to 
be submerged in Lake Huron and Geor- 
gian Bay. 
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LIQUIDUS TEMPERATURES IN MIXTURES OF THE 
FELDSPARS OF SODA, POTASH, AND LIME' 


R. R. FRANCO? AND J. F. SCHAIRER 
Geophysical Laboratory, Carnegie Institution of Washington 


ABSTRACT 


Liquidus data are reported for a series of forty-two ternary feldspar compositions, including com 


itions 


at each ro per cent interval by weight of soda, potash, and lime feldspars. Only two fields, those of feld 
and leucite, appear on the liquidus surface. In addition to the liquidus data, the temperature of appearance 
of a second solid phase was determined for some of the compositions. Measurements of the index of refrac- 
tion of the feldspar glasses are recorded. These data are principally of value as a starting point for studies 
of these same compositions with water as an additional component. 


INTRODUCTION 


The various feldspars are among the 
most common rock-forming minerals. 
Naturally, petrologists and mineralogists 
have devoted considerable time to studies 
of the properties, occurrence, and inter- 
relations of feldspars, particularly those 

‘of soda, potash, and lime. One of the first 
studies made by Day and Allen (1905), 
at the U.S. Geological Survey just pre- 
vious to the founding of the Geophysical 
Laboratory, was concerned with the 
melting relations of the feldspars, and 
today, even though nearly a half-century 
has elapsed and many studies have been 
made, there is still much to be learned 
about their stability relations. 


THE INDIVIDUAL FELDSPARS 


The melting relations of each of the 
feldspars of lime, soda, and potash have 
been determined. Day and Allen (1905) 
found that lime feldspar melted congru- 
ently. Rankin and Wright (1915) in a 
study of the system CaOQ-Al,O,-SiO, 
showed its melting relations in the pres- 
ence of excess of each of its components. 
Greig and Barth (1938) determined the 
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congruent melting temperature of soda 
feldspar and showed its melting relations 
with the feldspathoid nepheline. Schairer 
and Bowen (19474), in a study of the 
system Na,O-Al,O,-SiO,, showed the 
melting relations of soda feldspar in the 
presence of excess of each of its compo- 
nents. Morey and Bowen (1922) estab- 
lished the incongruent nature of the 
melting of potash feldspar to crystals of 
leucite and a liquid more siliceous than 
feldspar. Schairer and Bowen (19472), 
in a study of the system K,O-Al,O,-SiO,, 
showed the melting relations of potash 
feldspar in the presence of excess of each 
of its components. 


PAIRS OF FELDSPARS 


The melting relations of each pair of 
the feldspars under consideration have 
been determined. 


SODA FELDSPAR-LIME FELDSPAR 


Bowen (1913) studied the plagioclase 
series and found that at high tempera- 
tures soda and lime feldspars form a com- 
plete and unbroken series of solid solu- 
tions with a continuous rise in the liqui- 
dus and solidus from pure soda feldspar 
to pure lime feldspar. Schairer (unpub- 
lished) in a study of the system plagio- 
clase-silica prepared feldspar composi- 
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tions at each 10 per cent interval by 
weight between soda feldspar and lime 
feldspar and one additional feldspar 
(soda feldspar 95, lime feldspar 5). Liq- 
uidus data on these agree perfectly with 
the data of Bowen on the plagioclases, 
and all data at and above the tempera- 
tures of beginning of melting indicate a 
complete and unbroken series of solid 
solutions. Tuttle and Bowen (1950) 
found two different crystalline modifica- 
tions of soda feldspar with an inversion 
at about 700° C. and a probable break in 
the plagioclase series at temperatures 
well below that at which a liquid phase 
appears in the dry binary system. They 
are continuing their studies of these feld- 
spar relations. 


SODA FELDSPAR~-POTASH FELDSPAR 


Schairer and Bowen (1935) published 
a preliminary paper on the system 
NaAlSiO,-KAISiO,-SiO, and gave a pre- 
liminary diagram for the join NaAlSi,Os- 
KAISi,Os. Recently Schairer (1950) pub- 
lished the data for this join and the re- 
vised diagram. There is a complete series 
of solid solutions between these two feld- 
spars, with a minimum on the melting 
and freezing curves and a further com- 
plication that all these feldspars with 
more than 49 per cent potash feldspar 
melt incongruently to leucite and liquid. 
A short time ago Bowen and Tuttle (1950) 
published on the system NaAlSi,Os¢- 
KAISi,Os-H,O. They showed (1) that the 
complete series of solid solutions between 
these two feldspars persists in the pres- 
ence of water at low and moderate pres- 
sures; (2) that, as the pressure of water is 
increased, the incongruent melting of 
potash feldspar and potash-rich feld- 
spars is eliminated; (3) that the position 
of the minimum on the melting and 
freezing curves shifts from 65NaAlSi,Ox 
35KAISi,Os in the dry system toward 
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compositions richer in soda feldspar; and 
(4) that thesolid solutionsof K AlSi,Osand 
NaAlSi,Os, although complete at higher 
temperatures, are only partial at lower 
temperatures and the feldspars prepared 
at higher temperatures unmix or, con- 
versely, a pair of feldspars (one soda-rich 
and the other potash-rich) stable at lower 
temperatures become one homogeneous 
feldspar at the higher temperatures. Ad- 
ditional relations between the low-tem- 
perature and high-temperature crystal- 
line forms of these pure feldspars and 
solid solutions are now being studied. 


LIME FELDSPAR~—POTASH FELDSPAR 


Schairer and Bowen (19476), in a 
study of the system leucite-anorthite- 
silica, determined the relations between 
lime feldspar and potash feldspar. All 
compositions between lime feldspar and 
potash feldspar become completely crys- 
talline only at and below 1,040° + 20° 
and then consist of a mixture of pure 
(or very nearly pure) lime feldspar and 
potash feldspar. At temperatures above 
1,040° + 20°, leucite appears as one of 
the solid phases because of the incon- 
gruent nature of the melting of potash 
feldspar. 


THREE-COMPONENT FELDSPARS 


We have prepared a series of synthetic 
feldspar compositions at 10 per cent in- 
tervals by weight of soda, potash, and 
lime feldspars with a few additional com- 
positions. All these compositions have 
been crystallized completely or nearly 
completely. We have completed and 
record here our determinations on the 
liquidus surface. Studies of the beginning 
of melting are still in progress. All these 
feldspar compositions will be studied 
with water as an additional component 
at a series of isobars. 
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LIQUIDUS TEMPERATURES OF FELDSPARS 


PREPARATION OF COMPOSITIONS 


In the preparation of silicate composi- 
tions containing alkalies, special care 
must be exercised because of volatiliza- 
tion of the alkalies and consequent de- 
parture from composition. Particular 
care is necessary when both Na,O and 
K.O are present. In the preparation of 
the feldspar compositions for this study, 
two separate glasses were made with the 
proportions Na,0-6SiO, and K,0-6Si0,. 
These were prepared by melting Na,CO, 
or KHCO,, respectively, with purified 
quartz in platinum crucibles. The loss of 
alkalies was minimized by first sintering 
at low temperatures (around 700°) and 
then gradually raising the temperature 
to about 1,300° or 1,400°, respectively. 
The small loss (nearly zero in the case of 
the Na,O-6SiO, and usually less than 
0.0100 g. in the case of K,0-6SiO, on 
20.0000-g. samples) was determined by 
weighing and was corrected by adding 
the appropriate amount of Na,CO, or 
KHCO,, respectively, before subsequent 
fusions. At least three fusions with inter- 
mediate crushings of the quenched 
glasses were necessary to obtain a ho- 
mogeneous product, which was then crys- 
tallized completely to obtain a product 
which remains dry in a desiccator over 
fused KOH sticks and is easy to weigh. 

Calculated amounts of “‘Na,0-6SiO,,”’ 
“K.0-6Si0,,”" pure CaCO,, pure AlL,O,, 
and purified quartz were weighed in- 
to platinum crucibles and fused to pre- 
pare the feldspar compositions. Depend- 
ing on the compositions and consequent 
viscosity of the glasses, at least four fu- 
sions and sometimes as many as eight, 
with intermediate crushing of the 
quenched glasses, were necessary to 
effect complete solution of the Al,O, and 
secure a homogeneous glass. The ho- 
mogeneity of each of the glasses was 
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carefully checked under the petrographic 
microscope by obtaining the index of re- 
fraction of the glass grains. Each of the 
feldspar glasses was crystallized com- 
pletely or nearly completely at appropri- 
ate temperatures (1,000° for most of 
them) before proceeding with quenching 
experiments. 


DETERMINATION OF LIQUIDUS 
TEMPERATURES 


A determination of the liquidus tem- 
perature was made for each of the forty- 
two ternary feldspar compositions pre- 
pared, and for some of the compositions 
the temperature of appearance of a sec- 
ond solid phase was determined. These 
studies were made by the method of 
quenching (Shepherd, Rankin, and 
Wright, 1909), which consists in holding 
a small charge of a given composition at 
a measured temperature for a period of 
time adequate for the attainment of 
equilibrium between crystals and liquid, 
and then chilling instantly to room tem- 
perature. Quenched charges were ex- 
amined with a petrographic microscope 
to determine the nature and number of 
crystalline and liquid phases, and crys- 
tals were identified by their optical prop- 
erties. Liquids quench to a glass. 

Equilibrium between crystals and 
liquid in the ternary feldspar composi- 
tions was reached in a few hours at tem- 
peratures near and above 1,450°, in 1 
day at temperatures between 1,350° and 
1,450, and in 2 days at temperatures be- 
tween 1,200° and 1,350°. Below 1,200° 
equilibrium between crystals and liquid 
was attained only after a week or longer, 
depending on the particular composition. 

Quenching experiments were con- 
ducted in platinum-wound electric fur- 
naces whose temperature was controlled 
(+ 2°) with a temperature regulator— 
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TABLE 1 
RESULTS OF QUENCHING EXPERIMENTS 


COMPOSITION IN 
Per Cent 


| 


| 
| 
Time 
| 


KAISi,Os “asin CaALSi,Os 


Points in the Feldspar Field 


1,510 3 hr. Very rare feldspar in glass 

1,515 3 hr. Lote ef teldspa 

1,490 E ts of fe r in glass 

1,495 ; All glass 

1,490 : Small amount feldspar in glass 

1,405 All glass 

1,470 ; Moderate amount feldspar in glass 

1,475 All glass 

1,470 ; Moderate amount feldspar in glass 

1,475 ‘ All glass 

1,465 , Lots of feldspar in glass 

1,470 4 All glass 

1,440 : Moderate amount feldspar in glass 

1,445 x All glass 

1,445 5 Rare feldspar in glass 

1,450 All glass 

1,440 : Moderate amount feldspar in glass 

1,445 All glass 

1,440 ‘ Moderate amount feldspar in glass 

1,445 ; All glass 

1,415 2 Rare feldspar in glass 

1,420 ’ All glass 

1,415 : Moderate amount feldspar in glass 

1,420 All glass 

1,415 , Very rare feldspar in glass 

1,420 All glass 

1,410 ‘ Moderate amount feldspar in glass 

1,415 , All glass 

1,405 3 Moderate amount feldspar in glass 

1,410 ; All glass 

1,380 ’ Rare feldspar in glass 

1,385 All glass 

1,380 1 Small amount feldspar in glass 

1,385 : All glass 

1,380 ’ Very rare feldspar in glass 

1,385 All glass 

1,370 : Small amount feldspar in glass 

1,375 All glass 

1,205 Lots of feldspar and rare leucite in glass 

1,210 Lots of feldspar in glass 

1,270 ; Lots of feldspar and small amount leu- 

cite in glass 

1,275 , Lots of feldspar in glass 

1,365 ; Rare feldspar in glass 

1,370 3 All glass 

1,315 3 Lots of feldspar and rare leucite in glass 

1,320 ; Lots of feldspar in glass 

1,365 ; Very rare feldspar in glass 

1,370 All glass 

1,335 Rare feldspar in glass 

1,340 All glass 

1,330 Small amount feldspar in glass 

1,335 ; All glass 

1,330 Rare feldspar in glass 

1,335 x All glass 

1,315 , Moderate amount feldspar in glass 

1,320 All glass 

1,220 2 Lots of feldspar and very rare leucite in 
lass 


g 
1,225 Lots of feldspar in glass 
1,310 ; Rare feldspar in glass 
1,315 : All glass 


tive Temp. Pu 
or GLASS | (CC) ASES 
10 | 80 | 
| | 
10 20, 7° 
| 
20 | 10 7° | 
ee 10 30 | 60 
g 1.536.... 20 | 20 60 
| 1.536... 30 10 60 
10 | 40 5° 
20 30 50 
1.531.. | 30 20, 50 
40 | 10 50 
20 40 40 
30 30 40 
40 20 40° 
5.506... 50 10 40 
10 60 30 
20 50 30 
{ 1.512.. 30 eo | 3 
40 30 30 
{ 
50 20 | 30 
60 10 | Rie) 
| | 
1. 505 10 7° 20 
1.504.... 20 | 60 20 
30 50 20 
40 40 20 
| 
5° 30 20 


TABLE 1—Continued 


CoMPOSITION IN 


Weicat Per Cent 


T 
KAISi,Os NaAlSi,Os | 


Points in the Feldspar Field —Continued 


1,255 48 hr. Lots of feldspar and rare leucite in glass 
1,260 24 hr. Lots of feldspar in glass 

1,305 ’ Moderate amount feldspar in glass 
1,310 . All glass 

1,265 " Small amount feldspar in glass 

1,270 All glass 

1,255 ; Lots of feldspar and rare leucite in glass 
1,260 ’ Lots of feldspar in glass 

1,240 Small amount feldspar in glass 

1,245 All glass 

1,235 Rare feldspar in glass 

1,240 a All glass 

1,230 : Small amount feldspar in glass 

1,235 All glass 

1,225 ’ Rare feldspar in glass 

1,230 All glass 

1,220 ’ Small amount feldspar in glass 

1,225 J All glass 

1,215 ‘ Smal amount feldspar in glass 

1,220 I All glass 

1,185 Lots of feldspar and rare leucite in glass 
1,190 Lots of feldspar in glass 

1,100 Small amount feldspar in glass 

1, 105 All glass 


Points in the Leucite Field 


1,305 24hr. | Small amount feld and small 

| amount leucite in g 

1,310 24 hr. | Small amount leucite in glass 

1,315 24 hr. All glass 

1,335 24 hr. Very rare leucite in glass 

1,340 24 hr. All glass 

1,315 24 hr. Moderate amount leucite in glass 

1,310 4 days | Moderate amount leucite and moderate 

amount feldspar in glass 

1,265 48 hr. Small amount leucite and small amount 

feldspar in glass 

1,270 48 hr. Rare leucite in glass 

1,275 48 hr. All glass 

1,210 72 hr. Moderate amount leucite and very rare 

feldspar in glass 

1,215 72 hr. All glass 

1,295 24 hr. Small amount leucite in glass 

1,300 6 hr. All glass 

1,225 48hr. Lots of leucite and small amount feld- 
spar in glass 

1,230 48hr. Lots of leucite in glass 

1,380 2ohr. | Small amount leucite in glass 

1,385 2ohr. Ali glass 

1,245 72hr. | Lots of leucite and very rare feldspar in 

lass 


re. 
1,250 48hr. — Lots of leucite in glass 
| 


Point on Leucite-Feldspar Boundary Surface 


1,145 72 hr. Rare leucite and rare feldspar in glass 
1,150 72 hr. All glass 


| 
rive 
or Giass |~ 
+ 0.003 
60 20 20 
1.500... ‘| 62 23 15 
10 80 10 
1.498.... 20 7° 10 
30 60 10 
1.498.... 4° 5° 10 
5.998.... 50 4° 10 
1.407...- 57 33 10 
1.493----| 47 3 
} 
1.504.... 67 | | 
1.504.... 7o | 10 | 20 
| | 
65 | 20 15 
60 | 30 10 | 
7° 20 10 | 
| | | | 
| 10 10 
| 
| | 
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the Geophysical Laboratory furnace 
thermostat (Roberts, 1925)—-somewhat 
modified from the original design but 
identical in principle. Temperature was 
measured with a thermocouple placed 
nearly in contact with (about 1 mm. dis- 
tant from) the small charge wrapped in a 
tiny platinum envelope, which is sus- 
pended in that part of the furnace deter- 
mined as the “thot point.”’ In a well-in- 
sulated furnace there is a range of about 
5 mm. on either side of the “thot point”’ 
where the temperature falls off about 
one-half of a degree. The small platinum 
envelope containing the charge hangs en- 
tirely within this zone. 
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The Pt-PtgoRhio thermocouple was 
calibrated frequently at several fixed 
points defined in degrees centigrade as 
follows: 


Palladium melting point.... 1,549°5 
Diopside (CaO+MgO- 2Si0,) 
melting point............ 


Gold melting point......... 


1,301° 5 
1,062°6 


RESULTS OF QUENCHING EXPERIMENTS 


The results of the quenching experi- 
ments are given in table 1 and are shown 
graphically in figure 1. Only two fields— 
that of feldspar (solid solutions between 
soda, potash, and lime feldspars) and 
that of leucite (pure K,0-AlL0,-4SiO, 
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Fic. 1.—Equilibrium diagram for the system soda feldspar—potash feldspar—lime feldspar 
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LIQUIDUS TEMPERATURES OF FELDSPARS 


or KAISi,O¢)—-were found on the liquidus 
surface. No data were obtained on the 
composition of the feldspar crystals at 
or below the liquidus temperature, owing 
to the small size of the crystals. No ade- 
quate discussion of crystallization is pos- 
sible unless the compositions of the feld- 
spars are known. 

Upon examination of quenched sam- 
ples with the petrographic microscope, it 
was very easy to determine the tempera- 
ture of the first appearance of feldspar 
crystals for compositions with leucite as 
the primary phase, because the leucite 
crystals lie near the glass in index of 
refraction and the feldspar crystals with 
a much higher index stand out clearly. 
On the other hand, it was very difficult 
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to determine the temperature of the first 
appearance of leucite crystals for com- 
positions with a feldspar liquidus, and we 
were able to do this for only a few com- 
positions where the leucite appeared as 
large crystals which became birefringent 
on quenching after passing through the 
rapid inversion of pure leucite at about 
600° (Bowen and Schairer, 1929, n. 13, 
p. 311). 


THE JOIN SODA FELDSPAR~POTASH 
FELDSPAR~—LIME FELDSPAR 


Strictly speaking, soda feldspar-potash 
feldspar-lime feldspar cannot be rep- 
resented adequately on a_ triangular 
diagram. Because of the incongruent 
melting of potash feldspar and the ap- 


Fic. 2.—Diagram showing the position of the feldspar join (shaded area) in the tetrahedron NaAlSiO,- 
KAISiO,-SiO,-CaAl,Si,O3, and the compositions of the feldspars and ieldspathoids. 
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pearance of crystals of leucite (and, as a 
consequence, the appearance of liquids 
more siliceous than compositions which 
lie in the triangle), the results can be ex- 
pressed only quantitatively in terms of a 
four-component system. The behavior of 
the ternary feldspar compositions might 


R. R. FRANCO AND J. F. SCHAIRER 


INDICES OF REFRACTION OF GLASSES 


Homogeneous glasses were prepared 
for the forty-two ternary feldspar com- 
positions studied, and the index of refrac- 
tion of each glass was measured at 25° 
with the petrographic microscope by 
comparison with standardized immersion 


LIME 
FELDSPAR 


0 


FELDSPAR 


WEIGHT PERCENT 


10 


FELOSPAR 


Fic. 3.——Refractive indices of glasses at 25° C. 


be represented as a triangular join in a 
tetrahedron’ representing the system 
NaAlSiO,-KAISiO,-SiO,-lime feldspar.* 
The relation of the triangular join soda 
feldspar—-potash feldspar-lime feldspar 
to silica and to the feldspathoids nephe- 
line, leucite, and KAISiO, (kalsilite and 
other polymorphs) is shown as figure 2. 


3 For a discussion of the use of a tetrahedron to 
represent relations in a quaternary system see 
Schairer (1942, pp. 246-249). 


liquids. The values obtained are given in 
the first column of table 1. The values are 
accurate to 0.003. No attempts were 
made to obtain close accuracy, the 
glasses were not annealed, and the meas- 
urements were made in white light. The 


4 Strictly, this system is not completely quater- 
nary because of the appearance of crystals of Al,O, 
as a primary phase in some compositions in the sys- 
tem KAISiO,-lime feldspar (unpublished data of 
Schairer). 
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LIQUIDUS TEMPERATURES OF FELDSPARS 


data are plotted in figure 3. The compo- 
sitions measured are indicated by black 
dots, and isofracts or curves of equal 
index of refraction are plotted in the 
triangle. 


VALUE OF THE DATA 


The liquidus data on ternary feldspar 
compositions reported here is of value 
principally as a starting point for studies 
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of these same compositions at lower 
temperatures with water as an additional 
component. In the presence of water at 
moderate pressures, even at relatively 
low temperatures, equilibrium between 
crystals and liquid is attained rapidly, 
and such studies should answer many 
questions concerning the amount of solid 
solution in feldspars as a function of 
temperature. 
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DISCUSSION 


ON THE BIAS OF GRAIN-SIZE MEASUREMENTS 
MADE IN THIN-SECTION: A DISCUSSION® 


NORMAN N. GREENMAN 
University of Chicago 


INTRODUCTION 


The problem of obtaining the characteristics 
of the size distribution of sediments from thin- 
section data has been attacked by Krumbein 
(1935), among others, who derived the equation 


where the integral on the left is the nth moment 
of the thin-section distribution, the integral on 


> 


Fic. 1 


the right is the nth moment of the actual dis- 
tribution, and C is a constant depending upon 
n and derived on the assumption of spherical 
grain shapes. He applied this theory to two sand 
samples—the St. Peter sandstone and a re- 
worked glacial sand—and showed that there 
was good agreement between the calculated 
and the observed values of the first and second 
moments of each. 


* Manuscript received April 27, 1950. 


In a recent paper Chayes (1950) has exam- 
ined Krumbein’s work and concluded that the 
“agreement between observed and expected 
values is independent of grain shape, however 
defined or measured” (p. 158), that this agree- 
ment “is equally uninformative about whether 
they [the grains] vary in size”’ (p. 159), and that 
it “tells us next to nothing about size itself” (p. 
159). His final paragraph carries the implication 
that, since this agreement is independent of both 
shape and size variation of the grains, one is not 
justified in assuming on the basis of Krumbein’s 
experimental results that the use of his correc- 
tions will lead to unbiased estimates of the 
characteristics of the original size-frequency dis- 
tribution. It is the purpose of this paper to show 
that this implication is not justified by Chayes’s 
presentation of the case and that his statement 
concerning shape, although partially true, does 
not follow from his attempted proof. As a mat- 
ter of fact, shape affects the results in a way 
that neither Chayes nor Krumbein has brought 
out. 

SHAPE AND THE BIAS CORRECTIONS 


Chayes’s statement concerning shape is 
based only on an analogy. In figure 1, a pointer 
of constant length, L = 204, is rotating in the 
BCD-plane about the origin O, with a certain 
angular velocity. At random intervals of time an 
observer looking along BO notes the length of 
the projection of ZL upon CD, and a distribution 
of projection lengths is obtained. This amounts 
to replacing the thin-section, in which a popula- 
tion of grains is measured at a “single instant” 
in time with reference to a single plane, the sec- 
tioning plane, by a single “grain,’’ which is 
measured at many instants in time with refer- 
ence to all possible sectioning planes normal to 
OB. lf the angular velocity is such that the 
length of the projection of the pointer on OB 
(i.e., length 20A’) varies linearly with time? we 

?Chayes erroneously states that the angular 
velocity of the pointer is proportional to the sine of 
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have an analogy to Krumbein’s system because 
Krumbein assumed that the probabilities of 
cutting between equal intervals along OB are 
equal. Therefore, we can estimate the length of 
the pointer by applying Krumbein’s correction 
to the observed distribution of projection 
lengths. 

Now, says Chayes, although it is true that 
the pointer describes a circle as it rotates, we 
can also consider the pointer as remaining at 
rest and the observer rotating with the same 
angular velocity around O, measuring the pro- 
jection lengths upon a plane normal to his line 
of vision as before. The length L, he continues, 
can still be estimated, although the pointer 
describes no circle at all and the radius of the 
circle described by the observer is not given. He 
concludes, therefore, that shape is not a factor 
and that we can also measure the average c-axis 
of platy plagioclase crystals quenched in a lava 
flow or the a- and b-axes of rodlike hornblende 
crystals, in neither of which cases would grain 
shape affect the results. 

This conclusion by no means follows from the 
analogy, for two reasons. First, Chayes is say- 


the angle it makes with CD in order for its projection 
on OB to vary linearly with time. If we denote this 
this angle by @ and the length of the projection of the 
pointer on OB by y, then 
v=Lsin 0 

and 

dy dé 

—-=Lcos @—-=K 
dt 


where K is a constant. Since angular velocity, w, is 
defined as time rate of change of angle, we have 


= 

dt Leoos@ L 

so that the angular velocity is proportional to sec @, 
not sin @. 


sec 6, 


3 That is to say, the pointer projection length on 
OB must increase equal amounts in equal time 
intervals for the analogy of the rotating pointer to 
fulfil the requirements of Krumbein’s probability 
assumption. Thus if, while the pointer is in the 
first and third quadrants, its projection on OB 
varies from o to 20A’ in the time interval A/, then 
in the interval 2A¢ this projection must be 404’; 
in the interval 34¢, 604’; and so on until the pointer 
is parallel to OB. Because of symmetry, the same 
reasoning applies when the pointer is in the second 
and fourth quadrants, except that here the projec- 
tion length is diminished by 204’, 40A’, 60A’, etc., 
in the intervals Af, 2A/, 347, etc. 
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ing, in effect, that the distance of the observer 
is the analogue of grain shape. But it is clear 
that the analogue of grain shape is the path de- 
scribed by the end points of LZ, not that de- 
scribed by the observer. At any instant in time 
the length of projection of the pointer depends 
only upon ZL and the angle it makes with the 
projection plane; this length is entirely inde- 
pendent of the position of the observer along 
the line of sight when he makes the measure- 
ment. Moreover, it must also be clear that, with 
a given angular velocity, the distribution of 
projection lengths will, in general, be different 
for different rates of change in length Z as the 
pointer rotates; that is to say, since the path of 
the pointer ends is analogous to grain shape, 
each path corresponds to a certain shape, and, 
in general, each will give rise to a different dis- 
tribution of projection lengths. 

For example, suppose the angular velocity is 
proportional to sec @ and L is constant; we then 
have the analogue of spherical grains, and the 
distribution function will be given by 


as Krumbein has shown. Now, however, suppose 
the pointer so varies in length as it rotates that 
its ends describe a straight-line path. For ex- 
ample, if a plane is passed through any four 
coplanar vertices of an octahedron, the outline 
of the cross section so formed is a path of this 
sort, which might then be considered the ana- 
logue of octahedral shape. The equations of the 
four straight lines forming this path are: 


+ y = in the first quadrant , 
— x+y = }Luux, in the second quadrant , 
— y = ox, in the third quadrant , 
* — y = $2 nox, in the fourth quadrant , 


where Lmax is the maximum length of the 
pointer. Because of symmetry, only the first 
quadrant need be considered, so that in this 
case the distribution function becomes 


P(x) dx = — 


| 
In order to estimate the length L from the ob- 
served distribution of projection lengths, we 
must find the first moment, or mean, of the dis- 
tribution. In the circular case we haye 


LA 
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and in the linear case we have 


1 


Lmax 


— xdx = 


Thus, to estimate Z we must multiply the ob- 
served mean by 4/7 in the circular case and by 
2 in the linear case. Still other bias corrections 
would have to be applied if we assumed the 
ends of L to describe a sinusoidal path or a para- 
bolic path or any other type of path. Therefore, 
agreement between observed and expected 
results in a measurement can occur only where 
the correction factors are derived from equa- 
tions that adequately describe the shape of the 
grains being measured. Such agreement could 
not be obtained by applying Krumbein’s cor- 
rections to measurements of platy plagioclase or 
rodlike hornblende crystals because the shapes 
of these crystals are best described in terms of 
lines and planes, not circles. 

The artificiality of rotating the observer and 
the projection plane about the stationary 
pointer in order to have the pointer ends 
describe no path at all proves nothing and by 
no means nullifies the discussion just presented 
because the determining factor is the path 
described by the pointer ends with reference to 
the observer and the projection plane, irrespec- 
tive of which system is in motion. If we fasten 
the minute hand of a clock to a piece of glass 
and allow the works to rotate freely about the 
fixed hand, we will still be able to obtain the 
correct time, although the hand does not move 
at all; the face will merely rotate about the 
fixed hand, and the apparent motion of the lat- 
ter relative to the former will provide us with 
the necessary information. If we wish to meas- 
ure the rate of lengthening of the shadow of a 
tree, we will obtain the same results whether we 
consider the sun as moving about a fixed earth 
or the earth as rotating in the light beam cast by 
a fixed sun. It is entirely immaterial which of 
the two systems, “‘observer-projection plane” 
and “pointer,” is rotating and which is fixed or 
even whether both are rotating in either the 
same or opposite directions. The important fact 
is that the path described by the pointer 
ends relative to the observer and projection 
plane is precisely analogous to grain shape. 
It has already been demonstrated that varia- 
tions in this path cause variations in the dis- 
tribution of projection lengths. Therefore, by 
analogy, variations in grain shape will result in 
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varying size distributions measured in thin- 
section, so that agreement between expected 
and observed results of such measurements will 
depend upon the selection of proper correction 
factors, namely, those based on equations that 
adequately describe the shape of the grains. 

The second flaw in reasoning from this anal- 
ogy has already been hinted at, namely, that 
when we measure projections of the pointer 
upon a plane we are, in effect, measuring the 
diameters of circular sections of a sphere pro- 
duced by all possible sectioning planes parallel 
to the plane of projection, whereas in the case 
of a thin-section we are measuring circular sec- 
tions of many grains produced by only one of 
these sectioning planes. It is apparent, there- 
fore, that rotating pointer measurements corre- 
spond to thin-section measurements only in 
those cases in which the shapes of the grains are 
such that the dimension measured in thin-sec- 
tion has the same length as the projection of the 
pointer upon the plane of the thin-section (i.e., 
the plane normal to the OB-axis) and in no 
other case. This condition is fulfilled only by 
the sphere generally and by special cases of 
other shapes. 

To illustrate this point, let us consider, first, 
the case of spherical grains. Here the analogue 
is a pointer of constant length ZL. When the 
pointer makes an angle @ with the plane of pro- 
jection, its projection length is exactly equal to 
the diameter of the circular section of the sphere 
produced by a sectioning plane that is at a dis- 
tance L sin 6 from the sphere center. This 
situation holds for all angles, and the analogy, 
therefore, holds also. Now, however, consider 
a different shape—a cylinder of length / and 
diameter d. The projections of / and d will be 
equal to / cos 6 and d sin @, respectively, for all 
angles, so that rotating pointer measurements 
will yield a distribution of lengths varying from 
o to / if we measure / or from o to d if we meas- 
ure d. In the thin-section, on the other hand, we 
will usually observe elliptical sections with 
minor axis of length d and major axis of length 
d/sin @ for most angles. Thus in every elliptical 
section the short dimension will be true d, and 
we have to make no bias correction whatsoever. 
The distribution of long dimensions, however, 
will have to be corrected by some factor other 
than that used in the case of projection measure- 
ments in order to obtain true /.4 


4 The foregoing discussion neglects (1) the situa- 
tion at the ends of the cylinder, (2) the case for 
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The distributions in the two cases are thus 
markedly different, the difference arising from 
the fact that, in general, projection measure- 
ments are entirely unrelated to measurements 
made on a single section precisely because of the 
effects of shape. Therefore, because Chayes’s 
analogy holds generally only for one special 
shape—the sphere—and breaks down for all 
others, we cannot with confidence draw any 
conclusions from it as to the effect of grain 
shape on the agreement between expected and 
observed results of thin-section measurements. 

The manner in which grain shape affects the 
bias corrections has been brought out in Wick- 
sell’s work on this problem (1925, 1926). Wick- 
sell applied a rigorous mathematical treatment 
to both spherical and ellipsoidal grain popula- 
tions, taking into account not only the distances 
of the grain centers from the sectioning plane 
but the probability of a grain’s being cut by 
this plane as well. Krumbein’s solutions are 
thus a special case of Wicksell’s, where, in the 
latter’s notation, 


F(r) =f(r), 


where F(r) denotes the size-frequency distribu- 
tion of the actual population and /(r) the size- 


frequency distribution of those grains cut by 
the sectioning plane. It is shown that in this 
simplified case the same corrections apply, 


@ = go’, and (3) the case for @ = o and small angles 
in the vicinity of @ = o. Under the first head, there 
is the complication that cuts near the ends will pro- 
duce section boundaries that are part ellipse and 
part straight line. In the second case, all sections 
will be circles of diameter d, and thus no bias cor- 
rections are required for this particular cylinder 
orientation. In the third case, cuts through the 
ends will produce rectangular sections of length / 
for 6 = o and of length slightly greater than / for 
small angles near @ = o, the limiting values of such 
angles depending upon the relative dimensions of 
the cylinder. For the special case 6 = 0 no bias 
corrections are required to obtain true /, whereas d 
values—i.e., the short dimensions of the rectangular 
sections—must be corrected by Krumbein’s cor- 
rection factors to obtain true d. However, if / is 
sufficiently large compared to d, these less usual 
sections will be in the minority, the majority con- 
sisting of the elliptical sections discussed in the 
text. In any case, these complicating situations do 
not destroy the illustrative value of the cylinder 
because it was selected only to show that, in general, 
the effect of particle shape is to produce thin-sec- 
tion size distributions that are completely unre- 
lated to pointer-projection size distributions. 
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whether the grains are spheres or ellipsoids of 
any eccentricity and of any orientation. This 
result comes about because of a special relation 
between the sphere and the ellipsoid, namely, 
that the ratio of the area of any section of an 
ellipsoid to the area of a parallel central section 
is the same as the ratio of the area of a simi- 
larly situated section of a sphere to the area of 
its parallel central section. This is not too sur- 
prising, in view of the fact that mathematically 
the sphere is a special case of the ellipsoid; the 
circle is often considered an ellipse with zero 
eccentricity. 

In Wicksell’s general solution for the ellip- 
soidal case the departure from sphericity results 
in different correction factors; however, he has 
calculated these factors for the special cases of 
the prolate and oblate spheroids, and it turns 
out that, even with eccentricities as high as 
0.7, the values for the first five moments do not 
differ much from those for the spherical case. 
Thus, whether one applies Wicksell’s general 
solutions or Krumbein’s simplified ones, there 
is little error involved in assuming the grains 
to be spherical. 

It should be observed, however, that both 
Wicksell and Krumbein obtained their solutions 
by using the equation of the ellipse; Krumbein, 
in fact, used the special case of an ellipse with 
zero eccentricity, namely, a circle. Thus 
Chayes’s statement concerning the effect of 
shape is true, but only with regard to grains 
that can be assumed to have ellipsoidal out- 
lines. As has already been shown, if we assume 
grains with octahedral shape, the distribution 
of “diameters’’ measured in thin-section will 
differ from that of the spherical case and the 
bias corrections in this case will be different 
from those of either Wicksell or Krumbein. 
Now, since an orthorhombic dipyramid bears 
the same relation to an octahedron that an 
ellipsoid bears to a sphere with respect to the 
property mentioned above (i.e., the ratio of 
areas of similarly situated parallel sections), the 
correction factors derived for the octahedron 
will apply equally well to a tetragonal or ortho- 
rhombic dipyramid. The same reasoning can be 
applied to other three-dimensional shapes: thus 
the corrections for the sphere will apply to the 
ellipsoid; those for the octahedron, while differ- 
ent from those for the sphere, will apply to the 
orthorhombic dipyramid; and a still different 
set derived for cubes will work also for all rec- 
tangular parallelepipeds. An important fact 
concerning the effect of shape on the bias cor- 
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rections consequently emerges, namely, that 
the corrections are independent of shape defined 
on the basis of sphericity alone but do depend 
upon shape defined on the basis of the equations 
describing the grain surfaces. 

It may be pointed out that here we have a 
case in which the mathematical model (i.e., the 
derivation of the bias corrections) must take 
account not only of the experimental procedure 
but of the material upon which the experiment 
operates. As Chayes indicates, agreement be- 
tween the expected and the observed results of 
a measurement can be expected where the ex- 
perimental procedure closely corresponds to the 
mathematical model we construct. However, it 
is important to note that we must first examine 
the material in order to set up the proper mathe- 
matical model. Thus if we wish to measure 
Chayes’s platy plagioclase crystals or rodlike 
hornblende crystals, we must first examine the 
material and select the proper equations describ- 
ing the grain shapes in order to construct the 
proper mathematical model. Otherwise, agree- 
ment of results cannot be expected, for applica- 
tion of ellipsoidally derived corrections will 
result in errors that may be great or small, de- 
pending upon the amount of error involved in 
ascribing ellipsoidal outlines to shapes defined 
by planes and straight lines. 


SIZE VARIATION AND THE BIAS CORRECTIONS 


In discussing the effect of size variation on 
the bias corrections, Chayes has satisfactorily 
demonstrated that the vaiues of these correc- 
tions are independent of whether or not the 
grains vary in size. This is indeed a fortunate— 
in fact, necessary—situation because, if the 
values did depend upon the size variation of the 
original grain population, we would be faced 
with the seemingly insoluble problem of obtain- 
ing information by applying corrections the 
values of which themselves depend upon the 
unknown information we are using them to 
obtain. However, in his discussion and in the 
implications that follow he seems to have over- 
looked an important point, namely, the exist- 
ence of a distribution of thin-section “diame- 
ters’’—Krumbein’s Q(x) curve. 

There are, after all, three elements to be con- 
sidered in a measurement of this sort. The first 
is the size distribution of an actual grain popu- 
lation; this is the unknown information we are 
seeking. The second is an operation acting on 
this population so as to produce a bias, the cor- 
rection for which we derive mathematically on 
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the basis of certain assumptions. The third and 
equally important element is an observed size 
distribution of grain sections measured in thin- 
section. Now it is this third element, not the 
second, that tells us how the grains vary in 
size, for the Q(x) curve bears a definite relation 
to the original F(r) curve. 

If we consider any hypothetical distribution 
of grain sizes and assume that all the grains in 
any class interval are of a size equal to the mid- 
point of that class interval (as is done in calcu- 
lating the moments of a distribution), then for 
each of these classes there will be a distribution 
of thin-section sizes that has the form of Krum- 
bein’s P(x) curve. The Q(x) curve will then be 
the sum of all these P(x) curves, and it is easily 
verified, by drawing a number of distributions, 
that to each F(r) curve there corresponds 
a unique Q(x) curve. Inasmuch as correspond- 
ing moments of the two curves have a constant 
ratio, a change in the moments of one will result 
in a corresponding change in the moments of 
the other. Thus, although the value of this ratio 
is independent of size variation, the Q(x) curve 
is not, and we can therefore obtain a reasonable 
estimate of the desired moments from the ob- 
served moments. 

The values of the constant ratios involved 
depend upon the assumptions of randomness 
and upon the equations assumed to describe the 
grain shapes; that is, both the experimental 
procedure and the material dealt with must be 
considered in deriving the bias corrections. 
Krumbein’s experiments show that his assump- 
tions sufficiently approximate the physical 
facts that his corrections can be used to give 
reliable results in measurements of similar 
sediments. 


SIZE DEFINITION IN THIN-SECTION 
MEASUREMENT 


The definition of size is a problem separate 
from that of correcting the bias of a thin-section 
measurement. The bias corrections, it is true, 
are not wholly independent of the size definition. 
For example, Wicksell found it convenient, in 
developing his corrections for ellipsoidal grains, 
to define size as the geometric mean of the long, 
intermediate, and short diameters, so that his 
corrections can be applied only to distributions 
where the size is so measured. With Krumbein’s 
simplified theory, however, the corrections can 
be used whether we measure the long diame- 
ters, the short diameters, or their arithmetic or 
geometric means. The point is that, assuming 
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size to be defined on the basis of a certain pa- 
rameter, the sectioning process displaces the 
moments of the size distribution of this parame- 
ter by a fixed amount to give the moments of 
the thin-section distribution of this same pa- 
rameter. 

Now all we demand of a thin-section meas- 
urement is that, given a certain size definition, 
we have a reliable set of correction-factor values 
that will enable us to estimate from the ob- 
served distribution the distribution we would 
obtain from a sample of loose grains. In other 
words, if we measure long diameters (or short 
diameters or any mean of the diameters) in thin- 
section, we would like our correction factors to 
give us a reliable estimate of the distribution of 
long diameters (or short diameters or any mean 
of the diameters) we would measure in the sam- 
ple of loose grains. This is all we demand, and 
this is all we get. We can bridge the gap between 
the thin-section and the loose grains, but, once 
having obtained the loose-grain distribution, we 
still face the problem of what is size itself the 
same as if we had started with a sample of loose 
grains in the beginning. The problem is thus 
outside the realm of the bias-correction problem, 
and it is unfair to ask that the latter solve the 
former, since the definition of size is a problem 
that is inherent in any and every kind of size 
analysis we make. 

For example, if we make a size analysis of a 
glacial till, we use three different size definitions. 
The coarse fraction can be measured with cali- 
pers, so that “‘size’’ may be the lengths of the 
long diameters of the cobbles; the arithmetic or 
geometric means of the long, intermediate, and 
short diameters; or the diameters of spheres of 
the same volume or some other such measure. 
The pebble and sand fraction is usually analyzed 
by means of sieves, and here “size’’ is defined by 
the dimensions of the sieve openings. The silt 
and clay fractions are frequently analyzed by 
sedimentation techniques, where “‘size’’ is de- 
fined in terms of the settling velocities of quartz 
spheres. The distributions of these three frac- 
tions are thus not strictly comparable; but, if 
we have a number of till analyses, we can obtain 
important information by comparing the respec- 
tive coarse fractions, the respective sand frac- 
tions, and the respective fine fractions. 

Similarly, if we measure long-diameter dis- 
tributions in thin-section and correct them to 
obtain the long-diameter distributions we would 
measure in loose-grain samples, we can obtain 
valuable information by comparing these results 
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not only with other thin-section measurements 
but also with measurements of long-diameter 
distributions of loose-grain samples. Thin-sec- 
tion size analyses may suffer from the problem 
of size definition, but no more so than any other 
kind of size analysis. 


CONCLUSIONS 


In summary, then, we can list the following 
conclusions: 

1. Because the bias corrections are derived on 
the basis of certain assumed grain shapes, to 
obtain reliable results it is necessary to select 
the mathematical curves that most nearly ap- 
proximate the shapes involved. The values of 
these derived corrections thus depend upon the 
material being measured. In the case of sand- 
stones Krumbein’s corrections based on ellips- 
oidal grain outlines give reliable practical re- 
sults, as his experiments have shown. 

2. Chayes’s statement that Krumbein’s ob- 
served agreement of experimental and expected 
results is uninformative about the variation in 
size of the grains is true, but misses the point. 
The experiment does not attempt to tell whether 
the grains vary in size; the information concern- 
ing the size variation of the actual grain popu- 
lation is obtained from the distribution of thin- 
section diameters, not from the bias corrections. 
The experiment was an attempt to establish the 
hypothesis that, in sectioning the grains, the 
moments of the original distribution are dis- 
placed by a constant amount, no matter how 
the grains vary in size. The fact that the ob- 
served ratios of the corresponding moments of 
the thin-section and loose-grain distributions 
are in close agreement with the mathematically 
expected ratios indicates that this hypothesis is 
valid. If the hypothesis is valid, we then have 
reliable values for these ratios—which are the 
bias corrections—and, by applying them to ob- 
served thin-section distributions of similar sedi- 
ments, we can obtain reliable estimates of the 
distributions we would observe if we measured 
samples of loose grains. Thus we are not elimi- 
nating the bias by knowledge about the grains 
(except to the extent that we derive the theo- 
retical correction values by making certain as- 
sumptions as to the grain shapes), but we are 
obtaining knowledge about the grains—i.e., the 
moments of their size distribution—by eliminat- 
ing the bias. We eliminate the bias by obtaining 
reliable values for the constant ratios, and 
Krumbein’s experiments show that his values 


| 

| 

j 
} 


274 


are sufficiently reliable to give good practical 
results. 

3. The problem of how to define size itself 
crops up in connection with thin-section size 
analyses just as it does with every other type of 
size analysis, but the fact that it exists no more 
destroys the value of this method than it does 
that of any other. At the present stage of the 
science of sedimentology every technique of size 
analysis requires the setting-up of an arbitrary 
definition of size to fit the procedure used and 
the computations based on such procedure. The 
thin-section method is no more guilty in this 
respect than is any other method. Although it 
is certainly desirable to have a more uniform 
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and fundamental definition of size, until the 
day arrives when such a definition appears we 
can obtain much useful information for com- 
parative and other purposes from the thin-sec- 
tion method as it stands at present. 
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ON THE BIAS OF GRAIN-SIZE MEASUREMENTS 
MADE IN THIN SECTION: A REPLY! 


F. CHAYES 
Carnegie Institution of Washington 


Mr. Greenman seems to overlook the fact that 
my paper was intended solely as a discussion 
of certain aspects of Dr. Krumbein’s, and 
where he charges me with errors of interpreta- 
tion, it is usually because he has not taken the 
trouble to compare our work with sufficient 
care. Granting that lead shot do not depart 
sufficiently from sphericity to affect the result, 
Krumbein’s first experiment has significance 
only as a test of whether or not the experi- 
mental situation justifies his choice of an “equi- 
distance” rather than an “equi-angle” (or some 
other) definition of randomness. Although he 
nowhere intimates that choice is involved—i.e., 
that other solutions are possible—it seemed to 
me that, except for this possibility, he might 
have rested with either the mathematics or the 
experiment and had no need of both. 

Carried away by this discovery, I was guilty 
—although not quite so guilty as Mr. Greenman 
suggests—of understating the difficulty intro- 
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duced by shape differences. Krumbein did seem 
to regard the agreement between experiment 
and arithmetic as in some fashion a verification 
of the sphericity of the grains in his sandstone 
specimens, for he went so far as to discuss the 
results in terms of roundness and sphericity in- 
dices. But if, with the same assumptions, shapes 
other than spheres or circles give the same agree- 
ment, the agreement is not ipso facto indicative 
of the sphericity of the grains. It was in this 
sense that I suggested the two were “‘independ- 
ent,” and I agree with Mr. Greenman that the 
word is too strong. Krumbein’s Minnitaki sand 
grains do not look particularly round to me, yet 
evidently they are “sufficiently round.” It seems 
quite possible that many acicular crystals might 
also be “sufficiently linear.” At any rate, the 
core of my argument was not that one should 
measure acicular crystals in this fashion but 
rather that one should not regard agreement be- 
tween observed and expected values as neces- 
sarily related to roundness or sphericity. 

Mr. Greenman is correct about my erroneous 
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specification of angular velocity for the equal- 
distance criterion in figure 1, but his analogical 
discussion of the pointer analogy is misleading. 
There is hardly any need for a model of the 
sphere other than the sphere itself. If he chooses, 
he may regard either the path of the pointer or 
the “distance of the observer”—whatever this 
may mean—as analogous to grain shape. But 
what about the pointer itself? It is intended as 
a model of a crystal very much longer than it is 
broad or thick, and its length may be taken as 
either the maximum or the minimum dimension 
of such a crystal. Its path is intended as a model 
of rotational movement in metamorphism or, 
less adequately, of movement in turbulent flow. 

I should like to assure readers of the Journal 
that I did not overlook “the existence of a dis- 
tribution of thin-section ‘diameters,’”’ for it was 
the existence of several such that first drew my 
attention to the problem. The question was 
whether one should apply the correction factors 
before comparing results. Where there is some 
possibility of an independent check, as in un- 
consolidated sands or lead shot, there may be 
some excuse for multiplying observed values by 
constants. Where there is no possibility of such 
a check, as in the “grain-size” measurements 
made on most metamorphic and igneous rocks, 
this is just playing with numbers. The observed 
distribution of x contains just as much informa- 
tion as that of kx, and, ‘n dealing with it, one is 
much less likely to overestimate the importance 
of small differences. 

Finally, this is a good opportunity to correct 
two errors Mr. Greenman has overlooked. In my 
note (p. 157) it is asserted that from the meas- 
urement of random chords the proper estimate 
of d would be 2m. This was based on a guess that 
the correction for distances of centers from sur- 
face would be of the equal-distance type sug- 
gested by Krumbein, while that for displace- 
ments from the center in the plane of the section 
would be of the equal-angle type. From experi- 
mental work on briquetted lead shot it now ap- 
pears that this was a bad guess. Both corrections 
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are of the Krumbein type, and d = (16/x?) m 
gives a much better approximation than d = 
2m. 

In a discussion of the insensitivity of the 
correction factors to variation in grain size (p. 
159), it is alleged that the observed second mo- 
ment is to be multiplied by §, “whether d be 
regarded as fixed or as variable,” but it is obvi- 
ous that, if d is fixed, its second and higher 
moments will be zero. Now the second moment 
of measurements made in section will certainly 
be greater than zero, and it follows that from 
data obtained in this fashion we will never de- 
termine that the grains are all of the same size; 
in fact, it will always appear that they differ, 
though, of course, the variation may be greater 
in one batch than in another. Mr. Greenman 
overstates the case when he agrees with some- 
thing I did not say, viz., that agreement of ob- 
served and expected results is “uninformative 
about the variation in size of the grains.’’ The 
thing we cannot find out is whether there is 
variation in size. If it exists, we shall get some 
measure of it, but from measurements made in 
section it will always appear to exist. 

Mr. Greenman’s decision that this matter is 
inconsequential again reflects his failure to com- 
pare my work with Krumbein’s. I guessed that 
Krumbein “used a number of lead shot, all of 
the same radius,” because he suspected that 
variation in grain size would upset the check 
between theory and experiment. The guess may 
be a poor one, but, unless Mr. Greenman has 
some better explanation of this choice of mate- 
rial, it is he, not I, who missed the point. 

Mr. Greenman is certainly correct in con- 
cluding that “shape affects the results in a way 
neither Chayes nor Krumbein has brought out,” 
and if he too has failed to bring it out in sim- 
ple usable form, no one familiar with the prob- 
lem will wish to cast the first stone. His work 
shows promise, but it might have shown more 
than that if he had not spent so much time 
treating a highly specific argument as if it were 
a statement of general principles. 
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Le Modelé périglaciaire. By J. Tricart. Paris: 
Centre de Documentation Universitaire, 
1950. Pp. v+ 267; figs. 84. Planographed. 


Le Modelé périglaciaire is an outline pre- 
pared for students as part of a more extensive 
course in geomorphology. Two volumes dealing 
with structural geomorphology—(1) Le Relief 
des cétes and (2) Les Types de bordures de massifs 
anciens—have already appeared. This volume 
is the first of a proposed series dealing with 
climatic geomorphology. In its entirety this 
proposed course reveals courage and vision on 
the part of its authors, A. Cailleux and T. Tri- 
cart. It represents an approach to the study of 
geomorphology which, though apparent in the 
attitudes of many geomorphologists in both 
America and Europe, has heretofore not been 
presented at the student level. It breaks sharp- 
ly away from the trends of American textbooks 
and offers, instead, an approach which is more 
likely to render the subject applicable to re- 
lated fields of study. It implies that geomor- 
phology can no longer be studied at the broad 
sweep of the arm if, indeed, this were ever pos- 
sible. This system of geomorphology, as de- 
scribed by Tricart, consists of three parts: 

Part 
1. Azonal processes for the sculpture of relief 
Pluvial erosion 
Littoral erosion 
Eolian erosion 
2. Climatic geomorphology 
Topography of the cold lands 
Periglacial topography 
Glacial topography 
The topography of arid regions 
The topography of humid, forested lands 
3. Structural geomorphology 
The classification of rocks and the mechanism 
of their attack by erosion 
The structural types of relief 
Relief of hills 
Relief of depressions 
Relief of folds 
Relief of volcanoes 
The major structural units 
The geomorphology of the old massifs 
The sedimentary basins 
The recent mountain chains. 


It is still too early to see how this outline will 
be treated in volumes proposed but not yet 


written. Broadly defined, these units could en- 
compass the subject. Perhaps, however, a 
fourth part should be included, such as: 
4. Historical geomorphology 

The effect of tectonic changes 

The effect of climatic changes 

Unconformities. 

This approach, particularly as it is exempli- 
fied in Le Modelé périglaciaire, shows how de- 
tailed is the knowledge of climatology and how 
minute are the observations of soils required in 
modern studies of geomorphology. They are 
fully as important as structural geology and 
petrology. All these topics have become equally 
essential tools in current geomorphic thinking. 
Geomorphic studies described here show how 
very local climatic conditions may play an im- 
portant part in weathering and erosion. Micro- 
climatic studies, particularly when related to 
the climate of the soil, must therefore become 
increasingly important. 

The discussion of periglacial topography is 
presented here in a simple, clear, and direct 
style. The outline form which Tricart has used 
makes for easy reading and easy reference. It 
has enabled him to condense a remarkable 
amount of information into a relatively small 
volume. In some instances, however, there are 
statements or points of argument which merit 
amplification and further discussion. Peri- 
glacial topography is thoroughly considered in 
three chapters. In the first chapter the peri- 
glacial domain, the climatic environment of 
periglacial regions, and the behavior of ice in 
the soil are treated. The second chapter is de- 
voted to cryopedology, or periglacial soils. Here 
Tricart, in order to avoid commitment to a 
genetic hypothesis, has introduced a descriptive 
classification of periglacial soils. This classifica- 
tion is generally applicable and is therefore 
given below in translation: 


CLASSIFICATION OF PERIGLACIAL SOILS 


A. Soils with geometric form 
1. Soils with closed form 
a) Simple forms with arrangement or sorting 
(normal polygons and stone rings) 
6) Simple forms without sorting 
(tundra polygons) 
c) Complex types 
(subdivided polygons) 
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2. Soils with open form 
a) Striped soils 
6) Cryoturbation folds (i.c., solifluction ter- 
races) 
B. Amorphous soils 
1. Soils with protuberances 
a) Turf-covered mounds 
(Thufur, earth mounds, Erdhiigelbéden) 
6) Tundra hillocks 
(ground-ice mounds, pingos, palsen, hy- 
drolaccoliths) 
c) Twisted peat moss 
(Streifenstimpfe, Strangmoore) 
1a) Sloping soils with folds 
2a) Soils with steps 
(sheep tracks, terracettes) 
2b) Soils with swelling by cryoturba- 
tion 
Level soils without sorting 
2a) Mud fields 
2b) Stone fields 
2c) Stone pavement 
1c) Block streams and sorted rubble 
2a) Sorted rubble 
2b) Rubble flows 


In the course of the discussion a new word, 
“cryosol,”’ is introduced. This word is asynonym 
of “congeliturbate” and of “frost soil.” It ap- 
plies to any body of material modified by frost- 


action, both that which is orderly and that 
which has an amorphous arrangement. As a 
general term it has merit, but it is formed of an 
odd combination of a Greek and a Latin root. 
The third chapter deals with periglacial topog- 
raphy. Like the preceding, this chapter is com- 
plete in its content and logical in organization. 
Not only does it contain a discussion of such 
well-known topics as congeliturbation and 
rubble-mantled slopes, wind action and loess 
deposition, and sucn features of poorly drained 
conditions as the thaw-lakes and tundra hil- 
locks, but it also includes a discussion of much 
less thoroughly developed concepts, such as 
differential congelifraction and the asymmetry 
of valleys. The argument would, however, have 
been strengthened if, while discussing the forms 
produced by frost, more attention had been 
given to criteria for distinguishing similar non- 
periglacial forms produced under the combined 
action of chemical decomposition, gravity, and 
viscous soil flowage. At the end of each chapter 
there is an excellent annotated bibliography of 
important papers and books related to peri- 
glacial geology. In all, 262 references are cited. 

Two concepts included in this discussion are 
noteworthy for their novelty. Certainly they 
merit further discussion. One, attributed to a 


personal communication from Cailleux, is the 
concept that chemical weathering can, by form- 
ing an impervious residuum, render a surface 
immune against the action of frost. The second 
is the concept that peneplanation is impossible 
in at least some periglacial environments. Tri- 
cart assumes the dominant periglacial process 
of denudation to be equiplanation, as described 
by Cairnes from Alaska in 1912, whereby a 
gentle surface is formed, by the destruction of 
the slopes and filling of the valleys, at an alti- 
tude intermediate between that of the summits 
and the valley bottoms. It seems to the reviewer 
that this is a late mature or early old age stage 
of the periglacial cycle of erosion. It may repre- 
sent the maximum development of this cycle, 
where periglacial processes have operated only 
within the relatively short span of Pleistocene 
time. Congelifraction and congeliturbation do 
not, however, cease at this stage. Wave action 
and congelifraction produced by sea water 
must continue to occur along the coast. Nick- 
points will therefore continue to be formed at the 
shore line and must retreat headward along the 
flowing streams. Base level is thus as significant 
in periglacial regions as in temperate regions. 
A frost-produced peneplain, or surface of 
cryoplanation, must therefore remain a possi- 
bility, as was pointed out by Bryan in 1946. 
This excellent outline provides us, then, with 
the most complete summary available of peri- 
glacial geology and with an approach to the 
subject which is both logical and effective. It is 
valuable to the student of soils and surface 
geology for its content and to the teacher for its 
educational approach. 
Louts C. PELTIER 


West Boylston, Massachusetts 


Historical Geology. By Cart O. DuNBAR. New 
York: John Wiley & Sons, 1949. Pp. xii+ 
567; figs. 380. $5.00. 


Seventeen years ago, when reviewing the 
third edition (1933) of “‘Schuchert and Dun- 
bar,”’ Dr. Carey G. Croneis, president of Beloit 
College, made the following statement (Jour. 
Geology, vol. 42, p. 335): 

From the sidelines it has been interesting to 
watch the metamorphosis of that fine old Yale text- 
book team of Pirsson and Schuchert. On the physical 
side Longwell, Knopf, and Flint have now been in 
the line-up in place of Pirsson for some time; and 
recently a new man, Carl Dunbar, has been added 
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to the historical roster . . . it is safe to say that the 
new men can play the game quite as well as their 
elders. Of course, one does not have to read far. . . 
to see that the veteran Schuchert is still captain, 
yet nevertheless there is readily discernible a new 
attack. ... 


In a reply to this review by Dr. Croneis, a 
little poem was soon forthcoming from a 
“down-easterner,” and I think it interesting 
and worth repeating here: 

I’m happy to see, says hoary C. C. [Carey Croneis] 
The veteran sports reviewer 

That the new Yale bunch still carries the punch 
That makes their competitors bluer. 


I’ve watched Yale play since that far-off day 
When I put on my first long trouser, 
They’re a fine lot of men, I'll encourage them, then, 
I’ll give them a good write-up, yowser. 
They face somewhat east and their planets at least 
From the only true gospel have run far 
But I'll praise without stint, 
Knopf, Longwell, and Flint, 
Cap. Schuchert and even young Dunbar. 


These interesting little side lights on earlier 
reviews of the fine textbooks that have con- 
sistently come from the “Yale team”’ give a 
humorous, yet accurate, picture of the his- 
torical development of this series. Another 
chapter must be written on the development of 
the “Yale team,” for now Dunbar is the coach- 
captain-player of the historical team and has 
the sole responsibility for every play. If Dr. 
Croneis should review another “game,” I am 
sure he would say that the style of play is the 
same, but the attack is more polished and the 
running game much smoother. The “Yale sys- 
tem” is so thoroughly entrenched that even 
Dunbar continues to use “we.” 

This new Historical Geology is divided into 
seven principal parts, namely, (1) “‘Prologue’’; 
(2) “Before the Cambrian”’; (3) ““The Paleozoic 
World”; (4) Mesozoic World’’; (5) “The 
Modern World Unfolds”’; (5) “Epilogue’’; and 
(7) “An Introduction to Animals and Plants.” 
There are 567 pages in all, with some 450 excep- 
tionally fine pictures and illustrations, many 
completely new, making up about 33 per cent of 
the book. Pictures repeated from the fourth edi- 
tion were considerably improved by increasing 
the size and reproducing them on a better grade 
of paper. By placing most illustrations at the 
top or bottom and extending them to the edge 
of the page, the space allotted pictures was 
more than doubled at little or no expense to the 
space for text. Many of the pictures are original 
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photographs by the author; these are supple- 
mented by some excellent air photos, photo- 
graphs of outstanding paintings and murals of 
various museums, as well as many taken by 
professional photographers. The book is printed 
with the aid of large and clear type on superior 
paper, both of which make it a better textbook 
from the student’s point of view. The sentences 
are more concise, yet convey the same meaning 
as the longer, more involved sentences of the 
fourth edition. The phraseology is modern. 

Sixty-seven pages are given to the four chap- 
ters of the “Prologue.” The titles of these 
chapters are (1) “Records in Stone,” (2) “The 
Scale of Time,” (3) “Fossils, a Living Record 
of the Dead,” and (4) “The Constant Change 
of Living Things.”” The basic material varies 
but little from the fourth edition, but it has 
been rearranged to advantage. 

The second part, “Before the Cambrian,” 
is divided between (1) “Origin of the Earth” 
and (2) “The Cryptozoic Eon: Pre-Cambrian 
History of the Earth.”” The words of the author 
express the objectives: 

To avoid the monotony inherent in the system- 
atic account of period after period, emphasis is 
varied from chapter to chapter. In discussing the 
Cambrian, for example, considerable space is devot- 
ed to paleogeography and to the bases for subdivid- 
ing the rocks into series and formations; in the 
Ordovician chapter, the Taconian orogeny is dis- 
cussed at length because it illustrates the principles 
used in recognizing and dating all later orogenic dis- 
turbances. In this way, general principles of inter- 
pretation are developed one after another in the 
early chapters, leaving room in later ones for greater 
detail. This we believe to be fitting since human 
interest in rocks and fossils increases as we approach 
the modern world. 


To the reviewer it is difficult to reconcile the 
author’s statement, “general principles of inter- 
pretation are developed one after the other in 
the early chapters,”’ with the fact that only 35 
pages, or the same number apportioned the 
Cretaceous, is assigned to the entire pre-Cam- 
brian. Some geologists believe that there are 
enough principles involved in a study of the 
pre-Cambrian to allow this section more space, 
for chronologically it is here we must first 
study. It is a matter of question when one thinks 
“human interest in rocks and fossils increases as 
we approach the modern world.”” Maybe the 
teachers could create an interest in the pre- 
Cambrian if more emphasis were placed on this 
study and the principles involved. The author 
states that “about three-quarters of all geologic 
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time had transpired before the Cambrian 
period”’; yet only 6 per cent of the book—35 out 
of 567 pages—is devoted to a discussion of the 
history of 75 per cent of all geologic time. It ap- 
pears that we continue to sell the pre-Cambrian 
“short.” 

Part 3, ““The Paleozoic World,” is discussed 
in the next 179 pages. Dunbar has done an 
excellent job of reorganizing the chapters, 
simplifying sentences, and reducing the material 
in the paragraphs without loss of value and 
meaning. The paleogeographic panels are from 
Schuchert’s latest unpublished maps. The maps 
of the panels are nearly three times the size of 
those in the fourth edition, but they are ar- 
ranged in reverse order, that is, the map showing 
the paleogeography of the earliest or older part 
of the period is placed at the top of the page, 
and the youngest is at the bottom. These maps 
are framed with clouds and in the words of the 
author, “the clouds have no meteorologic sig- 
nificance; they are merely a device to hide 
critical areas for which evidence is lacking or 
inconclusive.” This type of panel is used 
throughout the Paleozoic and Mesozoic, but the 
simpler, minus-the-clouds type is used for the 
Cenozoic, and I believe, aside from the aesthetic 
aspect, that it is much better for the student. 

Correlation tables were omitted in the first 
printing, since the author thought them too 
technical for the student and not complete 
enough for the teacher. Footnote references to 
the correlation charts of the Committee on 
Stratigraphy of the National Research Council 
are found in appropriate places. At the request 
of professors the correlation charts are now 
included in every copy. This should widen the 
scope of the book. 

Most of the plates of fossils are the same as 
those used in the fourth edition, but, by repro- 
ducing them on the better paper, the minor 
characteristics of the fossils are more easily 
deciphered. The enlarged photographs add evi- 
dence to the material of the text. Many of the 
old pictures and drawings depicting Paleozoic 
life are replaced by enlarged photographs of a 
great mural (by Rudolph Zallinger) in the Yale 
Peabody Museum. 

The story of the Mesozoic world is told in 
the next 85 pages. The material is essentially 
that of the previous edition. A change was made 
in the section showing the contact of the Tri- 
assic and Jurassic. Attention is called to the 
finding of a small bipedal dinosaur in the 
Navajo sandstone of Arizona, and tracks of 
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dinosaurs have been discovered in the slabby 
sandstone between the Wingate and Navajo— 
this may make the Wingate sandstone late 
Triassic and not early Jurassic, as indicated in 
the fourth edition. New pictures, aerial views, 
cross sections, and block diagrams add much to 
this section of the book. 

The four chapters of part 5, “The Modern 
World Unfolds,” are (1) “Physical History of 
the Cenozoic,” (2) “Ice Sculptures the Final 
Scene,”’ (3) “Mammals Inherit the Earth,” and 
(4) “The Coming of Man.” The physical history 
of the Cenozoic is written as a unit rather than 
being broken into individual periods. To some 
extent the chapter on ice sculptures has been 
rewritten. The Keewatin and Labrador centers 
of ice accumulation have been combined into 
one now called the “Laurentide Ice Center.” 
Recent studies and implications of pollen 
analyses of interglacial deposits with reference 
to interglacial climates and history are dis- 
cussed. A part is assigned the importance of 
fossil mammals as a guide to the history and 
climate of the Pleistocene. A more complete 
calendar of Pleistocene time (by Flint) replaces 
the calendar of the previous edition. It is inter- 
esting to note that Dunbar omitted the contro- 
versial subject, “causes of glaciation’”’; this topic 
was given several pages in earlier editions. 

““Man’s family tree’’ is covered more com- 
pletely with the anticipation of the chapter to 
follow, the “Coming of Man.” This is generally 
a very interesting part to the beginning student 
in geology, and often it is hurriedly passed over 
by some authors. In the “Coming of Man” a 
more complete calendar of human history (by 
Flint) has replaced the old. Recent discoveries, 
including those in Mexico, are discussed. 

The Appendix, “An Introduction to Animals 
and Plants,” was prepared for the beginning 
student in geology without previous training in 
biology. The chapter begins with a simple classi- 
fication of animal and plant kingdoms, following 
which the major groups of «nimals and plants 
are reviewed. Most teachers will find this a very 
important part of the book, and the students 
lacking a background of bio! »gy should find it 
extremely helpful. 

In summary, the simple and lucia writing; 
the carefully selected, beautifully reproduced 
photographs and illustrations; the well-arranged 
and clearly set-off paragraphs; the simple but 
adequate diagrams; the treatment and organ- 
ization of the up-to-date material; the rather 
complete list of references and collateral read- 
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ings; and the general excellence should make 
this work by Dunbar an outstanding textbook 
for the teacher and student of historical geology. 


B. 'T. SANDEFUR 
Michigan State College 


Pleistocene Terraces of the Susquehanna River, 
Pennsylvania. By Louis C. PELTIER. (Penn- 
sylvania Geological Survey Bull., no. G 23, 4th 
ser.) 1949. Pp. 158; figs. 49; tables 46. 
Concealed by this modest title is the first 

comprehensive American study of terraces from 
the point of view of Pleistocene climatic fluctua- 
tions in both glacial and periglacial areas. Re- 
cent European concepts of periglacial weather- 
ing and mass movements are supported by field 
evidence and are integrated with terrace inter- 
pretations and Pleistocene chronology in gen- 
eral. From it emerges the view that more land- 
scape features than have been generally realized 
are inherited from Pleistocene conditions and 
that climatic cycles provide a natural basis for 
correlating Pleistocene events. 

The subject matter can be divided into six 
general parts: (1) introduction; (2) summary of 
glacial geology; (3) weathering and soils; (4) 
periglacial deposits; (5) description of terraces; 
and (6) Pleistocene climate of central Pennsyl- 
vania. 

Pleistocene deposits in the area include pre- 
Illinoian gravels, Illinoian till and gravels, 
Binghamton (Tazewell) till and gravels, Valley 
Heads (Cary) gravels, and Mankato gravels. 
The gravels represent both glacial outwash and 
periglacial deposits resulting from intensified 
frost-action and mass movements. Only the 
Mankato gravels are exclusively periglacial. Se- 
quences of (1) marginal and frontal kame ter- 
races, (2) valley-train terraces, and (3) perigla- 
cial terraces are recognized in going away from 
the glacial boundaries. Deposits of various ages 
are differentiated on the basis of degree of 
weathering and soil development, lithology, 
stratigraphic position, and physiographic rela- 
tions. Correlations with events outside the re- 
gion are based on the principle that the climatic 
cycles of tne Pleistocene were contemporaneous 
and of world-wide extent. 

LELAND HorBERG 


University of Chicago 
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Structural Geology of Canadian Ore Deposits: A 
Symposium. Arranged by a COMMITTEE OF 
THE GEOLOGY Division, CANADIAN INSTI- 
TUTE OF MINING AND METALLURGY. Mon- 
treal: Mercury Press, Ltd., 1948. Pp. 948. 


This volume is divided into four sections. 
Chapter i, “General Papers,” contains “The 
Cordilleran Region of Western Canada,” “The 
Canadian Precambrian Shield,” “Mineral Oc- 
currences in the Appalachian Region of 
Canada,” and “Relationships of Minor Struc- 
tures to Gold Deposition in Canada.” 

Chapter ii treats of the Cordilleran region 
and has twenty-eight papers on the mines and 
mining districts of the coastal belt, interior belt, 
and Rocky Mountain belt. 

Chapter iii contain ninety-four papers on 
the ore deposits of the Canadian Pre-Cambrian 
Shield. It is divided into ten parts as follows: 
“Yellowknife-Great Bear Lake Area,” “Mani- 
toba and Saskatchewan,” “Red Lake—Patricia 
Area,” “Little Long Lac-Sturgeon River 
Area,” “Steep Rock—Michipicoten-Goudreau 
Area,” “Porcupine-Ramore Area,” “Sudbury 
Area,” “Matachewan-Kirkland Lake—Larder 
Lake Area,” ‘“Rouyn-Chiboughamau Area 
(Western Part),” and “Rouyn-Chiboughamau 
Area (Eastern Part).” 

In chapter iv there are five papers covering 
deposits in the eastern townships and Gaspe, 
and in Nova Scotia. 

Structural Geology of Canadian Ore Deposits 
is clearly and concisely written. It is illustrated 
with numerous geologic maps and sections and 
with photographs, all of which are well executed 
and reproduced. For a volume which empha- 
sizes structural features, however, there are 
many illustrations which lack pertinent struc- 
tural data. Many of the maps show no strike, 
dip, or lineation symbols which might indicate 
why an ore body occupies its present position. 
That changes in structural attitude are com- 
monly loci of ore bodies is known to every stu- 
dent of mineral deposits. Yet in many cases the 
structural data which would illustrate this are 
not given. This may be in part the result of 
company policy with regard to publication, but 
it is probably not true of all such omissions. 

The volume is an excellent reference work, 
covering the important operating mines and 
mining districts of Canada. Although structural 
relationships are emphasized, the treatment in 
many of the articles is much broader and in- 
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cludes sections on general geology, mineraliza- 
tion, alteration, and genesis. The papers were 
written by 121 contributors especially for this 
symposium and have been published to com- 
memorate the fiftieth anniversary of the found- 
ing of the Canadian Institute of Mining and 
Metallurgy. For the most part government 
geologists contributed the general papers, and 
mine or consulting geologists wrote the papers 
on individual mines. Some of the articles de- 
scribe relatively new areas, and others present 
up-to-date information on older mines. The 
volume contains data previously widely scat- 
tered through the literature. In some cases the 
information is new. For example, to this re- 
viewer’s knowledge, here is the first detailed 
presentation of the modern views of the Inter- 
national Nickel Company’s geologists on the 
origin of the Sudbury ore body. The information 
on the structure of the Sudbury irruptive and on 
the genesis of the ore is most welcome. It illus- 
trates once again the advantgaes of having de- 
tailed information, especially when this is com- 
bined with the ability to interpret and general- 
ize, before theorizing. The data on the Sullivan 
mine and on the West Bay fault are most wel- 
come, to mention only two of the areas which 
previously had been treated inadequately in the 
literature. 

In general, each district or region is treated 
as a unit and then the individual mines are de- 
scribed. It is to be regretted that some of the 
older nonoperating mines and mining districts 
have been omitted. An example is the Cobalt 
district, a discussion of which would have added 
to the understanding of the mineralization 
within the petrologic province which extends 
from the Duluth gabbro eastward into Quebec 
and includes the Michigan copper district, Sud- 
bury, Cobalt, Silver Islet, and a number of other 
silver, copper, and pyrrhotite deposits. 

It would have been helpful if the deposits had 
been classified according to the principal struc- 
tural control, as was done by Newhouse in Ore 
Deposits as Related to Structural Features. This 
was done for certain of the deposits in the first 
section of the volume, but a more comprehen- 
sive discussion and tabulation would have em- 
phasized the structural aspects with which the 
authors of the book were concerned. 

As stated in the Preface, the book “is not 
presented as a completed undertaking but 
rather as a ‘stock taking’ of existing knowl- 
edge.””’ The committee of the institute has 
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achieved this purpose admirably. A vast 
amount of information has been assembled in 
one volume. The book should furnish many 
facts necessary for further work on Canadian 
ore deposits, such as: correlations of the struc- 
tural relations of the various mining districts 
with the tectonic pattern of Canada, types of 
mineralization in the different metallogenetic 
epochs and provinces of Canada, types of ores 
related to wall-rock alteration, and the relation 
of structure to ore genesis. Much of the infor- 
mation required for such studies is not available, 
but a start has been made, and many references 
are given in this volume. Such studies should 
produce generalizations of considerable value to 
the science of mineral deposits. 

The committee of the institute is to be con- 
gratulated for this comprehensive compilation 
of information. The officers of the Canadian 
mining companies are also commended for re- 
leasing hitherto unpublished data on individual 
mines. 

A. F. HAGNER 
University of Illinois 


Geology: An Introduction to Earth History. By 
H. H. Reap. (“Home University Library of 
Modern Knowledge,” no. 198.) London, New 
York, and Toronto: Oxford University 
Press, 1949. Pp. 248; figs. 30. $2.00. 

On the dust jacket of this book is the state- 
ment that, although all the authors in the 
“Home University Library” series are experts in 
their subjects, “they are not selected for their 
academic distinction alone; they are also skilled 
in the art of presenting their material and able 
to unite learning with lucidity.” If this volume 
is typical of the series it is no idle boast. Almost 
every geologist, at one time or another, has been 
asked by intelligent friends to suggest to them 
a book on geology. Whether their aim is to be- 
come familiar with the history of the earth, to 
acquire a knowledge of the processes that shape 
the earth, or to obtain a picture of the aims, 
methods, and scope of the science, this volume 
appears to fill the bill. It does not smack of the 
average elementary textbook, with its multitude 
of facts and choppy, often pedantic style, which 
one hesitates to recommend as leisure-time read- 
ing to an intelligent person. Although many of 
the fundamental problems of earth history are 
dealt with, the book demands no special pre- 
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liminary acquaintance with geology. Yet Dr. 
Read skilfully avoids the pitfall of making ev- 
erything too pat. In urbane style he leads the 
reader gently, yet alertly, along, provoking and 
anticipating questions and adroitly introducing 
topics when they enter the reader’s mind as the 
logical consequence of what went before. 

The first of the five chapters, ‘““The Search 
for the Method,” gives the history of the succes- 
sion of ideas which founded the science of geol- 
. ogy and at the same time prepares the reader for 
the rest of the book It discusses the significance 
of fossils and of fossil-bearing beds, the growth 
of the ideas that rocks could be derived from 
pre-existing rocks and of correlation by fossils, 
and the realization that the principle of uni- 
formitarianism is the cornerstone of geology. In 
other words, it presents the entry of the con- 
cepts of space and time into the study of rocks. 

The second chapter, “The Sedimentary 
Rocks as Historical Documents,’ treats two 
main problems: (1) the sedimentary rocks as 
revealing the environments in which they were 
laid down and (2) the succession of these rocks 
as revealing a sequence of environments. Under 
the first heading is a flow-sheet and adequate, 
though brief, discussion of the various sorts of 
sediments and their modes of deposition. Under 
the second are such topics as the great inunda- 
tions, marine transgressions and regressions and 
their records, the question of the permanence of 
the ocean basins, a first statement of the divi- 
sions of the record of earth history, and a discus- 
sion of the methods used in estimating geologic 
time in years. 

The chapter entitled “Interventions from the 
Depths” begins with a discussion of volcanic 
rocks and proceeds to an exceilent presentation 
of the problems of regional metamorphism and 
granitization. Read divides the rocks into sedi- 
mentary, volcanic, and plutonic, the latter in- 
cluding granitic, migmatitic, and metamorphic 
rocks. The chapter continues with a discussion 
of the knowledge of the earth’s interior from 
earthquake waves and of the origin of the 
earth, in which there is an unfortunate omission 
of recent theories for which the astronomical 
arguments are strong. 

In chapter iv, “The Revolutionary Epi- 
sodes,”’ Read takes up what he considers “the 
short violent periods when the sediments were 
disturbed and land-masses extended and ele- 
vated,” a concept that will not be palatable to 
many geologists. However, his discussion of the 
way in which these “earth revolutions” are re- 
corded and how the record is read, and his pres- 
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entation of joints, faults, folds, geosynclinal his- 
tory, and isostasy is very good and is followed 
by an excellent discussion of the ultimate mech- 
anism of diastrophism, in which the contrac- 
tion, convection-current, and continental-drift 
hypotheses are examined pro and con, with no 
dogmatic final choice made. 

Chapter v, “The Pattern of Earth History,” 
weaves into a unified story the diverse strands 
which make up earth history. Here is a chapter 
that might have been merely a boring synopsis 
and repetition in “the seas came in, the seas 
went out” manner or a “Sunday supplement” 
recital of the wonders and oddities of the past. 
Instead, he presents a lucid and always interest- 
ing account of the true highlights of earth his- 
tory, selecting only sufficient episodes to pre- 
serve the continuity. Of these he examines a few 
in some detail as illustrative of profound and sig- 
nificant matters in earth history, among them 
the special methods developed for reading the 
pre-Cambrian record, the pattern of lower Pale- 
ozoic life, and the Coal Measures rhythm. It 
might be expected that an English book would 
be unsuitable for the American reader because 
most of the examples would be taken from the 
geology of Great Britain or at least from Eu- 
rope. However, the whole book and this final 
chapter in particular are not entirely the geol- 
ogy of Great Britain; they are world-wide in 
scope, with only a very slight, and wholly jus- 
tifiable, emphasis on Read’s own country. For 
instance, there is a long discussion of Finnish 
and Canadian pre-Cambrian records, of the his- 
tory of the Appalachian geosyncline, of the 
Gondwana series of India and Africa, and of 
Jurassic and Cretaceous events in North Amer- 
ica. 

Lov WILLIAMS PAGE 
University of Chicago 


Principles of Sedimentation. By W. H. TwEn- 
HOFEL. 2d ed. New York: McGraw-Hill Book 
Co., 1950. Pp. xii+673; figs. 81; tables 32. 
$6.50. 


The second edition of this well-known text 
and reference work follows the first after a lapse 
of eleven years. It is larger by fifty-six pages, 
and there are nearly twice as many illustrations, 
including diagrams, photographs, and micro- 
photographs. Among features newly illustrated 
are progressive and regressive marine overlap, 
soil erosion, microphotographs of sand grains, 
thin sections of coal, mud cracks, armored mud 
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ball, cone-in-cone, and stylolites. The pictures 
are well chosen and fairly well reproduced but 
still leave unillustrated many features which 
would make the book more understandable to 
the student. The bibliographic references are 
increased about 30 per cent with citations of 
more recent literature. The new illustrations, 
tables, and expanded bibliography occupy a 
little more than half the additional pages. The 
organization of the older book has been exactly 
followed, the principal parts of the book being 
(1) environments (2 chaps.); (2) origin of in- 
organic sediments (1 chap.); (3) interrelations of 
organisms and sediments (1 chap.); (4) trans- 
portation and deposition, including diagenesis 
(1 chap.); (5) classification (1 chap.); (6) de- 
scription of types of sediments (6 chaps.); (7) 
structures, textures, and colors (2 chaps.). Only 
one new minor heading on selenium as a minor 
constituent of sediments was observed. 

As the older edition of this book is familiar to 
most geologists, this review is largely concerned 
with the alterations in the new book. There are 
virtually no omissions of text from the first edi- 
tion, but in numerous places sentences have 
been rewritten to achieve greater clarity of 
statement. The work is well proofread, and very 
few typographical errors were observed. A care- 
ful comparative reading of the two books re- 
sulted in finding over a hundred separate addi- 
tional statements largely based on recent publi- 
cations and on the author’s personal researches, 
especially in the field of lake sediments. These 
are generally as brief as one to three sentences, 
and the largest single addition, a description of 
the lower Mississippi delta based on work by 
Russell and Russell, is two and a half pages long. 
Most of the new material is so well integrated 
into the text that it might be missed by the 
reader unless the two editions were being com- 
pared. 

The author throughout displays the caution 
in the acceptance of new interpretations for 
which he is well known. He expresses open dis- 
approval of some proposals, such as Russell’s 
theory of loess formation, Shaub’s theories for 
cone-in-cone and stylolites, Pettijohn’s proposal 
to return to Grabau’s endogenetic and exo- 
genetic classes of rocks, and his proposal to 
amend the definition of graywacke so that it is 
not the basic equivalent of arkose. He with- 
draws support earlier given to Galliher’s ex- 
planation of glauconite as an alteration product 
of biotite. He states (p. 549) that, although sus- 
pended loads of sediment may tend to settle in 
stratified arrangement, it is doubtful that this 
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would occur in natural waters. Recently (1950) 
Kuenen has shown that sediment transported 
by turbidity currents may settle naturally in 
stratified fashon. 

The interval since Twenhofel’s first edition of 
Principles of Sedimentation has been a period of 
prolific output of significant new literature in 
sedimentation and related fields. There has been 
a more searching analysis of the operation of 
transportational agencies, a shift in emphasis 
from simple environments of sedimentation to 
the tectonic framework responsible for the pat- 
tern of environments. Studies in paleoecology 
have expanded rapidly, and increasing attention 
is being devoted to postdepositional changes in 
sediments. True petrographic methods are being 
increasingly applied to the products of sedimen- 
tation, and these are being treated more and 
more by quantitative methods of analysis. Geo- 
synclines, facies, reefs, turbidity currents, peri- 
glacial phenomena, and phreatic versus vados 
origin of caverns have all become live subjects 
to which symposia, memoirs, and special num- 
bers of geological journals have been devoted. It 
seems unwarranted under these conditions to 
adhere so strictly to the old organization and 
even more so to omit from consideration numer- 
ous of these newer subjects. This manuscript 
was completed for the printer early in 1948, and 
since that date there have appeared such works 
as (1) Pettijohn’s Sedimentary Rocks; (2) 
Shrock’s Sequence of Layered Rocks; (3) Shep- 
ard’s Submarine Geology; (4) Trask’s Applied 
Sedimentation; (5) The Geological Society of 
America’s memoir on facies; (6) special issues of 
the Journal of Geology and the Bulletin of the 
American Association of Petroleum Geologists on 
reefs; (7) several papers by Krumbein, Sloss, 
and Dapples on tectonic controls on sedimenta- 
tion; and (8) papers by Kay on types of geosyn- 
clines and cratonic platforms. These and many 
other recent contributions point to the need for 
future revision of this book to bring it up to 
date. 

The reviewer searched in vain for references 
to, or excerpts from, numerous papers which he 
has considered outstanding contributions to 
sedimentation of the period from 1939 to 1947 
which were presumably available to the author. 

Thus in treating environments it is not indi- 
cated that glaciers may dissipate through stag- 
nation as well as by retreat (Flint); there is no 
citation or reference to the excellent treatment 
of valley flats by Happ and Rittenhouse in 
“Some Principles of Accelerated Stream and 
Valley Sedimentation”; and, though Fisk’s 
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name is used as a source of information on val- 
ley-flat sediments, there is no change in the text 
on this subject from the previous edition prior 
to the issuance of his monumental paper, Geo- 
logical Investigations of the Alluvial Valley of the 
Lower Mississippi River. Under the cave en- 
vironment no mention is made of Bretz’s re- 
searches into limestone caverns. The periglacial 
environment and permafrost escape mention. 
Under “Lagoons” no mention is made of 
Lucke’s work on Barnegat Bay and Inlet. 

In treating transportation and deposition of 
sediments the outstanding book of Bagnold, 
The Physics of Blown Sand and Desert Dunes, is 
mentioned only in connection with maximum 
height of dunes, though this includes a splendid 
treatment of saltation in air currents, its initia- 
tion and maintenance. Saltation is mentioned 
only in connection with currents of water. In 
glacial transportation Demorest’s contributions 
to glacial dynamics also escape mention. 

In the classification and description of the 
various types of sediments, the traditional 
method of size grades for the clastics and chemi- 
cal composition for the nonclastics is followed, 
and there is no mention of tectonic facies, con- 
trasting sediments of actively subsiding basins 
with those of stable platforms or cratonic areas. 
This has been the most significant new trend in 
the treatment of sedimentary rocks and is the 
underlying principle in much of Pettijohn’s 
Sedimentary Rocks. Undue space (p. 505) seems 
to be given to an interesting occurrence of 
doubly terminated quartz crystals in Permian 
gypsum for which the author offers no interpre- 
tation. 

In summary the reviewer regards this book 
as an excellent useful reference work in sedi- 
mentation which concisely digests the content 
of several hundred research papers generally, 
with appraisal as to their relative importance, 
but he believes that the work suffers because the 
organization and content of the earlier edition 
have been adnered to closely, without adequate 
modification in the light of the multitude of im- 
portant researches and new viewpoints of the 
past decade. 

HAROLD R. WANLEsS 
University of Illinois 
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Introduction to Theoretical Igneous Petrology. By 
ERNEstT E. WAnLstTROM. New York: John 
Wiley & Sons, Inc., 1950. Pp. ix+365; figs. 
150. $6.00. 

This book starts with an introductory chap- 
ter and a long description (38 pp.) of the ele- 
mentary concepts of heterogeneous equilibria. 
In my opinion this section is unnecessarily long, 
but it is in harmony with the spoon-feeding 
tradition of American teaching. 

The section, “Equilibrium in Silicate Sys- 
tems,” is useful but is séparated from that on 
“Igneous Minerals” for no good reason. Other- 
wise the chapters are well compiled. However, 
in the discussion of the feldspar system—the 
most important of all ternary rock-making sys- 
tems—the consequences of the determination 
by Schairer and Bowen (1947) of the equi- 
librium relationships between potash feldspar 
and anorthite should have been emphasized. 
The discussion, as it stands, is misleading, and 
so, in consequence, is that of the potash-rich 
rocks (pp. 306 ff.). 

Part 2 contains illuminating, but somewhat 
randomly arranged, chapters on the crust and 
the interior of the earth, on magmas, igneous 
differentiation, assimilation and syntexis, gran- 
itization, and origin of rock types, and closes 
with a brief survey of the classification of rocks. 
A short chapter on elementary physicochemical 
concepts is appended. 

There is great need for a text covering the 
formation processes of all types of rock. Al- 
though Wahlstrom’s textbook, as indicated by 
the title, covers only igneous rocks, it is the first 
of its kind in England. It is my hope that it will 
prepare the ground and induce American col- 
leges and universities to teach less of the purely 
descriptive petrography and present, instead, 
what to me seems the more useful alternative in 
the education of both the academic and the 
professional geologist, namely, an examination 
and interpretation of the geological processes re- 
sponsible for the formation of igneous rocks 
(and eventually of all rocks) in terms of 
physicochemical principles. 

T. F. W. Barta 


Universitetet i Oslo 
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vania Dept. Internal Affairs, Top. and Geol. 
Survey Bull. G 24. Harrisburg, 1949. 

Index to Well Samples. By Dan E. Feray and 
Jasper L. Starnes. Univ. Texas Pub. sors. 
Austin, 1950. 

McDonald and Adjacent Oil Fields, Allegheny 


and Washington Counties, Pennsylvania. By 
A. I. Ingham, E. M. Tignor, and W. M. Na- 
bors. Pennsylvania Dept. Internal Affairs, 
Top. and Geol. Survey Bull. M 29. Harris- 
burg, 1940. 

The Mining Groups of the Yilgarn Goldfield 
North of the Great Eastern Railway (Maps, 
Plans, and Sections). Western Australia 
Geol. Survey Bull. ror. Perth: William H. 
Wyatt, 1950. 

La Partie orientale du bassin de Paris: Etude 
morphologique. Vol. 1: La Genése du bassin. 
By Jean L. F. Tricart. Paris: Société d’Edi- 
tion d’Enseignement Supérieur, 1949. 

Report of the Geological Survey for the Year 
1947. Mines Dept., Western Australia, Geol. 
Survey Branch. Perth: William H. Wyatt, 
1950. 

The Sedimentary Rocks of Sokoto Province. By 
Brynmor Jones. Geol. Survey Nigeria Bull. 
18. Kaduna: Government Printer, 1948. 

Sierra Blanca Limestone in Santa Barbara 
County, California. By George W. Walker. 
California Div. Mines Special Rept. 1-A. 
San Francisco, 1950. 

Submarine Topography and Sedimentation in 
the Vicinity of Mugu Submarine Canyon, 
California. By L. D. Inman. Beach Erosion 
Board Tech. Mem. 19. Scripps Inst. Oceanog- 
raphy Contr. 473, new ser. 1950. 

Subsurface Reconnaissance of Glacial Deposits 
in Northeastern Kansas. By John C. Frye 
and K. L. Walters. Univ. Kansas Pub., State 
Geol. Survey Bull. 86, pt. 6. Lawrence, 1950. 

Summarized Record of Deep Wells in Pennsyl- 
vania. By Charles R. Fettke. Pennsylvania 
Dept. Internal Affairs, Top. and Geol. Sur- 
vey Bull. M 31. Harrisburg, 1950. 

Utilization of Texas Serpentine. By Virgil E. 
Barnes, D. A. Shock, and W. A. Cunning- 
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The Water Resources of Greene County, Ohio. 
By Stanley E. Norris, William P. Cross,and 
Richard P. Goldthwait. Ohio Dept. Nat. Re- 
sources, Div. Water Bull. 19. Columbus, 
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The West Indies and the Mountain Uplift 
Problem. 3 pts. By C. T. Trechmann. Pri- 
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COMMUNICATIONS AND ANNOUNCEMENTS 


AMERICAN ASSOCIATION FOR THE ADVANCEMENT OF SCIENCE, PHILA- 
DELPHIA, PENNSYLVANIA, DECEMBER 26-31, 1951, GENERAL 
GEOLOGY PROGRAM, SECTION E, TITLES OF PAPERS 


Participants are urged to file titles and ab- 
stracts by June 15, 1951, in order to be assured 
of a place on the program. The final date is Sep- 
tember 15. Abstracts are limited to 250 words 
and should be submitted in triplicate. No ab- 
stract will be printed unless the paper is pre- 
sented orally, and no author will be allowed 
more than 15 minutes for presentation except by 
special request. 

In addition to the general sessions in geology 
and geography, special programs on “The Crys- 
talline Rocks of the Older Appalachians” and 


“Foreign Petroleum Areas” and possible joint 
meetings with other sections are being planned. 
The vice-presidential address will be given by 
Kenneth K. Landes. 


Your co-operation in calling this meeting to 
the attention of your colleagues or other pos- 
sible participants will be appreciated. Participa- 
tion is not limited to members. 

Titles and abstracts should be sent to: Leland 


Horberg, Secretary, Department of Geology, 
University of Chicago, Chicago 37, Illinois. 


ERRATUM 


In “Niagaran Reefs of the Great Lakes Area,” by Heinz A. Lowenstam, vol. 58, no. 4, p. 443, 
l. 10 from bottom of first column, read The cavity fillings are therefore recent and diagenetic 


instead of The cavities and filling. . . . 
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AMERICAN PAULIN SYSTEM 


The World’s Finest... 


pe 7 MICRO SURVEYING ALTIMETER 


MODEL M-1 . . . Range 
6,000 feet (—1000’ to+ 
5000’) in intervals of 1’ 


MODEL M-2 . . . Ronge 
10,000 feet (0’ to 10,- 
000’) in intervals of 2’ 


MODEL M-5 . . . Range 


15,000 feet (0’ to 15,- 
000’) in intervals of 5’ 


MODEL MM-1 . . . Range 
5,000 meters (0 to 5,000 
meters) in intervals of 1 
meter 


ay $300 EACH with leother 
cose, Thermometer, Mag- 

nifier, and Operational 
Procedures. 


DEPENDABLE Only American Paulin 
LIGHTWEIGHT System Altimeters are 
graduated in easily 


ad 1 foot divisions. 


The World’s Standard... 


TERRA SURVEYING ALTIMETER 


MODEL SA-1 . MODEL SA-2 . . . Range | MODEL SA-5S . . . Range | $200 EACH with leather 

4,360 feet to+ 10,600 feet (—900" to-+ 15,000 feet (—500" to+ cose, Thermometer, Mog- 

3,600’) in intervals of 2’ | 9,700’) in intervals of 5’ | 14,500’)in intervals of 10’| nifier, and Operational 
Procedures. 


Literature and Technical Publications Available on Request 


AMERICAN PAULIN SYSTEM 


1847 S$. FLOWER © LOS ANGELES 15 


If You Need 


MINERALS OR ROCKS 


FOR YOURSELF OR SCHOOL 
Don't Pass This Up 


We have individual specimens 
sonably priced, and complete Mineral 
collections from 50 boxed ores at $4.00 
to museum collections of 550 3 x 4” 
superb specimens for $2,000.00. If you 
are looking for any mineral to complete 
+. = or a series consult us. We have 

scarce items, books on GEOL- 
OGY, MINERAL CHARTS, ULTRA- 
VIOLET LAMPS, GEOPHYSICAL 
EQUIPMENT, GEIGER COUNT- 
ERS, MICROSCOPES, MAGNI.- 
FIERS, GONIOMETERS, THEODO- 
LITES, etc. Send for large free Cata- 


logue. 


ECKERT MINERAL RESEARCH 
(Dept. G) 
110 E. MAIN ST. FLORENCE, COLO. 


GEOLOGY APPLIED 
TO SELENOLOGY > 
By J. E. SPURR 


“This is the first time that lunar features 
have been studied carefully by one trained 
and experienced in modern structural and 
igneous geology.’’— Journal of Geology. 
Vols. I and II combined, Fea- 


TURES OF THE Moon, 1945, 430 pp., 
95 text figures : 


Vol. III, Lunar 
ree 1948, 253 pp., 47 text 
Vol. IV., THe SuHrunken Moon, 


1949, 207 pp., 36 text figures 
$4.00 


Complete Set... ....... $12.00 


ROBERT A. SPURR 
College Park, Maryland 


THE THEORY OF 
GROUND-WATER MOTION 


By DR. M. KING HUBBERT 
Past Editor of Geophysics 


ASSOCIATE DIRECTOR, EXPLORATION AND PRODUCTION RESEARCH 
SHELL OIL COMPANY 


160 


pages 48 figures $2.00 
Postage: free in U.S.; Canada, $0.04; foreign, $0.10 


A fundamental treatise on a subject of great interest to hydrologists, petro- 
leum exploitation engineers, geologists, and all others concerned with the 


motion of fluids through porous media. 


Reprinted from Journal of Geology, vol. 48, no. 8, pt. 1 


THE UNIVERSITY OF CHICAGO PRESS 


5750 ELLIS AVENUE 


CHICAGO 37, ILLINOIS 
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First English Translation 


Silicate Melt 
Kquilibria 


By WILHELM EITEL 
Translated by J. G. Phillips and T. G. Madgwick 


@ This book was originally published in Germany during the 
last war under the title DIE HETEROGENEN SCHMELZGLEICHGE- 
WICHTE SILIKATISCHER MEHRSTOFFSYSTEME, and despite war 
conditions sold out an edition of 4500 copies. When Dr. Eitel 
arrived in this country, he arranged to have the book trans- 
lated into English and it is this translation, with the author's 
revisions, which Rutgers University Press has just published. 


Profusely illustrated, its 200 diagrams furnish an abundance 
of material applicable to ceramics, process metallurgy, mineral- 
ogy, and geology. It should also appeal to a large circle of 
workers in high-temperature industries, as well as to students 
in petrology of eruptive rocks. $5.00 


At your bookstore 
RUTGERS UNIVERSITY PRESS 


| 
New Brunswick, N.J. 


THE CARTER 
HYDRAULIC ROCK TRIMMER 


Height over 
all 30 inches 


Maximum 
cutter open- 
ing 7 inches 


Clear width 
between 
columns 
12 inches 


Extra foot (shown) for use when used 
as a portable 
Hydraulic capacity 3 and 5 tons 


Price, 3 ton $100 Price, 5 ton $120 
F.O.B. Cars Buffalo 


For further information write 


THE CARTER MACHINERY SALES CO. 
172 PARKWOOD AVE. KENMORE 17, NEW YORK 


WANTED TO PURCHASE 


Your Sets, Runs and Volumes of 
GEOLOGICAL PERIODICALS 


and other Scientific and Scholarly 
journals in all fields and 
all languages 


Also please send us your list of wants 


ABRAHAMS MAGAZINE SERVICE 
Dept.K 56 East 13th Street’ New York 3, NLY. 


Back Number Magazine Specialists since 1889 


SUBSCRIBE NOW TO THE 
EARTH SCIENCE DIGEST 


The Earth Science Di- 
gest is an international 
illustrated 
issued monthly, 
voted to the Wr 
sciences. It contains 
articles and features of 
lasting interest to the 
professional and ama- 
teur geologist, mineral- 
ogist, and paleontolo- 
gist. Earth Science Ab- 
stracts and the New 
Books column are now 
monthly features. 


SUBSCRIPTION RATES: 
1 year—$3.00; 2 years—$5.00 
(Foreign: 1 vear—$3.50; 2 years—$6.00) 


Single copies—z25¢ each. Vol. 5 began with the 
August 1950 issue. Special subscription rates (10 
or more subscriptions: $2.00 each per year) al- 
lowed to educational institutions and societies. 


THE EARTH SCIENCE DIGEST 


JEROME M. EISENBERG, Editor 
Box G-28, Revere, Massachusetts 


Still Available! 
THE 


FELDSPAR ISSUE 


of the 


JOURNAL OF GEOLOGY 


(Vol. 58, No. 5) 


$0.50 per copy 
Canadian Postage $0.05 
Foreign Postage $0.10 


THE UNIVERSITY OF 
CHICAGO PRESS 
5750 ELLIS AVE. - CHICAGO 37, ILL. 
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Available for fall classes 


PRINCIPLES OF 
GEOLOGY by JAMES GILLULY 


U.S. Geological Survey 
AARON C. WATERS 
Stanford University 


A. O. WOODFORD 
Pomona College 


A glance at the Table of Contents reveals that the 
authors have used an approach emphasizing physical 
principles. A scrutiny of the book itself will show the 
vigor of the treatment and the insistence that the reader 
understand how the basic conclusions of the science 


have been reached. 


CONTENTS 


. Introduction 14. Groundwater 

The Earth's Broad Pattern 15. Deserts 

. Gravity, Isostasy, and Strength 16. Oceans 

Minerals 17. Igneous Activity 

Rocks 18. Earthquakes and the Earth's Interior 
Climates, Weathering, and Soils 19. Mountains 

. Erosion 20. Mineral Resources 


. Geologic Maps, Fossils, and Time Appendixes 


Movements of the Earth's Crust 

. Records of Earth Movement I. Techniques of Topogr aphic Mapping 
of Sell II. International Atomic Weights 

. Rock a III. Identification of Minerals 

. Stream Erosion and Deposition IV. Identification of Rocks 

. Glaciers and Glaciation Index 


wn 


Dr. Robert R. Compton of the Stanford geology staff 
has done all the drawings for this text. 


To be published in June 1951 Probable price $5.75 


W. H. FREEMAN &f% AND COMPANY 
549 Market Street [E)J} San Francisco 5, Calif. 
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NOW AVAILABLE— 
A general index to 


THE JOURNAL OF GEOLOGY 


Volumes XXXVI through LV—1928 through 1947 
Price $3.00; prepaid in U.S.; in Canada add 8¢ for postage; foreign postage 20¢. 


SPECIAL OFFER— 


A general index to the JOURNAL OF GEOLOGY, Volumes 
I through XXX V—1893 through 1927 (regular price $5.00) 
PLUS the 1928 through 1947 index described above, BOTH 
for $5.00. 

Prepaid in U.S.; in Canada add 20¢ for postage; foreign postage 50¢. 


THE UNIVERSITY OF CHICAGO PRESS 
5750 Ellis Avenue + Chicago 37, Illinois 


TO ORDER CLIP AND MAIL THE CONVENIENT COUPON BELOW 


The University of Chicago Press 
5750 Ellis Avenue 
Chicago 37, Illinois 
copies of the JOURNAL OF GEOLOGY INDEX, 1928-1947 


sets of the JOURNAL OF GEOLOGY INDEXES, 1893-1927 
and 1928-1947 @ $5.00 per set. 


(.) Check or money order enclosed. 
(] Please bill me. 
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STRUCTURAL 
GEOLOGY 

CANADIAN 
ORE DEPOSITS 


A Symposium 


This volume, published by The Canadian Institute of pag» 
and Metallurgy in commemoration of its Jubilee Year. . . is 
now available to all. 


It presents the result of investigations by more than 100 
Canadian geologists and mining engineers. 


Its 930 pages describe the main structural features of almost all 
Canadian metalliferous mines and areas. 


Its 930 pages are illustrated by over 300 four-colour reproduc- 
tions, halftones and line drawings. 


Its 930 pages contain 132 separate papers written by more than 
100 mining executives. 


PRICE TO The Special Volumes Committee, 


Canadian Institute of Mining and Metallurgy, 
811, Drummond Building, 
QUALIFIED Montreal, Quebec, Canada 


GENTLEMEN: 


STUDENTS Enclosed is my cheque for $10.00 in payment for one copy of 
“Structural Geology of Canadian Ore Deposits.” Please ship pre- 


’ $8.00 paid to 


Address 
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An A.A.P.G. Publication! 


TECTONIC MAP 


OF THE 


UNITED STATES 


Prepared under the Direction of the Committee on Tectonics 
Division of Geology and Geography, National Research Council 


CHESTER R. LONGWELL, Chairman, PHILIP B. KING, Vice-Chairman 
CHARLES H. BEHRE, WALTER H. BUCHER, EUGENE CALLAGHAN, 
D. F. HEWETT, G. MARSHALL KAY, ELEANORA B. KNOPF, A. I. 
LEVORSEN, T. S. LOVERING, GEORGE R. MANSFIELD, WATSON H. 
MONROE, J. T. PARDEE, RALPH D. REED, GEORGE W. STOSE, W. T. 
THOM, JR., A. C. WATERS, ELDRED D. WILSON, A. O. WOODFORD 


A Geologic Map of the United States and Adjacent Parts 
of Canada and Mexico 


Geologic structure, as evidenced and interpreted by a combination of out- 
cropping areas, bedrock, surface disturbance, and subsurface deformation, 
is indicated by colors, symbols, contours, and descriptive explanation. 
Igneous, metamorphic, and selected areas of sedimentary rock are mapped. 
Salt domes, crypto-volcanic disturbances, and submarine contours are shown. 


The base map shows state boundaries, rivers, a pattern of cities, and | 
degree lines of latitude and longitude. 


The scale is 1:2,500,000, or 1 inch = 40 miles. Printed in 7 colors, on 2 
sheets, each about 40x50 inches. Full map size is about 80x50 inches. 


PRICE, POSTPAID 


$2.00 rolled in mailing tube 
$1.75 folded in manila envelope 
$1.50 in lots of 25, or more, rolled or folded 


THE AMERICAN ASSOCIATION OF PETROLEUM GEOLOGISTS 
Box 979, Tulsa 1, Oklahoma, U.S.A. 


IGNEOUS AND METAMORPHIC PETROLOGY 


By Francts J. Turner and Jean Vernoocen, University of California. 490 pages. $9.00 


Presents a rigorous, unified treatment of the origin and evolution of rocks that have crystallized, 
or have been profoundly modified at high temperatures. The material is based on field and 
petrographic data, and on the physico-chemical behavior of rock-forming minerals as determined 
by experiment and deduced from thermodynamic theory. 


PRINCIPLES OF PETROLEUM GEOLOGY 
By Wriu1aM L. Russet, Texas A. and M. College. Ready in May 


An advanced petroleum geology text, this book provides a discussion of the principles, methods, 
and techniques important in petroleum geology. Emphasizing the practical applications of 
principles, the text covers the origin and accumulation of oil, important structures, classifica- 
tion of fields, stratigraphic features such as reefs, geophysics, and the various well-logging meth- 
ods. Gamma ray and neutron well logs and electrical well logs are treated for the first time. 


NON-METALLIC MINERALS. New 2nd edition. 
Occurrence-Preparation-Utilization 
By Raymonp B. Lapoo, Consulting Engineer, and W. M. Mrers. 605 pages. $10.00 
This text has been completely rewritten in order to give a contemporary picture of present-day 
non-metallic mineral industries. Coverage of the subject is complete, giving an accurate and 
concise picture of the properties of the mineral, where it can be found and how it must be treated 
to be made commercially useful, and in what markets it may be sold. Presentation is technical 
but clear, with emphasis on practical considerations. 


PRACTICAL OIL GEOLOGY. New 6th Edition 


By Dorsey Hacer. 589 pages. $7.50 


The new up-to-date revision of this successful text, the sixth edition covers thoroughly the essen- 
tial facts relating to the occurrence of oil and its extraction, prepared for the driller, operator, 
and mining engineer. Revised material is included on petroleum, stratigraphy, structural 
geology, prospecting and mapping, drilling, natural gas and gasoline, and geophysics. 


Send for copies on approval 


=—— McGRAW-HILL BOOK COMPANY, INC. 
= 330 WEST 42no STREET, NEW YORK 18, N. Y. 
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How the Mining Geologist Works 


MINING GEOLOGY 


By Hugh Exton McKinstry, Harvard University 


Combining all that is common knowledge in the ession with up to-date information 
fom the techni ofthe Bey txt Gecrbes precel the hinds of 
work a geologist must do and how he does it. 

The material is arranged in four parts: 


e Part I—Assembling Geological Data—by mapping and other means of investigation. 
Part Il—Geological of Ore-Search and Ore chiefly the 

touching ages the of ore-search. 
veloping, and valuing ore bodies the 
states from prospecting to opera- 


whether 


Published 1948 680 pages 6”x9” 


Another Popular Geology Text... 


ELEMENTS OF MINERALOGY 


By Alexander N. Winchell, University of Wisconsin 


e Completely modern in content and treatment 
e Presents minerals as natural substances which in many cases vary in composition and 


Gives illustrations of the atomic patterns of the various metals 
e Groups minerals in families 


Makes available at the beginning student’s level a definition of minerals, commonly 


Glossary, index, abbreviations and symbols, 468 charts, diagrams and photos 


Send for today! 


